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Annotation 


In this thesis we discuss how the brane world scenario can be realized dynamically 
within the field theoretical framework using topological solitons. As a playground 
we consider a bosonic sector of a (4+1)-dimensional supersymmetric gauge theory, 
which naturally supports soliton of co-dimension one, a domain wall. We first 
discuss separate localization of matter fields and gauge fields on the world-volume 
of the domain wall and then we present two explicit five-dimensional models, where 
both matter fields and gauge fields are localized together with minimal interactions. 
We show that matter fields localize in the adjoint representation of the non-Abelian 
gauge group and we calculate the effective interaction Lagrangian of these matter 
fields up to the second order in derivatives. We discuss similarities of our models 
with effective models describing pions in QCD and with D-branes from string 
theory. 

Key words: Brane world scenario, topological soliton, supersymmetry. 


A not ace 


V teto disertacni praci se zabyvame otazkou dynamicke realizace scenare branovych 
svetu pomocf topologickych solitonu. Nejprve diskutujeme lokalizaci hmotovych a 
kalibracnich poll na domenove stene oddelene a potorn uvazujeme explictni modely 
ve (4+l)-dimenzich, ve kterych se hmotova a kalibracni pole lokalizuji soucasne s 
minimalm interakcf. Pote vypoci'tame efektivni interakcnf Lagrangian hmotovych 
poll' do clenu druheho fadu v derivaci'ch. Nakonec diskutujeme podobnosti nasich 
modelu s efektivnimi modely pionu v QCD a s D-branami ve strunove teorii. 


Klfcova slova: Scenar branovych svetu, topologicky soliton, supersyme- 
trie. 
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Preface 


In this thesis we discuss how certain aspects of the so-called brane world sce¬ 
nario can be realized within the field-theoretical framework. In accordance 
with an old proposition, and for simplicity as well, we use a domain wall as a 
model of the brane. Thus, we consider the case of single, infinitely long extra 
dimension, along which domain wall varies. We show that, in the presence 
of the brane, certain modes of fluctuations of background fields are trapped 
inside the domain wall, representing particles localized on the brane. Then, 
according to specified rules, we derive an effective four-dimensional theory 
describing nonlinear interactions of these modes. The ultimate goal of our 
work is to localize a full spectrum of Standard model particles in this way. 
Here, we present an important step towards this goal, a model where var¬ 
ious scalar particles (Nambu-Goldstone bosons) with non-trivial quantum 
charges are localized on the brane along side with non-Abelian gauge fields. 
We show that the corresponding effective theory shares many aspects with 
low-energy theory of pions. In addition, we explore the full structure of non¬ 
linear interactions between pions up to second order terms in derivatives. 
We also show that our model naturally incorporates the Higgs mechanism 
which is, in contrast to the Standard Model, based purely on geometry. 
Lastly we address the issue of stability and typicality of our construction. 

This thesis is organized as follows. In Ch. [I] we make short historical 
overview of extra dimensions and their phenomenological implications for 
high energy physics. Ch. [2] has an introductory purpose too, as we consider 
there a simple toy model realizing the brane world scenario. Along the 
presentation we encounter several important concepts, which we utilize in 
the later chapters, such as a BPS soliton, zero modes and their localization 
and description by an effective Lagrangian. In Ch. [3] we make a detour and 


discuss another key concept, supersymmetry. First we establish connection 
between BPS solitons with supersymmetry in (l+l)-dimensions as a natural 
extension of our toy model of Ch. [2} Then we discuss basic aspects of super- 
symmetric theories in four and subsequently in five space-time dimensions. 
As we will demonstrate, supersymmetric gauge theories in five-dimensions 
are very convenient laboratories to explore the brane world scenario based 
on domain walls. The culmination of this discussion is represented by chap¬ 
ter [4j where we address the crucial issue of localization of gauge fields. At 
the end of this chapter all tools will be prepared to make a next logical step 
of Ch. [5j with which we start to review original results. 

In Ch.[5]we introduce simple, yet robust field-theoretical model of the brane 
world scenario. We demonstrate how the matter fields and gauge fields can 
be localized on a domain wall with nontrivial interactions. This chapter is 
based on the paper I listed below. In Ch.[6]we present a similar model to that 
of Ch. [5} but with more economical charge assignments. This subtle change 
in turn ensure questionable stability of the gauge sector of the previous 
model and also results in new features in the effective Lagrangian. This 
chapter is based on the paper II listed below. Finally we conclude in Ch. [7] 
by discussing possible extensions of our work. 

The author of this thesis took an active part in this research and main 
aspects of his contribution were twofold. First, he found closed formulas for 
low energy effective Lagrangians presented in Ch. [5] and Ch. [6j And second, 
he investigated the mechanism, which removed the potential instability of 
the gauge sector and which represents the main part of Ch. [6j 

For the reference we list below all publications to which the author of this 
thesis contributed during his PhD study. 
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M. Arai, F. Blaschke, M. Eto, and N. Sakai. Matter fields 

and non-Abelian gauge fields localized on walls. Progress of 

Theoretical and Experimental Physics , 2013(1):010003, January 2013. 
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matter and gauge fields localized on walls. ArXiv e-prints , 
March 2013. 
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Masato Arai, Filip Blaschke, Minoru Eto, and Norisuke 
Sakai. Localization of matter fields and non-Abelian gauge 
fields on domain walls. Journal of Physics: Conference Series , 
411(1):012001, 2013. 
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Introduction 


The Standard Model (SM), developed in early 1970s, successfully describes three of the 
four fundamental interactions (namely the strong, weak and electromagnetic interac¬ 
tions) of elementary particles on energy scales accessible to particle accelerators such as 
the Large Electron Positron Collider, the Large Hadron Collider at CERN in Switzer¬ 
land and the Tevatron at Fermilab in the United States. The SM predicts phenomena 
observed in these particle accelerators quite well and actually there is no inconsistency 
between the SM and current experimental data. However, it is widely recognized that 
SM is not a complete theory by a number of theoretical and phenomenological reasons. 
First of all, the SM does not provide answers for questions such as why the spectrum of 
particle masses is hierarchical or why the number of generations is three. It also does 
not explain the relic abundance of dark matter in the universe, as none of the particles 
in the SM can be a candidate for dark matter. Furthermore, the SM has a conceptual 
problem called the gauge hierarchy problem, stating that there is a huge hierarchy (17 
orders of magnitude) between the electroweak scale and the gravity (Planck) scale. Fi¬ 
nally, the SM does not describe gravity at all, as the correct quantum theory of gravity 
is still missing. The common expectation is that there is a more fundamental theory 
beyond the SM, in which these problems are naturally resolved. 

A very popular and one of the most promising ideas beyond the SM is an intriguing 
notion of extra dimensions. The possibility, that there might be more dimensions of 
space than three, has become a part of serious science more than a hundred years ago 
and since that time it influenced physics in many direct and indirect ways. Conse¬ 
quently this subject, that one could call “extra-dimensional physics”, is very vast and 


1 



1. INTRODUCTION 


progressively more richer in content as we follow its development up to the present 
day. Therefore, given limited space of this chapter, our presentation of this topic is 
accordingly coarse-grained. In fact, from the plethora of different models/theories, 
based on extra dimensions, we will focus only on the four most influential ones. Each 
of these corresponds to a shift in thinking about the role of extra dimensions in high- 
energy physics and also marks the beginning of a fruitful and ongoing investigation. 
These four milestones reflect in the organization of this chapter. Its goal is to provide 
enough background material to properly set the remaining content of this thesis in the 
landscape of extra-dimensional scenarios. 

The beginning of extra dimensions, which we describe in the first section, started 
from attempts to unify different forces in nature. Even though first investigations in this 
direction were made by Nordstrom, we are going to introduce the more famous work of 
Kaluza and Klein BE!- As we will see, in the Kaluza-Klein scenario the mathematical 
consequences of additional, fifth dimension result in remarkable unification of four¬ 
dimensional gravity and electromagnetism. However, the presence of the fifth dimension 
is considered here as an unwanted feature and, in order not to alter phenomenology, it is 
compactihed to an unobservable size. We will also very shortly describe developments 
of Kaluza-Klein ideas within string theory. 

In the second section we discuss a proposition of Rubakov and Shaposhnikov Gam, 
in which the role of an extra dimension is completely different. The central idea of 
their work, which later becomes identified as a brane world scenario, is based on an 
assumption, that our Universe is a three-dimensional defect evolving in a multidimen¬ 
sional space, called a bulk. All matters are trapped on this defect and, if the energy 
is not sufficiently high, it cannot propagate into other dimensions. Rubakov and Sha¬ 
poshnikov in their famous paper [4J discussed a simple model, in which this scenario is 
realized. We will revisit this model in great detail in chapter [2] In fact, the work we 
are going to present in this thesis, might be considered as an extension of this model 
towards to its logical conclusion. 

The canonical form of the brane world scenario, however, was fully realized later in 
the work of Arkadi-Hamed, Dimopoulos and Dvali J5J, which is described in the third 
section. Their model, known as the ADD model, rather than hide additional dimensions 
under the carpet, used the bulk as an instrument to explain phenomenological truths 
about the nature. In ADD scenario the extra dimensions are compactihed, but in 
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contrast to the Kaluza-Klein scenario, not to the Planck scale Mp 10 19 GeV, but 
rather to experimentally more accessible Weak scale M\y ~ 10 3 GeV. As a consequence, 
they are called large extra dimensions , even though from human perspective they are 
still rather small (around micrometer or less). The reason for this is to account for 
apparent weakness of gravity. In the brane world scenario, gravitons are not, like 
ordinary matter, trapped on the brane, but they are free to roam in the bulk. This 
results in weakening of gravitation from a point of view of a brane-bound matter. We 
also briefly describe different approach to the brane world scenario, which is the work of 
Randall and Sundrum [SE]. To localize particles on the brane, they use warped extra¬ 
dimension, compactified on a S 1 / Z 2 orbifold and a negative cosmological constant in 
the bulk. In contrast to ADD scenario, the fundamental scale in models of Randall 
and Sundrum is the Planck scale. However, the warping of the extra dimension causes 
exponential screening of the fundamental scale down to My/. 

The dichotomy between gravity and other forces in the brane world scenario ac¬ 
tually originated in string theory. Indeed, the very word “brane” comes from string 
theory. In 1995 Polchinski [8] showed that string theory possesses a remarkable degree 
of freedom: a D-brane. The D-brane naturally realizes a brane world scenario, since 
open strings (matter and gauge particles) ends on the D-brane. On the other hand, 
closed strings (gravitons) cannot be confined to the D-brane and must propagate into 
extra dimensions. Introduction of D-branes led to new ideas in phenomenology such 
as the above mentioned large extra dimensions. Moreover, it is widely recognized that 
topological soliton such as a domain wall, vortex and instantons can be seen as an 
effective low-energy analog of a D-brane. This opens up the possibility to realize brane 
world scenario dynamically and to directly investigate the spectrum of particles local¬ 
ized on the brane in the field-theoretical framework. We discuss this issue in the four 
section, where we also set the goals of this thesis. 

1.1 Kaluza-Klein theory 

In 1919 Theodor Kaluza came up with an elegant way how to put gravity and electro¬ 
magnetism into one field-theoretical frame. Considering that at that time those were 
only known forces, there is no wonder that his theory raised serious attention in the 
scientific community. Indeed, upon submitting his paper to Einstein himself, Kaluza 
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had to wait two years (an unusually long time even by the standards of those slow-paced 
days) to receive a confirmatory response and recommendation for publication [I]. 

In his work he noticed, that if one starts with pure five-dimensional gravity one can 
reduce it to Einstein-Maxwell theory in four dimensions. Curiously, the only ad hoc 
assumption he had to make, was that no fields depends on the fifth direction (which 
we from now on denote as y). He discovered that this assumption alone naturally 
decomposes five-dimensional metric tensor into three four-dimensional fields of spin 
2,1 and 0 and, upon plugging the appropriate ansatz into the five-dimensional Einstein 
equations, one obtains separate sets of equations for each field0 which are precisely the 
correct equations of motion for each spin. It worked like a miracle. 

Let us briefly expose the mechanism of this miracle from a modern point of view. 
The spectrum of fields in effective four dimensional theory is easily shown to be degrees 
of freedom (DoF) consistent. Indeed, we started with five-dimensional metric tensor 
5 gMN with M, N = 0,1, 2,3, y which contains 5(5 + l)/2 = 15 real DoF. With respect 
to y we can split 5 gMN into the four-dimensional metric tensor g yu = with 10 DoF 
and a vector field A y = 5 c/ w of altogether four DoF. The last degree of freedom is 
associated with the scalar field <j) = 5 g yy , which sums up back to 15, as required. 

The spectrum itself, however, can be deduced from symmetry considerations. The 
independence of metric tensor 5 gMN on extra-dimensional coordinate y is formulated 
by the so-called Cylinder condition d y 5 gMN = 0, which is obviously not a covariant 
statement. However, there are two kinds of transformations which are compatible with 
this condition, namely 

I type: x M -A f^(x)df? + yd™ , (1.1) 

II type: x M -» x^6ff + (y + e(x))5™ . (1.2) 

The first one represents a general change of coordinates in four dimensions with no 
change to extra dimension (for simplicity we ignore here global transformations such 
as rescaling and constant shifts). Under this type of transformations four-dimensional 
fields g yu , A^ and <fi transform respectively as a rank-two tensor, vector and a scalar. 
The second type can be respected as a gauge transformation. For infinitesimal shifts 

1 Apart from scalar field which is coupled with both gravitational and gauge kinetic terms. 
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|e(x)| <C 1 we get 


9[iv t Q[iv Q^eA^ , 

(1.3) 

y A y d^scj) , 

(1.4) 

4>- 

(1.5) 

This leads to introduction of the following quantities 


9[iv — 9/iv — A y A v /cj) , A ^ — A [i/(j ), <f> — cj). 

(1.6) 


As one can easily check, these have correct transformation properties under the gauge 
transformation and enable us to identify them as a four-dimensional metric tensor, 
U( 1) gauge field and a real scalar field. 

At the time of Kaluza, the scalar field <j> was an unwanted feature, which was 
forced out of the picture. In the present time, however, it bears the name radion (or 
graviscalar or also dilaton ) and it has a significance to modifications of Einstein theory, 
supergravity and string theories. 

The Cylinder condition of Kaluza d y 5 gMN = 0 was a main source of criticism of 
his otherwise stunning discovery. Another obvious objection was the manifest lack 
of the fifth direction one supposedly should be able to see, if one takes the Kaluza 
idea literally. Kaluza himself, however, did not believe in actual existence of this fifth 
dimension, but considered it as a useful mathematical construction. Such a view, 
however, was changed few years later when Swedish physicist Oskar Klein [2] presented 
an idea, that the invisibility of the fifth dimension is due to the fact that it is compact , 
more precisely a circle of extremely small radius. This construction not only settled 
most of the qualms but also gave an unexpected bonus: a quantization of charge. 

If the fifth coordinate is periodic, the Cylinder condition is not necessary. To obtain 
a desired four-dimensional effective theory, one decomposes all fields into a series of its 
Fourier modes: 


OO 


9nv(x,y)= 9${x)e lky/R , 

k=— oo 

(1.7) 

OO 

A,(x,y)= Y AP(x)e- iky/R , 

k=— oo 

(1.8) 

oo 

<t>{x,y)= Y 0 (fc) (*)e“ lfc ^ /jR , 

(1.9) 


k=—oo 
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where R is a radius of the circle. 

It can be seen most easily from five-dimensional Klein-Gordon equation 

OO 

(d 2 - d 2 )(p(x,y) = 0 ^ (d 2 ~ m 2 k )^ k) {a:) = 0 (1.10) 

k =—oo 

that all modes with k ^ 0 are massive with the mass m|. = k 2 / R 2 . Only zero modes 
k = 0 remains massless. These zero modes are then identified with the known four¬ 
dimensional fields our universe is filled with (i.e. gravitons and photons). The presence 
of an extra dimension is only visible through the infinite Kaluza-Klein (KK) tower of 
massive modes. In addition, if the radius R is sufficiently small, even the lightest of 
massive modes will be too heavy to affect the physics at our energy scales and thus the 
whole theory is effectively four-dimensional. 

Let us also point out, that spectrum of massive modes consists only of massive 
gravitons and not of massive photons nor scalars. The reason is that while massless 
graviton has only two polarizations, massive one has five. Similarly, massive gauge 
field must have three polarizations. In order to accommodate these extra degrees of 
freedom, a graviton in the process very alike the Higgs mechanism eats a scalar field 
and a gauge field to form a proper massive mode with five physical degrees of freedom. 

The quantization of the charge is a direct consequence of a periodic nature of the 
extra dimension. The charge is identified with the momentum along the y direction, 
which is a conserved quantity. The compact nature of this dimension then forces the 
momentum to be quantised as multiples of \ZSttG 2 /R. Upon identification of this 
with elementary charge e one can estimate the radius of the extra dimension to be 
R ~ 1CU 33 cm. In other words, the Klein-Gordon picture demands that characteristic 
length scale of extra dimension is close to the Planck length. If this is correct, there 
is really little hope that we would be ever able to probe such tiny distances by our 
experimental devices directly. 

1.1.1 Compactification and string theory 

Despite the lack of direct phenomenological consequences of the Planck length sized 
extra dimensions, the idea was pursued further, mainly to explore its indirect conse¬ 
quences. In the 1920’s it was Schrodinger, Gordon and Fock, while in the 1930’s it 
was Einstein [9j, Mandel m and Pauli HU. Most notably in 1953, a year before the 
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introduction of non-Abelian gauge symmetries into the quantum field theory by Yang 
and Mills, Pauli discovered all its essential features in the work based on Kaluza-Klein 
ideas where he compactified a sphere rather then a circle m • However, since he could 
not avoid having a massless vector meson in the spectrum, which was not observed 
experimentally, he decided not to publish his results 0 

After this, there was a long pause before the compact extra dimensions were revived 
in the string theory framework. It started with the work of Scherk and Schwarz in 
1972 [13]. They proposed a superstring model based on a product space of our four¬ 
dimensional space-time with six-dimensional compact manifold. This observation lead 
to the discovery of Candelas [l4j that, if the compact manifold is the so-called Calabi- 
Yau manifold, in the low energy limit one recovers E( 8) x E( 8) gauge group which 
contains the Standard Model gauge symmetry and three generations of chiral fermions. 
Actually, this discovery (published in 1985) triggered what is usually referred to as the 
superstring revolution. Since the typical size of a Calabi-Yau manifold in superstring 
theory is the Planck length, this picture can be regarded as a culmination of the Kaluza- 
Klein ideas. 

Systematic research of superstring-based models of the SM type then started in 
1990’s. With the turn of the century, however, the D-brane concept |8j dramatically 
increased plausibility of such models naina. In these models, intersecting D-branes in 
low energy limit naturally provides the right structure of gauge interactions and give 
the correct number of quark/lepton families such that they are anomaly-free. It also 
gives Dirac neutrino masses and stable proton quite naturally. This development based 
on superstring/D-brane engineering, although very promising, does leave zero room for 
experimental observations of extra dimensions, being all confined to the Planck size. 

1.2 Rubakov &; Shaposhnikov model 

In 1983 two Russian physicists V. A. Rubakov and M. E. Shaposnikov asked a simple 
question, which also makes the title of their famous paper: “Do we live inside a domain 

2 Today, we base our model of electro-weak interaction on almost the same gauge group as in Pauli’s 
work. However, massless vector mesons (gauge bosons in modern terminology) are not present in the 
spectrum due to the Higgs mechanism, where the symmetry is spontaneously broken and these particles 
become massive. 
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wall?” Their three-pages long text [3j marked the beginning of new thinking about 
extra dimensions, nowadays commonly referred to as the brane world scenario^ 

Their idea is based on the following observation. Assuming there are additional 
dimensions than three, the phenomenology dictates that these extra dimensions must 
be hidden. One possible way how to do this is a compactification as in the Kaluza-Klein 
theory. Another possibility, as Rubakov and Shaposhnikov argue, relies on dynamics. 
In their scenario all particles are trapped inside a potential well, sufficiently narrow 
in the direction of extra dimensions and flat in ordinary space directions. This means 
that propagation in extra dimensions is energetically unfavored resulting in effective 
three-dimensional appearance of our Universe. Notice, however, that particles are not 
forbidden to climb out of the well, provided they have enough energy. Such an event 
would seem to violate energy and momentum conservations from a point of view of 
observes inside the well. 

The origin of the well-like potential can be purely dynamical. In their work [2] 
Rubakov and Shaposhnikov illustrate this fact in a toy model, where the potential is 
realized by a background solitonic solution: a domain wall. As we are going to give all 
the details in Ch. [2j let us here only touch the most important points. 

Consider a five-dimensional theory defined by an action 


S 


d 5 xL 



U(^) , 


( 1 . 11 ) 


where M = 0,1, 2, 3, y and where 

V(4>) = . 


( 1 . 12 ) 


Notice that this model has two degenerate vacua (j) = ±v. A domain wall solution is 


given by 


4>k{v) = v tanh 



(1.13) 


3 It turned out, as so many times before, that Rubakov and Shaposhnikov were not the first one to 
propose a brane world scenario. In 1982 Japanese physicist K. Akama published a paper El, where he 
discuss a possibility that our world is an object called a vortex, a soliton of co-dimension 2, immersed 
in a six dimensional flat space-time. He then argues, using rather involved path-integral methods, that 
small fluctuations of the vortex background induce an Einsteinian gravity on the vortex’s world-volume. 
Although the approach is very different, the idea exactly match what we now perceive as brane world 
scenario, with its huge phenomenological implications. This paper, however, was not widely known 
and so its origin is consistently connected with a year older paper of Rubakov and Shaposhnikov. 
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1.2 Rubakov &; Shaposhnikov model 


The non-trivial topology of this solution can be seen from the fact, that 4>k{v) ap¬ 
proaches different vacuum value ±v as y —> ±oo. Since both vacua are true vacua, it 
costs no energy for a field to be near them. Therefore almost all energy of a domain 
wall is stored near the point dq, where the transition between vacua takes place. Since 
most of the energy is localized at the wall’s position, it does act as a barrier, separating 
the space into two pieces. 

Let us investigate small fluctuations around this solution. First we denote 0(x^, y) = 
(/>A'(y)+e(a: /i , y), where y, = 0,1, 2, 3 and then we plug this into corresponding equations 
of motion. Retaining only linear terms in e we have 


d\e — 



(1.14) 


To find normal modes we employ the ansatz e = 6(y)e _lkx+lajt . Putting this back into 
the above equation we obtain a Schrodinger-like eigen-value problem for the “wave 
function” b(y) 

~dyb(y) + U(y)b(y) = {u 2 - k 2 )b(y ), (1.15) 

where U(y) is a potential given as 


U{y) 


2s\uh 2 (vJ\{y - y 0 )) -1 
Xv 2 -- 7 =- — - 


cosh 2 (vy | (y - y 0 )J 


(1.16) 


This is the well-like potential Rubakov and Shaposhnikov had in mind (see Fig. 2.2). 
In Ch. [2] we provide explicit formulas for the entire spectrum of eigenmodes of this 
potential. We will also argue that a few lowest modes are trapped inside the well. 
These are identified as particles localized on the brane. 

The above analysis, however, considers only scalar fields. Rubakov and Shaposh¬ 
nikov also showed that fermionic field can be localized as well. One simply adds fol¬ 
lowing terms to the model (|1.11[) 




(1.17) 


which describes a Dirac spinor coupled to a scalar field cj) via Yukawa term with the 
coupling constant h. The equation of motion for fermionic fluctuations y) is given 
as 

i$i/)(x,y) + h0 K (y)ip(x,y) = 0, (1.18) 
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with the solution 

i/>(x,y) =e-^ dy, ^V 0) (*), (1-19) 

where ^°\x) is a (3+l)-dimensional left handed massless spinor. Moreover, the profile 
^(x, y ) along the extra dimension shows that ip(°\x) is, indeed, localized on the domain 
wall. 

Obviously, the toy model of Rubakov and Shaposhnikov cannot be considered as 
a realistic realization of the brane world scenario, as there is only single scalar and 
single fermion localized on a domain wall. One can, however, continue to improve the 
methods how to simultaneously localize more fields and how to make them interact. 
Indeed, developing such a “domain-wall engineering” is the aim of this thesis. The 
true value of the work of Rubakov and Shaposhnikov lies in the “existence proof” 
that the brane world scenario can be realized within the framework of (quantum) field 
theory. In other words, they dispelled the Kaluza-Klein paradigm and opened a new 
line of investigation of extra-dimensional physics. It took, however, another 15 years 
before people realized the full power of the brane world scenario with respect to the 
SM phenomenology. We will describe this shift in thinking in the next section. 

1.3 Brane world scenario 

One of the yet unsolved problems of the SM is the gauge hierarchy problem. Simply 
speaking, this problem points to the stupefying mismatch between what we think should 
be a fundamental scale of the SM and what it actually is. The fundamental scale of 
the SM should be the Planck scale. Recall that the Planck scale is a unit independent 
energy scale obtained as a combination of fundamental constants of physics, namely the 
gravitational constant, speed of light and the Planck constant Mp = 10 19 GeV. 

This is the energy at which effects of quantum gravity cannot be neglected and where 
the SM breaks down. This statement is supported from analysis of running of coupling 
constants of all forces, which shows that they meet (roughly) at Mp. The fact that the 
natural size of Higgs dimensionful parameter, which controls masses of W and Z bosons 
and fermion masses, is around the Weak scale M\,y ~ 10 3 GeV, 17 order of magnitudes 
lower than our expectation is the crux of the gauge hierarchy problem. 


10 



1.3 Brane world scenario 


There is, however, a deep assumption at work here. While non-gravitational forces 
has been experimentally probed reasonably close to the Weak scale, precision measure¬ 
ment of gravity is, at best, approaching 1 MeV or a micrometer distances m- Can 
we really be confident that our low-energy concepts of how gravity works holds down 
to the weak scale distances 1 TeV ~ 10 -19 m or even down to the very Planck length 
1(T 35 m? It might be the case that the Planck scale is only a low-energy illusion and 
that the true fundamental scale of gravity is much lower. 

This possibility is realized within the brane world scenario quite naturally. There 
are two typical models of the brane world scenario. One of them is proposed by Arkani- 
Hamed, Dimopoulos and Dvali (ADD) [19] . The mechanism in this model is the follow¬ 
ing. One can argue that gravity, connected with the geometry of space-time, cannot be 
confined to the brane and must “leak out” into the bulk. This will lower the apparent 
strength of gravitational force from a point of view of brane-bound observers, setting 
the effective four-dimensional Planck scale Mp^ much higher than what the actual 
(4 + n)-dimensional Planck scale Mp +n ' > is. This “diluting effect” of n compact extra 
dimensions can be easily quantified. Let us assume that typical size of extra dimensions 
is R. Then, for distances r<Cfi, the effective four-dimensional gravitational force may 
be written either in terms of effective Planck scale Mp ^ or in terms of actual Planck 
scale Mp 4+n ^ as 

F(r ) ~ 

This is giving us the relation 


mi?7i2 1 
(. M i p ) ) 2r2 


m\vri2 


(. M ( p +n) ) 2 + n R n r2 


( 1 . 20 ) 


(Mp ] ) 2 


(Mj 


(4+n) \2+nj^n 


( 1 . 21 ) 


If we set Mp ^ 10 19 GeV, then the radius of extra dimensions must be of the order 


R 


— _i q 

10 n 9 m X 


1 TeV 


M 


(4+n) 


1+2 

n 


( 1 . 22 ) 


In the ADD model, most interesting case of Mp + "' 1 = Mw is discussed, which 
simplifies the above formula to 


R 


ln ——19 
~ 10 n m . 


(1.23) 


For n = 1 we obtain R ~ 10 11 m implying deviations from Newtonian gravity over 
Solar system distances, which is clearly unacceptable. For n = 2, however, we have 
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R ~ 1CP 4 m, which is experimentally accessible by today’s instruments, promising an 
exiting new physics just behind the current experimental boundaries. Compared to 
Kaluza-Klein models, where R is typically around the Planck length, we see that in 
ADD scenario, R is quite large. For that reason, such a medium scale extra dimensions 
are called large extra dimensions (LED), even though from a human point of view, they 
are still quite small. 

In the work [19] the authors presented a particular realization of n = 2 scenario, by 
considering a six-dimensional theory with our Universe localized at the core of a vortex. 
We are not going to describe the details of this model, since the phenomenological 
consequences of LED idea are not sensitive to them. In fact, the same trio of authors 
later discussed the phenomenology in the paper [20] and, together with Antoniadis [2T], 
they outlined a way of embedding the ADD model in superstring theory. 

The ADD scenario, which can be seen as peculiar marriage of Kaluza-Klein ideas 
with the concept of brane-worlds, sparked a paradigm shift in thinking about the extra 
dimensions, known as the LED paradigm (a term, which we borrow from [22]). After 
ADD many other studies followed, giving rise to a fruitful new direction in the high- 
energy physics, which continues to the present day. The lasting popularity of LED- 
based theories is a result of a wish to have new physics at the Weak scale. Moreover, 
LED paradigm opens up possibilities to explain some of the properties of SM in a 
new way. Let us, as an example, mention the work of Dvali and Shifman [5], where a 
hierarchy of fermion masses is explained without invoking broken family symmetries. 
Instead, Dvali and Shifman consider a case, where each generation lives on a separate 
brane with gauge fields propagating freely in the bulk. The Higgs fields are localized on 
another brane with an exponentially decaying profile in the bulk. In this way the closest 
brane has much larger overlap than the brane further away and even larger than the 
furthest one. Since the overlap of profiles in the bulk determines the interaction strength 
of fermions with Higgs field, this arrangement leads to a large hierarchy between masses 
in each generations, even though there is no hierarchy of distances between the branes. 

Many old problems were revisited within the LED approach. Let us, again, pick 
up a few examples. The issue of proton stability, for one, must have been completely 
rethought. The reason is that with strong gravity at Mw scale, any global symmetry, 
such as baryon number, is violated by virtual black hole and wormhole production, 
leading to a rapid proton decay in LED-based theories. Many clever ideas how to 
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circumvent this catastrophe were invented. In the original ADD model, the authors 
resolved this by invoking some new physics above M\y scale, which suppresses these 
dangerous processes. More concrete possibility is based on an idea to turn the baryon 
number conservation into a discrete gauge symmetry [23], which is protected by the 
gauge principle. The phenomenological implications of this method was shown to be 
compatible with experimental data (see |23J for a review). More geometrical approach 
to stabilize the proton was put forward by Arkani-Hamed and Schmaltz (25] . Their idea 
relies again on exponentially decaying profiles of localized particles on a brane. They 
considered a case of a thick wall with quarks and leptons localized on different ends. 
Because of this, interaction terms between quarks and leptons, which are proportional to 
the overlap of their profiles, are exponentially suppressed, leading to the experimentally 
save lifetime of the proton. 

The conventional approach to the issue of neutrino masses is based on the so-called 
seesaw mechanism, which relies on the enormity of the fundamental scale of the SM. 
This is an example of another problem, conventionally regarded as “solved”, which 
requires a reconsideration within LED paradigm. Its possible solution was actually 
proposed by the authors of ADD scenario themselves |[26). Their idea was to introduce 
a right-handed neutrino, freely propagating in the bulk. The neutrino mass term is 
proportional to the overlap of the left-handed neutrino, confined to the brane, and the 
right handed neutrino (see [22j for detailed formula). Since the right-handed neutrino 
propagates in the entire bulk, its wave function is suppressed by a factor V n 7 , where 
V n is the volume of extra dimensional space. This effect places the neutrino mass in 
the right ballpark. Thus, from a point of view of LED-based theories, light neutrinos 
and the weakness of gravity are both consequence of the same phenomenon: large extra 
dimensions. 

The other model was proposed by Randall and Sundrum Ed, which is called a 
warped extra-dimensional model. This is a five-dimensional model, where one extra¬ 
dimension is compactified on a S' 1 /Z 2 orbifold and a negative cosmological constant is 
introduced in the bulk. Two D3-branes are placed at fixed points of the orbifold (f> = 0 
and <j) = 7T (cj) is an angle of S 1 ) with opposite brane tensions. A brane at r/> = 0 with a 
positive tension is called the hidden brane and the other one at <f> = n with a negative 
tension is called the visible brane on which the SM fields localize. Solving the Einstein 
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equation of this system, the five-dimensional bulk geometry is found to be a slice of 
anti-de Sitter (AdS) space, 

ds 2 = e~ 2nrc ^ rj^dx^dx 1 ' - rldcj) 2 , rj^ = diag(l,-1,-1,-1), (1.24) 

where k is the AdS curvature in five dimensions and r c is a compactification radius. This 
background geometry allows us to take the Planck scale as a fundamental scale. Indeed, 
in the effective four-dimensional description an effective mass scale on the visible brane 

— /CjN _ 

is warped down to A n = Mp e~ nKrc due to the effect of the warped geometry, where 
Mp is the reduced Planck scale in the five-dimensional space-time. Therefore, with a 
mild parameter tuning, nr c — 12, we can realize A^ = 0(1 TeV) and obtain a natural 
solution to the gauge hierarchy problem. 

Phenomenological aspects of the warped extra-dinrensional model have been vigor¬ 
ously investigated, as was performed in the ADD model. Apart from that, the warped 
extra-dinrensional model has an interesting connection to the AdS/CFT correspondence 
mmm- This correspondence tells us that four-dimensional theories with strong 
interactions is related to weak-interacting five-dimensional supergravity theories. In 
other words, one can understand the properties of five-dimensional fields as those of 
four-dimensional composite state m- This has opened up new directions for tackling 
questions in particle physics such as the flavor problem, grand unification and the origin 
of the electro-weak symmetry breaking or supersymmetry breaking. 

1.4 Purpose of this thesis 

The brane world scenario explained briefly above has a lot of interesting theoretical 
and phenomenological features, but there is a big assumption that the SM fields are 
localized on the three-dimensional hyper-surface, a brane. As mentioned above, such 
a situation is naturally realized in the string theory. One of the promising candidates 
of the fundamental theory is the superstring theory, where the number of space-time 
dimensions is ten. Its low-energy effective theory is the so-called supergravity that is 
described as the ten-dinrensional field theory. If the SM is realized as a low energy 
effective theory of the supergravity/the string theory, it is necessary to make a gap of 
the dimensions between the SM and the supergravity. One of ways is to compactify 
the space-time dimensions, but it is usually performed by hand and therefore it is not 
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dynamical. The other way is to use a topological soliton, where fields naturally localize 


around that, as explained in subsection 1.2 


However, there is an obstacle that the mechanism of Rubakov and Shaposhnikov 
works only for bosons and fermions. In order to realize the SM model on a topological 
solitons, it is necessary to localize not only bosons and fermions, but also massless gauge 
fields associated with the SM gauge groups SU (3) X SU (2) x[/(l). Naive application of 
the mechanism by Rubakov and Shaposhnikov to bulk gauge fields tells us that gauge 
fields localize around a topological solution, but inevitably become massive, leading 
to the violation of the four-dimensional Coulomb’s law (details will be discussed in 
Ch. [4]). It has been a long-standing problem in construction of a brane world scenario 
in the field theory framework. However, recently, a very simple mechanism to localize 
massless gauge fields has been proposed by Ohta and Sakai [31]. The key ingredient 
of this mechanism is a position-dependent gauge coupling that is a coefficient of the 
gauge kinetic term. Consider the five-dimensional space-time. If the position-dependent 
gauge coupling has a profile such that it goes zero at infinities of extra dimension and 
has a peak around a topological soliton, it can be shown that the massless gauge held 
localizes and that the four-dimensional Coulomb’s law holds. It is a great step towards 
a construction of a realistic brane world scenario. 

For a realistic model building, we further need to proceed along this line. In this 
thesis, we will propose two five-dimensional models where gauge fields and charged 
matters localize on a domain wall solution. The first model, introduced in Ch. [5j 
successfully makes matters and gauge fields localize, but it has an instability that 
the position-dependent gauge coupling can be negative. The second model of Ch. [6] 
improves this shortcomings. We will also derive the low energy effective actions on a 
domain wall, based on a technique called the moduli approximation |32] (we discuss 
it in Ch. [4]). Resultant actions are similar to the chiral Lagrangian being the effective 
theory of QCD. 
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2 

Brane world scenario and solitons 


The key ingredient how to realize the brane world scenario in the field-theoretical 
framework is a topological soliton. In order to illustrate this fact, in this chapter we 
discuss in detail a simple model containing a 3-brane. This is a (4+l)-dimensional 
theory with single scalar field and a quartic potential with degenerate vacua. Non¬ 
trivial vacuum structure allows us to construct a solitonic solution, varying only in 
one dimension, which we identify as a brane. Then, studying small fluctuations on 
the background of this solution, we discover that in the direction of extra-dimensional 
coordinate the effective potential has a well-like shape, centered at the position of the 
brane. Inevitably, some modes of the fluctuations will be trapped in this well and we 
speak of them as of particles localized on the brane. The culmination of our discussion 
will be the construction of the effective Lagrangian describing dynamics of these trapped 
modes on (3+l)-dimensional world-volume of the brane. 


2.1 A brane-like solution 


Let us consider a theory defined by an action 


S 


d 5 x£ 



7,d M 4>d M (l) — V((f>) , 


( 2 . 1 ) 


where </> is a scalar field and M = 0,1,2,3, y is a coordinate index of (4+l)-dimensional 
spacetime with mostly negative signature In such theories, the question of 

existence and properties of solitonic solutions is determined by the topological structure 
of the so-called vacuum manifold. The vacuum manifold is simply a set containing all 
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homogeneous (i.e. space-time independent) solutions to a equation of motion (EoM) 
with zero energy E = 00 where: 


E = 



d 4 x 


^ 1 2 + ^V0 2 + E(0) 


( 2 . 2 ) 


These solutions are called vacua of the theory. Since they are space-time independent, 
they correspond to the global minima of the potential. If the potential has a single 
minimum, such as the potential of the </> 4 theory 

V 0) = ^rn 2 4? + ^ 4 , (2-3) 


there is only a single solution to EoM with E = 0 (as one can always arrange things, 
by adding physically unimportant constant to the Lagrangian, so that this condition is 
fulfilled), which is 0 = 0. In this case the vacuum manifold Vi contains just a single 
point and it is, therefore, topologically trivial V\ ~ {0}. In consequence, no topological 
solitons exist. 

However, in the case of a potential 

v{4>) = \{v 2 -<!> 2 ) 2 , (2-4) 


there are two vacua <j) = ±u, giving topologically nontrivial vacuum manifold V 2 . It 
may seem strange to think about a set of two points as a manifold, but geometrically 
speaking, Vo is equivalent to a zero-dimensional circle S° of a radius v. At this point, 
it is hardly necessary to plunge ourselves into technical details of the homotopy theory 
to see, that set of two points is (topologically) different than a set containing just one. 
This topological non-triviality is usually a strong signal that solitonic solutions may 
exist0 And, indeed, in this case there exists a well-known solitonic solution, namely 


4>k(u ) = v tanh 



(2.5) 


What makes this solution topological? In general, any physically acceptable solution 
must approach a vacuum value at a spatial infinity sufficiently fast, otherwise it would 

1 We will always assume that the potential has at least single global minimum with energy E = 0. 

2 It is not a priory guaranteed, however, whether solitonic solutions are physical. For example, a 
topologically non-trivial structure of the so-called 0(3) non-linear sigma model gives rise to a solitonic 
solution called a lump. Upon inspection of its dynamics, it turns out that this solution has a scaling 
instability, meaning that it can evolve into a constant zero or a delta peak in a finite time. This makes 
it unphysical (see [32] for details). 
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2.1 A brane-like solution 


X = 2, v = 1, y 0 = 0 



Figure 2.1: Profiles of the kink solution (2.5) (solid red line) and of its energy density 
(green dashed line). 


have an infinite energy. A non-topological solution approaches the same vacuum ev¬ 
erywhere, which makes it diffeomorphic to it. In other words, it can be transformed 
to a constant by smooth deformations only. On the other hand, topological solution 
takes different values in vacuum manifold as we go along spatial infinity and it is, by 
definition, impossible to deform it to a constant (vacuum). 

Indeed, as y —> — oo function ^>k{v) tends to — v and to v in the opposite infinity 
y —> oo. Therefore, as far as we concentrate only on the spatial infinity in the y- 


direction, the solution 4>K(y) is topological. The solution (2.5) is known as kink or, if it 
is embedded in higher-dimensional model like in our case, a domain wall. The reason 
to call (f>K a “kink” its clearly identifiable from the Fig. |2.1| (red solid line), which shows 
a rapid transition around a point y o and otherwise nearly constant profile. The reason 
to call it a (domain) wall is connected with the fact, that in other coordinates (fixiv ) is 


uniform. Moreover, if we plot its energy density (Fig. 2.1 green dashed line) we observe 
that it is clearly localized around y$, where the transition takes place. This promotes 
the intuition that a kink is indeed a wall. It creates a barrier, which for example, affects 
incoming radiation, etc. 

Another important fact about is its stability. Once a soliton exists, it is very hard 
to destroy, for example by some external forces. In fact, within the formal boundaries 
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of our model, it is impossible. One can quantify this by introducing a conserved charge, 
carried by solitons, called a topological charge: 

OO 

Q =iS^ =m ^ — ■ (2 - 3 * * 6) 

—oo 

where the normalization has been chosen so that Q^k) = 1- Due to the Z 2 symmetry 
4> —> —(/> of our theory, there also exists an anti-kink solution q= —4 >k with the 
corresponding charge Q{4 > k) = — 1- This solution behaves in many ways as an anti¬ 
particle to a kink. For example, a gentle collision of kink and anti-kink results in their 
annihilation into radiation^ 

One can easily construct associated (two-dimensional) conserved current = 
£ llv d u (f>. where ji, u 6 {0, y} and is a completely antisymmetric symbol with e 0y = 1. 
The curious feature of this current is that unlike the Noether current, which is connected 
with some continuous symmetry of the theory, is conserved identically, without a 
need of EoM: 

d p r = = 0 . (2.7) 

Given its wall-like character and its stability, ensured by topological properties, a 
kink provides an excellent field-theoretical model of a brane. 


2.2 BPS argument 


The solution (2.5) has strictly speaking infinite energy with respect to all dimensions, 


but if restricted only to y-coordinate, we obtain a finite number, called a tension 




J dy£(<^r) 


— OO 


2V2X , 

- v 6 


( 2 . 8 ) 


It turns out, that T^ n k is a minimal amount of tension for any static solution of EoM 
with topological charge Q = 1. Let us prove this assertion. In doing so we will employ 

3 This is not, however, the only possibility. In general, the outcome of a collision between kink and 

anti-kink depends very sensitively on their initial velocity (in the center-of-mass frame). Numerical 

investigations of kink and anti-kink collisions in (l+l)-dimensions revealed delicate pattern of reso¬ 

nances, which are formed after the initial annihilation, and which then decay back into kink-anti-kink 
pair. See i.e. [33] for details. 
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a famous trick due to Bogoml’nyi, which will also lead us to the solution to the question 
we have been avoiding so far: How one finds a soliton solution such as (pK? 

Let us write down a tension of a static field <f>(y) 


T = 


d y 


\<t >' 2 + vffl 


(2.9) 


Here <f>' denotes y-derivative. We can rewrite this expression as follows 


T = 

T > 


d y - 
y 2 


OO 

<P'~ J d yy/Wtf)#, (2.10) 


OO 

J d y y/2V{cp)(p'. 


( 2 . 11 ) 


Since the first term in Eq. (2.10) is always positive, one can create a lower bound on 


the tension T, called after its inventor the Bogomol’nyi bound , which is Eq. (2.11) 
Furthermore, if the potential V(cp) can be recast as 


this bound is simply 




T > W(^ 0O ) - W(0_ c 


( 2 . 12 ) 


(2.13) 


where (poo and <p-oo are values of cp at respective infinities. 

The topological nature of this bound is clear. If the field cp interpolates between 
same vacuum values we have (poo = (p-oa and T > 0. For topological solutions, however, 


there is a minimal value of a tension they must have. In the case of the potential (2.4) 
one finds, up to an irrelevant integration constant 


W(cP) = - - <P 2 ) . (2.14) 

Using the fact that a topological solution (up to Z 2 transformation) <^±00 = ±v, we 


obtain for the bound (2.13) 


T > 


2\/2A 3 


= T, 


kink ■ 


(2.15) 
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This completes the proof that any static solution of EoM with Q = 1 (or equally as 
well as Q = —1) must have the tension at least Think- It also provides us with a tool 


how to construct cj>K , since from Eq. (2.10) we see that the bound will be saturated 


(T = Think) if the following first order ordinary differential equation (ODE) is fulfilled: 


-^(v 2 -6 2 ) = 0 . 


(2.16) 


(2.17) 


Eq. (2.16) is called BogomoTnyi-Prasad-Sommerfield (BPS) equation. One can check 


that, by taking its derivative, BPS equation is fully compatible with EoM, meaning 
that any solution to BPS equation is also a solution to EoM. But there is also an 
indirect way how to see this. By saturating the bound we have clearly found a field 
configuration in the minimum of the tension functional T\<p\. However, since for static 
fields a minimum of a tension must also be a minimum of an action, solution to BPS 
equation must also be a solution to EoM. 


It is a straightforward exercise to check that integration of Eq. (2.17) leads to the 


solution (f>k of Eq. (2.5). Solutions of BPS equation are called BPS solitons. We will 


have more to say on BPS solitons in the chapter [3j where we establish important 
connection between BPS solitons and supersymmetry. 


2.3 Trapped modes 

In this subsection we will show that the domain wall solution <f>K can be regarded as a 
field-theoretical model of a brane. This can be established by investigating a spectrum 
of fluctuations of (j>K and showing that few lowest-lying modes are trapped in a well-like 
potential along the extra dimension centered around the position of the wall. Thus, 
these trapped modes can freely propagate only in the four-dimensional world-volume 
of the wall, functionally playing a role of particles living on the brane. 

Let us denote 4>(x fJ, ,y) = 4 >k(u) +e(x M ,y), where n = 0,1,2,3 and |e(x^,y)| <C v. 
Upon pluging this into EoM and retaining only linear terms in e we get 

d\e - A e(v 2 - 3<&) = 0. (2.18) 
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To find normal modes we employ an ansatz e(x^,y) = 6(y)e~ lkx+iajt and plug it back 
into the previous equation to obtain a Schrodinger-like eigenvalue problem for a “wave 
function” b(y) 


-d 2 b(y) + U(y)b(y) = (u 2 -k 2 )b(y), 


where U(y) is a potential given as 


U(y) = Xv 2 


2 sinh 2 (v (y - y 0 )) - 1 


cosh 2 


fty 


y o 


(2.19) 


( 2 . 20 ) 


One can easily check that the so-called zero mod^] 


bo{y) 


Uk(v) = 3 1 / 4 -^ T 

V 2 cosh 2 (v-\f\{y — 2/o)^ 


( 2 . 21 ) 


is a solution to the eigenvalue-problem (2.19) with Wg = k 2 and m\ = 3Au 2 /2. Here 
“zero mode” simply means zero mass. Indeed, this solution corresponds to a freely 
propagating wave in (3+l)-dimensions with a usual dispersion relation; the factor bo(y) 


determines spreading of this wave into the extra dimension. From the Fig. 2.2 (green 
dashed line) we can clearly see, that the zero mode is localized very closely to the 
domain-wall’s position. 

The existence of zero modes is actually dictated by symmetry. Background solution 
4>k, although invariant under translations in four dimensions, spontaneously breaks 
translation invariance in the y-direction. According to the Goldstone theorem, for 
every spontaneously broken symmetry generator, there must be a massless particle in 
the spectrum. Hence a particle associated with the zero mode is often called a Nambu- 
Goldostone (NG) boson. 


We can easily see that our zero mode (2.21) is indeed connected with the breaking 


of the translational invariance in the y-direction. The quantity, parameterizing this 
breaking is, of course, the position of the wall yo. An infinitesimal shift yo —> yo — $yo 
gives 

4>K(yo ~ Syo) = 4>K(yo ) + b 0 (y)dy 0 + o((5y 0 ) 2 ^ . (2.22) 


4 We adjusted overall constant factor so that bo is properly normalized: 


oo 

J Aybliy) = 1 . 
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X = 2, v = 1, y 0 = 0 



Figure 2.2: Profile of the potential U{y) (red solid line) together with unnormalized 
profiles of the zero mode 2 b 0 ( y ) — 1 (green dashed line) and the massive mode 2 bi (y) + 2 
(blue dotted line). 


Therefore, since <pK is a solution to EoM for all values of yo, the zero mode (2.21) must 


be a solution to (2.19). 


Another solution to the eigenvalue-problem (2.19) is a massive mode 


V 2 cosh 2 \vy — yo)J 


(2.23) 


with ujf = k 2 + |Aw 2 , which corresponds to a particle with mass m 2 = 3Av 2 /2. From 
Fig. (2.2) (blue dotted line) we see that, although somewhat wider than in the case 
of the zero mode, profile of the massive mode is definitely localized around the wall’s 
position. 

Being massive a massive mode is hard to excite. If one is interested only in low- 
energy processes it is a common practice to discard all massive modes from the picture, 
in the same spirit as one neglects the existence of massive Kaluza-Klein (KK) modes 
in the KK theory. Unlike the KK scenario, however, in our case there is not an infinite 
tower of massive excitations. The massive mode b\(y) is the only one in the eigenvalue- 


problem (2.19) (as the number of bound states depends very non-trivially on the form 
of the potential U(y)). 
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2.4 Effective theory on a domain-wall 


When the frequency reaches a threshold = 2Xv 2 the spectrum becomes continu¬ 
ous and corresponding particles are no longer localized around the position of the wall. 
Surprisingly the “wave function” of such an unbound state can be found explicitly^ 

b q (y) = e 1TO ^3tanh 2 y — 1 — q 2 — Sigtanhy^ , (2.24) 

The corresponding energy is u q = {2+q 2 /2) Xv 2 . This type of excitation is a “radiation”, 
which can escape from the brane into the bulk or other way around. From the point of 
view of a brane-bound observer an excitation of these particles corresponds to processes 
with missing energy and momentum. 

2.4 Effective theory on a domain-wall 

Investigation of small fluctuations around the domain wall solution revealed two modes 
trapped on its world-volume. We can capture the effective (3+l)-dimensional dynamics 
of these modes (interpreted as particles living on the brane) by constructing an effective 
Lagrangian. Let us again employ the following ansatz 

0(^,2/) = <t>K(y) + e(x lx ,y ). (2.25) 

Let us decompose the general fluctuation field e(x^,y) into the trapped normal modes: 


e(x M , y ) = 6 0 (y)u 0 (^) + h(y)ui(x M ) 


(2.26) 


with bo(y) and b\(y) given in Eqs. (2.21) and (2.23). Now we plug this into the La¬ 
grangian and integrate the result over the extra-dimensional coordinate y 


^eff 


d y 


\d^d^(t) - ^(d y (j)) 2 - 4(v 2 - 4> 2 f \ = Leg + Leff + °( e3 ) > ( 2 - 27 ) 


1 


A 


2 ^ 2 ' Ur) 4^ ^ ) eff r eff 


where 


r(0) _ 
^eff — 


d y 




— OO 
OO 


^eff = j d y ^ed^e- ^(0„e) 2 - ^V"((j) K )e 2 


(2.28) 


(2.29) 


s The wave function b q (y), however, is not a normalizable state and in practice one should consider 
some appropriate wave packet of these states. 
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Notice that first order contributions from the kinetic term and the potential term 
exactly cancel each other, as it should be since the background solution fix lies in the 
minimum of the action. 

The first term can be easily evaluated as 

= v 3 = -T kink , (2.30) 


which is (classically) just unimportant constant equivalent to minus the tension of the 


kink. In the second term, using the decomposition (2.26), equation of motion for the 


normal modes and the fact that they are orthonormal 

OO 

/ d ybi(y)bj(y) = S itj , i,j = 0,1, 


(2.31) 


we obtain 

= ^df.uod^uo + ^d^uid^m - ^m\u\ . (2.32) 

Thus, we found that at the lowest order of approximation, the effective Lagrangian 
contains two noninteracting scalar fields, one of which is massive, with mass m\ = §Au 2 . 
This result is in accordance with our interpretation of trapped modes as particles living 
on a world-volume of a domain wall. Up to this point, however, we ignored continuous 
modes. Such approximation is correct if the total energy of a massive mode does not 
exceed 2\v 2 . If so, radiation would carry the excess away from the brane. Notice, 
however, that similar discussion does not apply to a zero mode, which can have an 
arbitrary velocity (lower than speed of light). This follows from the fact, that a wall 
with an exited zero mode can be interpreted as having a non-zero speed in the y- 
direction. Due to the Lorentz invariance, this should be physically equivalent to a 
static wall. Therefore, no radiation should occur. 

(3) 

We could go on and calculate remaining higher order contributions, namely 
and . These turn out to be 


£ iff = J dyV"'(<l> K )e 3 = - A?r 3 g 4v ^ 1 ( 3 u o Ul + 2 u i) > 

— OO 

r (4) _ 1 f 1 T 7-1111(1 \ 4 _ A3\/3mi / 4 2 2 , 1 4^ 

^eff "p / WKfZ — [Uq + ^ U 0 U 1 + 2 U l) ■ 


(2.33) 

(2.34) 
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2.4 Effective theory on a domain-wall 


No higher corrections exists. Whether we should take these corrections seriously, how¬ 
ever, is a highly non-trivial question. 

Apart from the similar issue with continuum modes, there are deeper reasons to 


be cautious. In Eq. (2.26) we have artificially restricted a space of all possible held 
configurations to a very small subspace. But it is not clear whether such a subspace 
contains all typical held configurations, which emerges if we turn on the full dynamics. 
And indeed, since we ignored continuous modes, we know that for sufficiently high 
energy, it does not. On the other hand, if there is a reason to believe that the held 
stays close to (2.26) in the configuration space, interaction terms in L^ and should 
be reliable. This has to be checked a posteriori, however, for example by numerically 
solving the exact EoM in an energy range of interest. If the energy of fluctuations 
is small, the intuition supports the view, that uo and u\ should be free helds, so the 
leading effective Lagrangian seems to be most reliable. 


In conclusion, for low energy of fluctuations the effective Lagrangian (2.32) is our 
best bet. Moreover, in the case at hand, we can reasonably assume that the mass of 
massive mode m\ is very high and, therefore, we can ignore the held u\ altogether. All 
we are left with, then, is a free zero mode fluctuation propagating along the wall. 

As we will see, in a more complicated theories it is very laborious to derive an 
effective Lagrangian in the direct-fashion way we employed here. Moreover, due to 
reasons similar to those we have just discussed, it is not guaranteed that such a direct 
approach yields correct higher order corrections, as the question of beyond leading 
terms is notoriously difficult. For this reason in this entire text we will concentrate 
on leading order terms only. In fact, we will be satished with an effective description 
of dynamics of zero modes. Fortunately, there is a very beautiful and straightforward 
way, how to obtain it. As we will rely on this method, called the moduli approximation 
invented by Manton [32] . quite heavily in the subsequent discussions, we will briefly 
explain it below. 

Let us first stress once more, that existence of zero modes in the spectrum is closely 
connected with the symmetry breaking by the background solution. Geometrically, 
zero modes corresponds to flat directions in the configuration space, meaning that 
deformation of the field in these directions costs no energy. Therefore, an associated 
mode of fluctuation (particle) is massless. Interestingly, what Manton found is that 
a low-energy dynamics of zero modes can be understood in geometrical terms too. It 
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turns out, that leading order effective Lagrangian is generally given as a sum of kinetic 
terms of zero modes, i.e. 

£ e ff = ^(uo)3 m Uo<9 m Uo , (2.35) 

where i , j runs from all types of zero modes Uq and where gij(uo) is a metric tensor. 
In other words, low energy dynamics of zero modes is equivalent to geodetic motion 
in a curved space, called a moduli space. It is even possible to write down an explicit 
formula for the metric fjjj . Geometrically, if the starting theory has a global symmetry 
group G and the background solution break it down to some subgroup H , the zero 
modes should take values in the coset space G/H . If H is a normal subgroup, then the 
coset G/H itself is a group and its continuous part corresponds to a manifold. Since 
breaking of discrete symmetries does not lead to NG particles, it is only this continuous 
part of G/H which we identify with a moduli space. The metric g t j should be then a 
metric of a continuous part of G/H. 

The formula itself can be written down quite easily, using the following observation. 
Since the background solution breaks a part of the global symmetry of the theory, 


it must contains a number of parameters (such as yo in (2.5)) called moduli , which 
describe the breaking. Let us denote all moduli in our generic example as yj, where 
* = !,..., dim(Cr/R). The metric is then given as (up to an overall factor) 


Sij 


d y 


d(j>b d(/>b 


dyi dyj 


(2.36) 


where is a background solution. Going back to the case of a domain wall solution 


(2.5) we have a single molulus yo an d a metric is just number 

oo 

d<f>K d(j)K 


g ii 


dy 


dyo dy 0 


= Ti 


kink : 


(2.37) 


which can be absorbed into the normalization of the field uq. 

Putting geometry aside, Manton’s method relies on a intuitively obvious picture, 
that low-energy fluctuations of the background solution are deformations of its moduli 
parameters. In other worlds, it is possible to derive effective Lagrangian of zero modes 
just by promoting moduli to fields y* —> yj(x^) on the soliton’s world-volume. Let us 
demonstrate this in our toy model. 
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2.4 Effective theory on a domain-wall 


As the first step we promote the only modulus in the background solution (2.5) 
to a four-dimensional held yo —> yo(x fl ). Now we plug this promoted solution back to 
the Lagrangian and integrate the result over extra-dimensional coordinate y. For the 
moduli approximation, however, it is sufficient just to consider kinetic term, while the 
potential term is dropped: 


^eff 


d y yoix^d^faiy, y 0 (x)) 


-OO 

OO 


= f - 7 =- d^yod^yo = ^.T^kd^yod^yo • 

-L cosh^uJKy-yo)) 2 


1 /2 

If we denote uq = T k - nk yo we obtain 


(2.38) 


7j e ff ^d^uod^uo 


(2.39) 


which is just a Lagrangian of a single massless scalar held, in accordance with our 
previous result. In our toy model, the effective Lagrangian contains single zero mode 
and as we see, the metric is trivial, making the moduli space just M. This indeed 
corresponds to a coset G/H, where G = M is a translational invariance with respect to 
the y-coordinate, which is broken down to the identity H = lj^] 

In chapters |4j [5] and [6j where models we consider becomes progressively more com¬ 
plicated, the moduli approximation method of obtaining effective Lagrangian will be 
one of the keys for the realization of brane world scenario. 


6 One may be tempted to identify G with 0(4,1) the group of space-time symmetries in (4+1)- 
dimensions and H with 0(3,1) the Lorentz group in (3+l)-dimensions, which is preserved by the 
domain wall. Since the number of generators of 0(n — 1,1) is ; it follows that the dimension of 

the coset is dim(0(4, l)/0(3,1)) = 4. Indeed, the domain wall breaks translational symmetry in the 
direction perpendicular to the wall and three rotational symmetries around the three axes tangent to 
the wall. Based on this naive counting, there should be four zero modes, not just one. The reason why 
that is not the case can be intuitively understood by realization that small fluctuation on the domain 
wall can be equally well respected as a local translation or a local rotation. Therefore, zero modes 
connected with the breaking of the rotational invariance are not independent of the translational zero 
mode and there exist really only one independent mode of fluctuation. More informations about this 
issue can be found in (34) . 
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3 

Some aspects of supersymmetry 


Since BPS solitons obey first order rather than second order differential equations, 
namely the BPS equations, they can be often expressed in a closed form. This fact 
greatly reduces the labor needed for derivation of an effective Lagrangian. In particular, 
with a closed formula at hand, we can employ the moduli approximation introduced in 
the previous chapter to obtain the low-energy effective Lagrangian almost effortlessly. 
Therefore, it seems practical to prefer BPS solitons over their non-BPS analogues as 
field theoretical models of a brane. 

However, as we want to localize particles other than scalar fields, presumably we 
need to deal with much more complicated theories than the toy model of Ch. [2] But how 
can we guarantee that solitonic solutions in those theories are going to be BPS solitons? 
More precisely, the question we need to answer is: what is the proper matter content 
of a five-dimensional theory to allow for a BPS domain wall rather just a domain wall? 

The key to answering that question is both simple and surprising and it forms 
the content of this chapter. It is a well-known fact, that if the underlying theory is 
topologically non-trivial and supersymmetric (SUSY), BPS solitons exist. Therefore, 
before we can discuss more advanced model building of Ch. U Ch.[5]and Ch.§ we need 
to introduce few concepts from SUSY and explain their connection with BPS solitons. 

Let us first, however, clarify several things. At the most basic level, SUSY can be 
regarded as an invariance with respect to mixing of bosonic and fermionic degrees of 
freedom. The toy model of Ch. [2] is seemingly not supersymmetric, since it contains 
only bosons. But we can view it as a bosonic part of a supersymmetric theory. As 
far as classical solutions are concerned (i.e. no macroscopic fermionic fields) we can 
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ignore the other part completely. Therefore, our toy model can be enlarged to SUSY 
theory without changing any result of the previous chapter. As we will see repeatedly 
later on, our models are chosen in such a way that they can be enlarged, by adding 
appropriate fermionic and bosonic fields, to be supersymmetric. This property allows 
us to construct BPS solitons, as we will explain below. 

Also, the fact that we use BPS solitons may seem as a simplifying assumption 
about our model (in gauge theories, for example, this corresponds to a special relations 
between gauge coupling and interaction couplings of charged matter fields), but we will 
see that, at least at the lowest order of approximation, this will not affect the result. 
Therefore, even though we will rely on SUSY concepts at the motivational level, our 
final model need not be supersymmetric in order to arrive at the same low-energy 
effective Lagrangian. 

Let us now proceed with a brief introduction of supersymmetry. The first objective 
of this chapter is to clarify the connection of SUSY with BPS solitons. This can be 
most easily/pedagogically done in (l+l)-dimensions. We devote the first section to 
that. The second and ultimate goal of this chapter is to describe general properties of 
SUSY theories in (4+l)-dimensions. This is described in the second section. It will 
turn out, that several objects, namely Fayet-Iliopulous (FI) terms and the so-called 
prepotential, which arise in five-dimensional SUSY gauge theories, are necessary to 
localize gauge fields on the domain wall. We will address these issues in chapter |4j 

3.1 BPS solitons and supersymmetry 

The main goal of this section is to demonstrate that generic SUSY theories in (1+1)- 
dimensions naturally contain BPS domain wall solution. In fact, we will see that their 
bosonic part is essentially the same as the toy model of Ch. HE 

3.1.1 Superspace in (l+l)-dimensions 

There are two ways how to represent SUSY. Either we use ordinary fields, in which 
case we end up having different multiplets of bosonic and fermionic fields, all of which 
fall into the same irreducible representation of SUSY, or we use superfields, which 
are objects transforming naturally under SUSY transformation. The advantage of the 

1 In this section we employ the notation and conventions of [35] . 
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former approach is its clear interpretation, since it involves only known objects, namely 
fields on space-time. On the other hand this representation is SUSY non-manifest and 
it heavily relies on the representation theory, producing rather cumbersome expressions 
which are equally cumbersome to manipulate with. If the theory is formulated in terms 
of superfields, SUSY is manifest and resulting formulas are elegant and concise. The 
cost, however, is the need to introduce new objects, namely superfields, which are fields 
on the so-called superspace. 

A superspace is a combination of ordinary space-time dimensions and additional 
Grassmann dimensions. Thus a superspace in (l+l)-dinrensions has not only the usual 
commuting coordinates x^ G (x^x 1 ), but also two real anti-commuting coordinates 
9 a G (9 1 , 62 ), where a is a spinor index, with the property 

{9 a ,9pj = 8 a p , 6 ^ = 6 a , (3.1) 


where {A, B} = AB + BA. Given their anticommuting nature, the inner product 
between Grassmann coordinates cannot be symmetrical, i.e. 9 2 = 9 a 9 a = 0. In other 
words, the metric tensor must be antisymmetric as well: 


99 = 


A = -2WA , 


(3.2) 


where we introduced Dirac conjugation 9 = O' 7 0 . With this convention Grassmann 
cordinates 9 a are Majorana spinors, if the gamma matrices are in the Majorana repre¬ 
sentation 7 0 = (72 and 7 1 = id 3 , with Pauli matrices given as 


(Tl = 


0 1 
1 0 


o 2 = 


0 -i 
i 0 


0-3 = 


1 0 
0 -1 


(3.3) 


It is easy to check that for two different spinors p’ and (p it holds: 


pxp = 1 pip , ip = —(pj^ip ■ 


(3.4) 


As a consequence we have 9^9 = 0, which is true only in two dimensions. 

Since in superspace we treat both commuting and anti-commuting coordinates 
equally, we must also invent differentiation and integration with respect to Grassmann 
coordinates. Derivatives can be introduced quite naturally: 


d a — J 

d9 a 

* = _ 0 d 

a ^d9 a laf3 d9A 


da@/3 — fia/3 ■> 

d a 9p = 7°^, 

(3.5) 

da0(3 = fia/3 ? 

II 

1 

P O 

"CD 

(3.6) 
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Integration, on the other hand, is defined as 


1 = 0 , 

J d 9 a 9p 

7a 

II 

J d 9 a 9p = 7 °^ , 

(3.7) 

1 = 0 , 

f d9 Q 9fj 

7a 

ii 

7a 

1 

II 

7a 

Qb 

(3.8) 


which may not be everyone’s first guess. In particular, why is f d 9 a -1 = 0 and not 

oo 

f d9 a ■ 1 = 9 a l The reason is that we think of f d9 as a direct analog to f dx rather 

— OO 

than to an anti-derivative. One of characteristic properties of a line integral in ordinary 
space is that a constant shift of an integration variable does not change the result. If we 


demand that the analog holds for 6 integration, we obtain exactly the properties (3.7) 


(3.8). Moreover, comparing (3.7)-(3.8) with (3.5)-(3.6) we see that fomally f d 9 a = d a , 
from which f d 9 a ■ 1 = d a ■ 1 = 0 follows immediately. 

Let us also define a volume integral as 


d 2 9 = / d9 a d9 a = 2i / d9od9\ = —dd = 2i<92<9i 


(3.9) 


And in particular 


d 2 9 99 = -3899 = 4. 


(3.10) 


3.1.2 SUSY generators and superfields 

Let us look at symmetry transformations of superspace coordinates which leaves the 
volume measure d 2 xd 2 # invariant. These obviously include Lorentz transformations 
and constant shifts of space-time coordinates x^ and in analogy, transformations that 
mix together anti-commuting coordinates 9 a . Both of these are not particularly impor¬ 
tant for the following discussion. A physically interesting symmetry mixes both type 
of coordinates together. Let us introduce a SUSY transformation 5 £ in the following 
way 

5^0 a = £« , = -\9^i , (3.11) 

where £ is a constant Majorana spinor parametrizing the SUSY transformation. From 


Eq. (3.11) we can extract the formulas for fermionic operators Q a , which generate the 


SUSY transformation 5^ = £ a Q a = Qafia 


Qa = d a + i, Qa = -d a - i Oplfadu 


(3.12) 
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It can be easily checked that these generators satisfy an algebra 


{QonQp} = - 2i ($) a/3 ! 


(3.13) 


where the Feynman slash notation is expanded as $ = 'y^d^. All other anti-commutators 


are zero. Moreover, one can check that adding Qi,Q 2 to the list of generators of 
Poincare group of symmetries of (l+l)-dimensional space-time results in a well-behaving 
graded Lie algebra (or simply superalgebra ), which is closed under repetitive use of 
(anti)commutators between all its elements. In other words, SUSY generators extends 
the usual space-time algebra in a non-trivial way (i.e. not as a direct product). 

Back in 1967 a famous no-go theorem of Coleman and Mandula (36| stated that 
exactly this kind of mixing of internal and space-time symmetries is impossible. To 
be more precise, the Coleman-Mandula theorem claims that any realistic theory (non¬ 
trivial S-matrix) with a mass gap in (3+l)-dimensions can only have symmetry Lie 
algebra which is always a direct product of the Poincare group and some internal 
group. As it stands, the theorem is indeed correct. However, it turned out that 
assumptions of the theorem are too restrictive. In 1971 Golfand and Likhtman m 
explicitly constructed a non-trivial extension of the Poincare algebra by adding to it a 
set of fermionic operators (forbidden in Coleman-Mandula analysis), which marked the 
beginning of SUSY. Four years after this breakthrough the search of possible extensions 
of the Poincare algebra culminated in the so-called Haag-Lopusanski-Sohnius theorem 
[35] , showing that under fairly general conditions SUSY is the only extension beyond 
the four-dimensional Poincare algebra. 

After this historical sketch let us now proceed to introduce another vital concept of 
SUSY: a superfield. As far as we are interested in supersymmetric theories in (1+1)- 
dimensions, there is only one type of superfield <F = <F(x, 8 ), which is simply an arbitrary 
function of superspace coordinates. In more dimensions, where superspace description 
is in general redundant, there are more types of superfields. 

Given their anti-commuting nature, any function of Grassmann coordinates 9 a can 
be Taylor-expanded into a finite series of terms. Thus, a superfield can be always 
decomposed into component fields as 



(3.14) 
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where cj)(x),F(x) are scalar fields and ip(x) is a Majorana spinor field. These compo¬ 
nents are various projections of a superfield 


= 4 * 


6»=0 


Ipa = d a $ 


, F = - l dd§ 
0=0 2 


0=0 


(3.15) 


and they are independent of each other. Under SUSY transformation, a (scalar) super¬ 
field only changes by change of its arguments: 

= <h(x + S^x. 9 + 5^9) — <h(x, 9) = ^ + 9{F — ■ (3.16) 


Finally, matching up the powers of 9 we can express the SUSY transformation in terms 
of components as 


$t<i> = i'P , 

(3.17) 

5^ = -i+ F£ , 

(3.18) 

. 

(3.19) 


3.1.3 SUSY theories in (l+l)-dimensions 


Notice that the 99 component of a superfield transforms into a total derivative: 6 ^F = 
i5 At ('07 /i ^). This observation allows us to construct supersymmetric theories in a very 
simple fashion, since its immediate consequence is that integration of the F-term over 
space-time is a SUSY scalar. Indeed, 

<5g J d 2 xF = i j d 2 xd tl ('ip'y fJj p t ) = 0 , (3.20) 


where the last equality follows from the Stokes theorem, provided that F vanishes at 


the boundary. Using (3.10) and (3.14) we can reformulate this observation in a SUSY 
invariant way as 


i2 * F =\ 


d 2 xd 2 9Q = SUSY scalar. 


(3.21) 


This result is rather pleasing. It tells us that if we integrate a superfield over the 
whole superspace we obtain a SUSY invariant quantity (which is by construction also 
a Lorentz invariant). Let us stress, however, that the above relation holds for any 
superfield. As a consequence, the way how to write down an action for SUSY invariant 
theory in (l+l)-dimensions is simply 


S 


d 2 x£ = - / d?xd 2 9Ls, 


(3.22) 
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where £5 is a superfield. Factor ^ is added for convenience. 

Let us in the following reserve the symbol to represent a superfield, whose com¬ 
ponents are fundamental degrees of freedom of our theory, namely <f>, if) and F. We 
can use to construct other superfields. In particular, any function of <L, say /(<!>), is 
also a superfield. We would also like to construct terms, which on the component level 
contain quantities such as d^cf) and and are superfields at the superspace level. It 
turns out, that individual derivatives and <9 a <f> do not transform in a desired way 
under SUSY transformation. However, if we combine them into a covariant derivative 


D a = d a - i( 7 ^) 0 ^ , D a = -d a + i(6y%<9 M , 


(3.23) 


the associated quantities D a < f> = 7 ° Q ,.D/ 3 < 1 > are superfields. Therefore any Lagrangian 
density of the form £s(<I>, D a &), with properly contracted Lorentz and spinor indices, 
is a viable candidate for supersymmetric theory in (l+l)-dimensions{^] In addition, if 
we demand that component fields are to have canonical kinetic terms, we end up with 
a general expression 

£s = \D a $D a $ + W($), (3.24) 

where the so-called superpotential W(<3?) is an arbitrary function of <f>. 


We can now Taylor-expand the supersymmetric Lagrangian density Lg of Eq. (3.24) 
with respect to 8 coordinates and perform the integration f d 2 9, which picks out the 
88 component of Lg. The result is 


1 


£ + ^F 2 + 


1 


dW(</>) 1 d 2 W(</>) - 


-F - 

2 


d (j) 2 


ipip. 


(3.25) 


Since the F component is not dynamical it can be eliminated through its EoM 


dL 

~dF 


= 0 . 


F = — 


dW($) 

d4> 


(3.26) 


which leads to the final form 


£ = -d^d 11 ^ + - 




2 V d<L J 


2 


2 V’V’- 


(3.27) 


2 Arguments coming from dimensional analysis and locality further reduce this general form into 

C s = ^X($)D a $D a $ + W(4>), 


with both X and W being arbitrary functions of <£>. 
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The bosonic part of the Lagrangian (3.27) (i/j = 0) is in fact the one, which we used 


as our toy model in the previous chapter, with the superpotential given in Eq. (2.14). 
Thus, we confirmed the claim we made at the beginning of this section. Let us now 
settle the remaining issue, namely of the interplay between SUSY and BPS solitons. 
We will discover that a BPS solution has a unique property of preserving a part of the 
SUSY. A property which will not only provide an alternative way of constructing BPS 
solitons, but also automatically localize fermions. 


3.1.4 BPS solitons and breaking of SUSY 


Let us discuss what kind of solutions the theory (3.24) can have, with respect to their 
behavior under SUSY transformation. The equations of motion can be written down 
either in the superfield formalism as 

-DD$- 
2 


or by using components 


d 2 (j) — 


d 2 W 


F + - 


dW n 
d<F “ ’ 

(3.28) 

1 d 3 W - 

(3.29) 

2 d^' W ’-°' 

d 2 W 

(3.30) 

’ d^- 0 ’ 

dW 

(3.31) 

h ~ do 


Both ways are, of course, equivalent, which can be checked by 9 expansion of the 
superfield EoM, with the help of the identity 


DD = 88 + i{9^8 + 8^0)8,, + 968 2 . 


(3.32) 


A generic solution of the above system of partial differential equations breaks 
SUSY, meaning that 

Qi<f>/0, Q 2 $ / 0. (3.33) 

Therefore, if we took <f> as a background solution, the effective theory of its fluctuations 
will not be supersymmetric (moreover, as the V’a component of is generically non¬ 
zero, even the Lorentz invariance will be broken). But is this true for all solutions? 
Certainly not. It can be easily checked that one type of solutions, let us denote them 
as 4>°, are SUSY invariant. This means that apart being solutions to EoM, they also 
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satisfy conditions Qi<f > 0 = Q 2 ^° = 0. Solving for those conditions we discover that 
<h° = </P, where (fP are vacua of the theory ^ = 0 . 

Since we have two SUSY generators, there exists a third possibility, that some 
solution (h 1 / 2 breaks only one generator, for instance the first one Q\Q l P ^ 0 , while 
the second one is preserved Q 2 <& l P = 0. Expansion of the last condition to components 
yields 

Q2$ 1/2 = (02 + 2 + H 1 ' 2 + \oOF l ' 2 ) 

= V>2 /2 + 0 2 F 1/2 + i^J a d^4> l/2 - ^OO^d^U 2 . (3.34) 

Setting each order of 9 a to zero, we obtain 

d t fi 1/2 = 0 , F 1/2 + d x (j) 1/2 = 0, (3.35) 

d t ^\ /2 = o, v4 /2 = 0 • (3.36) 


We have to check whether these conditions are compatible with EoM. In particular we 
have F 1 / 2 = — 4^, which together with the first and second condition gives 


8rfl/2 (u = % 


cf)=(f ) 1 / 2 


(3.37) 


This is nothing else but the BPS equation (2.16)! We can see this directly by recalling 
Eq. (2.12). For the ipl / 2 the 1/2 SUSY preserving conditions (3.36) demand that = 
0 and that i ^ 2 = ^ 2 (x) is an arbitrary function of x = x 1 . The exact x dependence 

of can be easily obtained from EoM (3.30) and it is uniquely determined through 
the superpotential by (jpP^x ) up to an integration constant: 




1/2 

1 


= e 


rx d 3 W 
J d<fr2 rjj 


(3.38) 


where is a constant Grassmann number. Thus, we can always consistently set i\) = 0, 
which we will always do for classical solutions of EoM. If we would allow for nonzero 
i / 1 / 2 as a part of the background solution, resulting effective theory will not be Lorentz 
invariant. Let us also point out that setting i / 1//2 to zero allows us to recast superfield 
solution T 1 / 2 into a particularly nice form 

& 1 ! 2 = <j> K (x - \06), (3.39) 
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where <j>K is a kink solution of Eq. (2.5). If we demand that the first generator is 


preserved Qi^ 1 / 2 = 0, we simply obtain the anti-kink solution and the corresponding 
superfield is d) 1 / 2 = —c/)k(x + ^00). 


In the background of the BPS solution <3? 1 / 2 of Eq. (3.39) a half of the SUSY 


generators is spontaneously broken. Therefore we expect that in the effective theory 
there will be a corresponding Goldstone fermion or goldstino. This new particle is a 
superpartner to a bosonic zero mode, which arises from breaking of the translational 
symmetry. This observation can be easily verified by direct calculations. If we write 


down linearized EoM for fermionic fluctuations S'ip a = e luJt c a (x), with use of Eqs. (2.5) 


and (2.14) and setting A = 2, v = 1 for convenience, we obtain 


d x c\ — \UC 2 — 2 tanh(x)ci = 0, 
d x C 2 — itaci + 2 tanh(x)c 2 = 0 . 


(3.40) 

(3.41) 


This system has a non-trivial solution only if u> = 0, for ci, C 2 are real. Moreover, from 


general conditions (3.36) we have ci = 0 and solving for C 2 yields 

1 


C2 = 


cosh 2 (a?) 


V, 


(3.42) 


where r) is a constant Grassmann number. The profile of C 2 (x) is exactly the same as 


the profile of a bosonic zero mode (2.21). Thus, we have found a localized fermionic 


fluctuations, which is the sought goldstino. Let us stress, that similar analysis can be 
applied to (4+l)-dimensional case in the same fashion as we did for the toy model 
of Ch. [2j The same goldstino localizes on (3+1 (-dimensional brane, together with 
previously found scalar zero mode. The resulting effective action, however, will not 
be supersymmetric in (3+l)-dimensional sense, but only in (0+l)-dimensional sense. 
This is so, because SUSY of the parent theory is a rather special one, as it completely 
ignores three spatial dimension. Indeed, in the next section, we will set up a theory 
with a genuine (4+1 (-dimensional SUSY, which will give us much richer spectrum of 
particles in the effective theory. 

In conclusion, we have arrived at very important fact: BPS solutions partially 
preserve supersymmetry. This characteristic property of BPS solitons is universal and 
it is not, for instance, restricted to (l+l)-dimensions. It also gives us an alternative way 
how to find solitonic solutions. The Bogomol’nyi trick we introduced in the previous 
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chapter is not so straight-forward to implement, as it involves some guessing work and 
it may become cumbersome as the number of field in the model increases. On the 
other hand, partial SUSY preserving conditions are just a set of first order differential 
equations, from which one can read off corresponding BPS equations upon inspection. 

It is also worth mentioning that on the quantum level topological charges become 
central charges of the superalgebra (see e.g. j39| for details). 


3.2 Supersymmetry in higher dimensions 


In order to make our discussion of localization of gauge fields of Ch. [4] unimpeded by 
technical details, it is convenient to develop basic understanding of SUSY theories in 
higher space-time dimensions. In this section, we are going to present basic facts about 
supersymmetry with four and subsequently eight supercharges. Our exposition will 
be focused on the model building aspects, namely on the matter content of generic 
N = 1 SUSY theories in (3+l)-dimensions and (4+l)-dimensions. This will help us to 
cement our notation used in the rest of this text and also to expose vital field-theoretical 
ingredients, required for localization of gauge fields (see section 4.3 in Ch. [4]). A more 
detailed discussion of supersymmetry in four and five dimensions can be found in 

55 , 52 , 53 ], @ 13 , to name just a few. 


3.2.1 A SUSY invariant theories in (3+l)-dimensions 

A superspace realization of N = 1 SUSY in (3+l)-dimensions is obtained by adding 
four Grassmann coordinates 9 a , 6 „ to the space-time coordinates x 12 , with // = 0,1,2,3. 
Here a = 1,2 and d = l,2 are both spinor indices, but while first one is reserved for 
spinors in (^,0) representation of the Lorentz group, the second one is used for spinors 
in (0, g) representation. The dotting notation makes explicit that indices from different 
representations should never be contracted together. This also means that although 9 a 
and 9a look formally conjugated to each other, they are in fact independent objects. 
Raising and lowering of spinor indices is performed by e symbols as 9 a = e a ^9p and 
Qa _ where conventions are set in such a way that 

9 2 = 9 a 9 a = -2 9 1 9 2 , 9 2 = M" = 2<9j 9 2 , (3.43) 
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holds. Following identities are frequently used 

e°e i3 = ^ e a/ V, e 2 , e^eea u e = \ g ^e 2 e 2 , (3.44) 

where 

0o li e = = -OaG^Oa = -e^e , (3.45) 

with a 11 = (1, (Tj) and a 11 = £ T a^£ = (1, —a j) being four-vector extension of the Pauli 


matrices given in Eq. (3.3). The sign convention for the metric tensor is (+, —, —, —). 

Under supersymmetry transformation 5^ superspace coordinates (x’ m ,0 q ,0q) un¬ 

dergo a translation 


5^ £ = i#cr M £ + i0<r M £ , 

(3.46) 


(3.47) 

Og 1" Qa = Ca ) 

(3.48) 


where £, £ are independent Weyl spinors parameterizing SUSY transformation. This 
transformation can be represented as an operator on the superspace ^ f = C a Qa+CaQ a > 
where the supercharges Q a and Qa are given as 

’ <3° = - Js + ie “<A ■ (3.49) 

Correspondingly, the covariant derivatives are 

A. = , A. = - J- - i»XA • (3-50) 


3.2.1.1 Chiral superfields 

A superfield is in general an arbitrary function of superspace coordinates f(x^, 8 a ,9a)- 
However, we are at liberty to impose covariant constrains on a superfield. It is known 
that a chiral superfield, denoted as 4>, defined by the constraint 


Da$ = 0 , 


(3.51) 


forms an irreducible representation. The above constraint can be easily solved, if one 
realizes that there are two primitive superfields, namely y^ + iOa^O and 9 a , which 

are chiral superfields 

D,ay^ = Da ( X» + Wa^e) = 0 , Da6 a = 0 . (3.52) 
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An arbitrary function of these primitive chiral superfields <f>(y, 0) is also a chiral super¬ 
field. We can expand <3?(y, 9) in terms of components as 

<%, 0) = <f>(y) + V29^(y) + 9 2 F(y) (3.53) 


and further expansion of Grassmann variables inside y yields 

<f >(x,9,9) = 4>(x) + ida tJ '9d^ l (j)(x) — ^9 2 9 2 d 2 cj>(x) 

+ y/20i[) H — l -=6 2 9a^d^(x) + 9 2 F(x) . (3.54) 

V2 

We see that a chiral multiplet corresponding to this chiral superfield consists of a 
complex scalar </>, Weyl spinor ^ and an auxiliary held F. Similarly, an anti-chiral 
superfield is defined by the constraint D a = 0 


d>V, 9) = + V29^(y*) + 9 2 F\y *). 

&(x,9,9) = <ft(x) — \9o^9d^ {x) — ^9 2 9 2 d 2 ^(x) 

+ — l j=9 2 9a ll d fl 'ijj(x) + 9 2 F\x). (3.55) 

v2 

A generic SUSY invariant theory constructed purely out of chiral and anti-chiral su¬ 
perfield can be written as 


L = 


d 2 9 d 2 9 + 


d 2 0W($) + 


d 2 0W t ($ t ), 


(3.56) 


where W(4>) is a superpotential. The reasoning is as follows. First, any function 
of a chiral superfield is obviously also a chiral superfield. Since the 9 2 component 
of a chiral superfield (or 9 2 component of an anti-chiral superfield) transforms into 


a total derivative under SUSY transformation, the superpotential part of (3.56) is 
indeed SUSY invariant. Second, a product of chiral and anti-chiral superfield T A 3 is a 
general superfield. Consequently its 9 2 9 2 component transforms into a total derivative. 
Moreover, at the component level (up to a surface term) 


$ 1 $ 


e 2 e 2 


d^d^cj) + + F^F , 


which gives canonical kinetic terms for the fields cj) and ■!/>. 


(3.57) 
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Let us note that is not the only possibility for acceptable SUSY invariant 
kinetic term. In fact, considering more chiral fields for a moment <3?j, with i = 1,... N, 
the most general structure is given by 

J d 2 9d 2 0X(<$>,ft), (3.58) 


where 3>t) is the so-called Kahler potential. The reason why it is called a “poten¬ 
tial” is that at the component level the second derivatives of %(<!>, <3?t) 




0 2 9 2 


d 2 X 


dftdcj) 


jd^ftd^ftj 


+ 


(3.59) 


second derivatives of X((f >, ft) determines the metric tensor in a target space of ft fields. 
Moreover, as for any potential, there is certain amount of arbitrariness in definition of 
3C(4>, <K). It is easy to see that adding any superfield of the form /(<h) + ft (4?t) to the 
Kahler potential A(<h,<fd) gives exactly the same actionj^] An invariance under such 
a transformation (incidentally called a Kahler transformation) hints to the underlying 
geometrical structure of IV = 1 SUSY theories in (3+l)-dimensions, that a target 
space of bosonic fields 4>i must be a special type of a complex manifold called a Kahler 
manifold [33] . 


3.2.1.2 Vector superfields 

A vector superfield V satisfies a reality constraint V = V'. Let us first consider the 
Abelian case. If we make the replacement 

V —>• V + i(A — A f ), (3.61) 


where A (A^) is a chiral (anti-chiral) superfield, we obtain another vector superfield. It 
can be shown that this transformation coincides for a particular component of V with 
the usual gauge transformation. Therefore V is considered as a SUSY extension of a 


gauge field and (3.61) as a SUSY generalization of a gauge transformation. In the same 


3 The reason is that a superspace integral of the type 

d 4 xd 2 6d 2 6 f($) 


/' 


(3.60) 


vanishes, since one can first make a change of variables x 11 —> y M and then use the fact that /($) has 
no 9 2 term. A similar argument can be applied to / 3 ($'). 
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way as we fix the gauge for a gauge field, we can use (3.61) to reduce the number of 
components of V. The usual choice is the so-called Wess-Zumino (WZ) gauge [ifij . in 
which V is decomposed as 


V 


= —Qa^QAn + y6 2 0\ - \6 2 6\ + -0 2 6 2 D , 

wz 2 


(3.62) 


where A M is a U( 1) gauge field, A is a Weyl spinor (gaugino) and D is an auxiliary field. 
Let us point out that the WZ gauge does not fix all the gauge freedom and the usual 
gauge transformations A M —> A p + d^a are still allowed. The SUSY generalization of a 
U( 1) field strength is defined as 

W a = -^D 2 D a V, W ( , = -^D 2 D 6l V . (3.63) 


Using the identity 

{D a ,Da} = -2ia£ & d lt (3.64) 

it can be easily verified that field strength W a (and correspondingly W^) is gauge 
invariant and chiral (anti-chiral) superfield. Thus, in terms of components in the WZ 
gauge (with y p = x p + i6 a p 6) 

W a = -iA a (y) + O a D(y) - ^(a p a u 9) a F fJU (y) + 0 2 (cr M 3 M A(y )) q , (3.65) 


where F pu = d p A u — d v A p is an Abelian field strength. 

Since iU Q is a chiral superfield, W a W a is also a chiral superfield and its 0 2 compo¬ 
nent can be used as a SUSY invariant term. In particular we have 

T(/ d Hw°w a + Jtfewtw*) = + F &_ Ua<,%a. (3 ,66) 

Also, let us note that inclusion of the so-called #-term (not to confuse with Grassmann 
coordinate!) can be achieved in a very compact way 

4lm(r J dWnU) + (3.67) 


where r = 6 /n + &iri/g 2 and F pi ' = e tJjl,prT F pa . 

Let us now write down SUSY Lagrangian with local U( 1) symmetry and interacting 


matter fields belonging to the chiral multiplet. We first need to modify (3.56) to be 


invariant under local 17(1) symmetry. Let us note that the Lagrangian (3.56) is already 
invariant under a global 17(1) symmetry, which acts on the superfields as 4> —> e lA< h 
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and —>• e~ lA 4>t, where A £ M^] It turns out that to make this symmetry local, it is 
necessary to promote A to a full chiral superfield and replace the kinetic term of (3.56) 
with 

$ f e“ v $. (3.68) 


Therefore, a simplest supersymmetric 17(1) gauge theory with a single chiral multiplet 
is defined: 

£ = _L (y d 2 0 + c.c.) + J d 2 0 d 2 * e &e~ v § . (3.69) 

The generalization of the above Lagrangian to contain more chiral superfields is straight¬ 
forward. 

Similar considerations are applied in the non-Abelian case. The vector superfield 
V = V a T a takes values in the adjoint representation of the gauge group and transforms 
under a gauge transformation as 


e -V^ e iAt e -V e -iA_ 


(3.70) 


The normalization of the generators T a is chosen as Tr[T a T & ] = ^5 ab . The chiral 
superfields are in the fundamental representation of the gauge group and transforms as 


4> ->• e lA <h , ->• $te lAt . 


(3.71) 


Notice that the kinetic term (3.68) is still invariant under this transformation, if the 
ordering of superfields is exactly as denoted. The non-Abelian field strength is defined 
as 

W a = ^B 2 e v D a e~ v , (3.72) 

which is still a chiral superfield, but it transforms as an adjoint representation 

W a -> e iA W a e _iA . (3.73) 


The Lagrangian (3.69) can be readily modified to non-Abelian case just by includ¬ 
ing trace in the kinetic term for the vector superfield. Notice, however, that now 
represents a column vector of chiral superfields. 


£ = d 2 6Ti[W a W a ] +c.c.) + J d 2 0d 2 9&e- V <$>. (3.74) 

4 It can be shown that requiring global 17(1) invariance and renormalizibility severely restricts pos¬ 

sible forms of a superpotential W($). If one allows for multiple chiral superfields the superpotential 

can be at most a cubic function [40] , In the case of a single chiral superfield, however, the superpotential 
is zero. 
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In terms of components we obtain 


+ - ii 


— Tr D + ix/ih/^A — iV^^X — FF t] 
dW dW t 1 d 2 W 1 d 2 W 

d^i d(j>\ 1 2 d4>id(j)j J 2 




(3.75) 


where index i labels the fundamental representation of the gauge group and are co¬ 
variant derivatives. One can easily include 9 term in the above Lagrangian by replacing 


^2 (J d 2 9Ti-[W a W Q ] +c.c.) -> ^Im^rTr j d 2 0IU"iy a ) . (3.76) 

The auxiliary Helds F t and D a can be eliminated through EoM. The resulting scalar 


potential is given as 


v{4>) = £ 

i 


dW 

dfa 


(^T a <t>) 2 • 


(3.77) 


3.2.1.3 Fayet-Iliopoulos term 

The 17-term of a vector superfield has some unique properties. It is gauge covariant and 
being 6 2 6 2 component of a superfield, it transforms into a total derivative under SUSY 
transformation. Thus, one can always consider the so-called Fayet-Iliopoulos (FI) term 


cjd9 2 d9 2 V = cD (3.78) 

as a valid term of a SUSY invariant theory. Notice that the non-Abelian analog would 
be cTr[I7]. This is non-zero only for those components of D , which corresponds to a 
17(1) subgroup the gauge group. Thus, there can be as many different FI terms as the 
number of different 17(1) gauge groups in the theory. 

Fayet and Iliopoulos have shown [46], that presence of such terms in (3+l)-dimensional 
SUSY gauge theories leads to the spontaneous breakdown of supersymmetry. In the 
Abelian gauge theories in the five-dimensional space-time dimensions, however, FI 
terms do not break SUSY and can be freely included into the theory. This turns 
out to be important for the construction of domain walls, since constants appearing in 
front of FI terms, incidentally called FI constants, are related to vacuum expectation 
values of Higgs fields. 
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3.2.2 A SUSY invariant theories in (4+l)-dimensions 

Our discussion of SUSY has been, so far, based on the superspace formalism, which is 
unarguably very powerful in (3+l)-dimensions or less. Interestingly, realization of the 
superspace idea in higher dimensions turns out to be rather non-trivial. To be precise, 
the practical usefulness of maintaining manifest SUSY via the superspace formalism 
becomes questionable, when the number of supercharges is larger than four. For most 
of these situations, it is just more efficient to work directly with components. However, 
there are interesting alternatives in the case of eight supercharges, like the harmonic 
superspace approach m or closely related projective superspace approach 32 ] 
(see [50] for a review). A more conservative possibility is to keep SUSY manifest 
only partially and use the language of N = 1 superfields in four dimensions, which 
we developed in the previous subsection. This has been achieved in the work m, 
where five-dimensional N = 1 SUSY gauge theories were described in terms of four¬ 
dimensional superfields, but only for the price of sacrificing manifest Lorentz invariance. 

In this section, we will adopt similar, but less direct approach. Our goal is to 
write down (a bosonic part of) a generic N = 1 SUSY non-Abelian gauge theory 
with minimally interacting matters in (4+l)-dimensions. It will be sufficient for our 
purposes to obtain the result in terms of component fields. However, instead of just 
citing the result, we will motivate it by drawing an analogy between N = 1 SUSY in five 
dimensions and the so-called extended N = 2 SUSY in four dimensions. Formally, there 
are many similarities between these types of models. Both posses eight supercharges 
and both have nearly identical structure of SUSY multiplets. We will use this similarity 
to our advantage and approach the problem in the following way. First, we establish 
how to write down SUSY invariant Lagrangians for N = 2 SUSY in (3+l)-dimensions, 
through N = 1 superfield formalism. In this part we rely heavily on the work ||42J. After 
that we write down these Lagrangians in terms of componens and comment about how 
to elevate them to five dimensions. 

There are two kinds of multiplets in N = 2 four-dimensional SUSY. The N = 2 
chiral multiplet, also called a hypermultiplet, consists of a N = 1 chiral multiplet 
3 Q and anti-chiral multiplet 3 Q\ where Q,Q t denotes cor¬ 

responding superfields. On the other hand, the N = 2 vector multiplet is made of 
a chiral multiplet (E,£, Y) 3 and a N = 1 vector multiplet (W^,,X,D) 3 V. Let 
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us now develop an IV = 1 SUSY manifest way, how to describe these multiplets and 
how to construct supersymmetric SU(N C ) gauge theory with Nf minimally coupled 
hypermultiplets in the fundamental representation. In order to achieve it, let us briefly 
study both kinds multiplets and their Lagrangians separately. 


3.2.2.1 Vector multiplet 

In addition to the space-time coordinates, the N = 2 superspace is made of eight 
Grassmann coordinates 9 a ,0a, 9^,9^- A generic superfield f(x, 9 , 9, 9, 9) contains many 
redundant degrees of freedom and one has to impose two or more constraints on it, 
in order to obtain irreducible representations. It turns out (see i.e. 02!, 03]) that an 
N = 2 vector multiplet is associated with the generalization of N = 1 chiral superfield, 
or in other words, a superfield T constrained by conditions 


Da^ = 0 , Da '£ = 0 . 


(3.79) 


The covariant derivatives D are the same as in Eq. (3.50) with 9 replaced by 9 every¬ 


where. Both constrains can be solved in analogy to the N = 1 case, by having T depend 
only on primitive chiral superfields (y, 9, 6), where y^ = y^+i9a^9 = x fJ ’+i9(j^ L 9+i9a fl 9. 
This leads to the decomposition 


T(y, 9, 9) = (y, 9) + V29 a ^ (y, 9) + O 2 ^ (y, 0) . (3.80) 


After additional constrains are imposed in order to remove certain unphysical degrees 
of freedom, T can be expressed in terms of N = 1 superfields as [43] 

*(y, 0 , 0) = 4(y, 9) + V29 a W a (y , 9) + 9 2 G(y , 9 ), (3.81) 


where W a is given as in Eq. (3.72) and 


G(y, 9) = J d 2 0 [$(y - i 9a9, 9, 0)] f e - v ('J-iOaO/J/J) (3.82) 


The d 2 0 integration is meant at fixed y with ^(x, 0, 9) given in Eq. (3.54) and V ( x , 9, 9) 
in Eq. ( |3.62 ). We see that N = 2 chiral superfield (also called a vector superfield) really 
contains (if in a rather complicated way) only two N = 1 superfields, that is the chiral 
superfield and the vector superfield V, as claimed in the beginning of this subsection. 
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The most general N = 2 Lagrangian containing SU(N C ) vector multiplet, can be 
written as 

£ = -j^—Im^Tr J d 2 0d 2 0J(^)^ , (3.83) 

where T('I') is the so-called prepotential. In four dimensions the prepotential can be for 
renormalizable theories at most quadratic function of 4>. If we integrate out 6 variables, 
we obtain the expansion 


£ = —Ini 
47T 


1 [ d 2 C &S 

2 / d<5> a dd> b 


W a 


*W b + 


d 2 ed 2 e 




-v 




(3.84) 


where a = 1,..., N 2 — 1 labels generators of SU(N C ) group. Let us write down the 
bosonic part of this Lagrangian. For that we choose the prepotential in the form 


The result is 


T(4*) = ^$ a 4> Q 

9 


(3.85) 


£ = + X -D 2 - E f [D, E] + F^j , (3.86) 

where the trace applies for the gauge index. To turn this expression into the bosonic 
part of a five-dimensional SUSY gauge theory, we have to elevate partial derivatives 
and gauge fields into their five-dimensional counterparts and also consider only real 
adjoint scalar fields E^ a = E a , since one DoF is reserved for one of components of the 
five-dimensional gauge field. Notice that this eliminates the fourth term in the above 
Lagrangian (by virtue of properties of the trace). We can now combine auxiliary fields 
F a and D a into SU(2)r triplet, which we call Y£, where a = 1,2,3. SU( 2)r is an 
internal symmetry of superspace, which mixes tilded and untilded quantities together. 
By taking the SU(2)r frame, so that only the third component a = 3 is nonzero, we 
can write = Y a and Yf 1 = Y£ = 0. Thus, in five dimensions we arrive at generic 
N = 1 SUSY Lagrangian for a vector multiplet 

£ S YM = ^ a6 (^) {-\F a MN F bMN + X -DuF a D M E 6 + W a Y b ^j , (3.87) 

where 3“ a 6 (4?) = d 2 3'/d& a d& b . But the situation in five dimensions is even more com¬ 
plicated than that. It follows from general principles that 3“ can be cubic function of 
chiral superfields m and careful analysis of SUSY transformation of 3'( < L) reveals that 
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also five-dimensional Chern-Simon (CS) term, which is proportional to third derivatives 
of the prepotential, might be included 


£cs = ^ a6c ($) 


-^ MNPQR W a M (F b NP F c QR + l - [W N , Wp] b Fq R 


+ ^[w N ,w P ] b [w Q ,w R Y 


(3.88) 


where T a6c (<h) = d 3 3 / d$ a d$ b d<&°. Lastly, as we discussed in the previous subsection, 
if the gauge group of the theory includes some t/(l)’s it is always possible to add 
corresponding FI terms to £>sym- 


3.2.2.2 Hypermultiplet 

As matter fields and gauge fields transform under different representation of the gauge 
group, they must belong to different multiplets. As discussed at the beginning of this 
subsection, hypermultiplet can be understood as a pair of N = 1 chiral multiplets 
Q, Q. In order to write down N = 2 SUSY Lagrangian with Nf minimally coupled 
hypermultiplets, we have to understand how kinetic terms and potential terms of both 
chiral superfields must be combined to ensure N = 2 SUSY invariance. It turns out that 


this is relatively easy. The kinetic term (3.68) for minimally coupled chiral superfield 


is unchanged, apart from ensuring that normalization of both Q and Q terms is the 
same. The superpotential, however, is restricted in N = 2 SUSY. The reason is that Q 
and Q transform as a doublet under SU(2) R symmetry, which is an internal symmetry 
of the superspace. The term such as W(Q) would obviously break it. Furthermore, the 
interaction between Q and Q are restricted only to quadratic (mass) term. Because 
in renormalizable N = 2 SUSY only possible interactions are gauge interactions, the 
coupling between the hypermultiplet and the vector multiplets enters the Lagrangian 
in form of a cubic interaction, which is the same for all flavors [12]. In conclusion, 
the Lagrangian for minimally coupled hypermultiplets, can be described in the N = 1 
language as follow^] 


£ = 


d 2 0d 2 # ( Q\q v Qi + Qie v QJ J + / d 2 0 ( \/2 Qi®Qi + rmQiQi ) + c.c. , (3.89) 




where i = 1,..., Nf. 


s Notice that SU(2)r indices are suppressed. 
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To obtain the bosonic part of the five-dimensional SUSY Lagrangian for minimally 
interacting matter fields we again elevate the bosonic part of the above Lagrangian into 
five dimensions by considering real adjoint scalars. The result is 


^ matter — Tt[D m H x D m h{] +Tv[HlY a (a a ) XK H K 

- (SR a - H x M)(XH x - H X M )t 


(3.90) 


where a a are the Pauli matrices (3.3) and A, k = 1, 2 label fundamental representation of 


SU (2)ij. Here, the Higgs fields H x are assembled into N c x Nf matrices for convenience. 
Thus, combining all our results, the N = 1 SUSY SU(N C ) gauge theory with Nf 
hypermultiplets in (4+1) dimensions is defined by the Lagrangian 


Ai — £/SYM 4“ £>CS T &matter ■ 


3.2.2.3 A supersymmetric gauge theory for domain walls 


(3.91) 


We started our discussion of supersymmetry out of interest what kind of theories sup¬ 
ports BPS domain walls as their classical solutions. Let us, therefore, conclude it by 
a brief account, what conditions must be imposed on a generic model such as ( 3.91[ ) 
in order to ensure existence of a domain wall(s) solution. As we discussed in Ch. [2j 
both existence and type of a topological soliton is determined by a structure of the 
vacuum manifold. Domain walls are bound to appear, when there are discrete vacua. 
Supersymmetry, on the other hand, ensures that at least some of the domain wall so¬ 
lutions are BPS domain walls. Therefore, we only need to clarify what conditions lead 
to discrete vacua. 

There are only two of them. First, the gauge group must contain at least one 17(1) as 
a subgroup, which allows for non-zero FI term. And second, masses of hypermultiplets 
must not all be zero and they must be (at least partially [52]) non-degenerate. Existence 
of FI term and non-degenerate masses guarantee that some vacua will be discrete. For 
example, in the U(N C ) non-Abelian gauge theory coupled to Nf flavors Nf > N c , with 
non-degenerate masses and the FI terms, it has been shown be that the number of 
disconnected vacua is 


N f \ 


(N f -N C )\N C \ 

To illustrate this, let us write down the bosonic part of U(N c ) gauge theory coupled 
to Nf hypermultiplets with the FI term and a simple prepotential, 


T($) = ^W($ 2 ) • 


(3.92) 
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The normalization is chosen such that = -J=<F 0 l/v + <h a T a , where a runs from 

\/2,N c 

1,... ,N 2 — 1 and Tr[T a ,T^] = ^5 a * > . The values of U( 1) and SU(N C ) gauge couplings 
are taken to be the same for simplicity. The Lagrangian is given as 


L = — 


4g'- 


Z7' 


F 


aMN | ^ [ j \^a r^M\^a . ^ v" tt V a dvd 

+ ~2g2 UM2 -‘ U ^ + 2(p r a r a~ c a r a 


+ Tr 


+ ^iy a (a a ) AK i^ K - (EH X - - if A M)t] . (3.93) 


Here index a = 0,1,. .. N 2 — 1 labels generators of both 17(1) and SU(N C ) subgroups. 
We take the mass matrix to be totally degenerate M = diag(mi,..., mpr f ). 

After using SU(2)r rotation to fix FI parameters as c a = (0,0, c), with c > 0, we 
eliminate vector multiplet auxiliary fields Y a , obtaining the D-term potential in the 
form 


V D 



(cl Nc - H,H\ 


H 2 h\) 2 + AH\H 2 H\H\ 


(3.94) 


The complete potential is given as a sum of Vo and the T-term potential Vp 


V = V D + V F , V F = (Eff A -ff A M)(Eff A -ff A M) t . (3.95) 


Notice that if the mass matrix M ~ 1 N f , the above model possesses SU(Nf) global 
flavor symmetry of the Higgs fields H\ —> H\Upr f . Since we are considering the 
non-degenerate masses, this symmetry is explicitly broken down to a maximal torus 

H(l) iV /- 1 p| 

To find vacua of this model amounts to finding configurations of the Higgs fields 
H\ and the adjoint scalar field £ for which V = 0. The second term of Vp vanishes 
if we set one part of the Higgs doublet to be zero. Let us take H 2 = 0 and H\ = « 0 
The potential now reads 


V = \g 2 Tr 


(cl Nc - HH^f 


+ Tr 


(EH - MH)(EH - MH ) 


(3.96) 


Let us for clarity consider the simplest nontrivial case N c = 2 and Nf = 3. The first 
part of the potential is given in components as 


cl Nc - HH ] = 



1 0 
0 1 


|77oo| 2 + |i7oi| 2 + |i7o 2 | 2 -c 

IT* TT | TT* TT \ TT* IT 
^ 00^10 + ^ 01-^11 + ^ 02^12 


TT * 

H oo 

H* 10 \ 

IT* 

**01 

H* n 

TT* 

.02 

H*J 

Hoo 

H* w + 


hxH\ x + H Q2 H* 2 
\H\o\ 2 + \Hn\ 2 + \Hi 2 ( — c 


6 We subtracted overall (7(1) which is gauged and therefore not a genuine global symmetry. 

7 Tha fact that H 2 does not contribute to domain walls can be proven in general, see the appendix 


in |54|. 
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This will be zero if we choose for every row index of H, labeled by r, a color index A r . 
such that H r A r = \J~c. All other components in r-th row of H are set to zero. The 
phase of H r A r can be always eliminated by using flavor symmetry. Thus, in order to 
specify a vacuum, we pair up every flavor index A with a particular color index r. The 
list of all flavors {A\ ... An c } for each color represents a vacuum. In our example we 
can write down six possibilities {0,1}, {0, 2}, {1, 0}, {1, 2}, {2, 0}, {2,1}. However, the 
ordering of colors in the list does not really matter as the global gauge symmetry can 
always interchange any two colors. Thus, there are only three physically distinguishable 
combinations, say {0,1}, {0, 2}, {1, 2}. Moreover, for each vacuum the adjoint scalar is 
uniquely given as S = diag(mA 1 ,mA 2 ,'mA 3 )- This is the color-flavor locking. We can 
easily generalize this argument to a general case to prove the claim, that the number 
of discrete vacua are ^ Nf N ^'y N i • 
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4 

Localization of Abelian gauge 
fields on the domain wall 


In chapter [2] we have seen how scalar particles are localized on the domain wall within 
a simple toy model. Although, we did not address the question of localization of 
fermions, it is worth stressing, that a level of sophistication required to localize fermions 
is comparable with the one for bosons. For example, as we saw in Ch. [3j if we use 
supersymmetric theory in (4+l)-dimensions and BPS domain wall as a background 
solution, fermionic zero modes are naturally localized together with bosonic ones as 
a manifestation of unbroken supersymmetry in the (3+l)-dimensional effective theory. 
But even if the parent theory is not supersymmetric, there are natural mechanisms how 
to localize fermions, for example by considering appropriate Yukawa coupling terms 
between fermionic and scalar fields {3j. 

In this chapter, we will focus on the remaining ingredient of any realistic model 
of elementary particles: a gauge held. To be more precise, our goal is to introduce 
a mechanism, which localizes massless gauge fields on a domain wall. As we will see, 
this mechanism will be an essential ingredient for more realistic model building of 
effective theories, which we are going to undertake in the next chapter. But before we 
plunge ourselves into technical details, it is necessary to put this mechanism in the right 
context, by discussing main developments and setbacks of the gauge fields localization 
issue. 

Compared with scalar and fermionic fields, whose zero modes localize quite natu¬ 
rally, localization of gauge fields, turn out to be a difficult problem, which remained 
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unsolved for a very long time. If we set its beginning with the original proposal to use 
domain wall as a model for a brane by Rubakov and Shaposhnikov [3], we count 30 
years. During that time, there has been many proposals how to solve it (see [31J for 
a brief description of this development), but they either turn out to be inapplicable in 
all cases or very hard to implement. The solution, we are going to present, has been 
developed quite recently [31] and as we will show, do not suffer from any of these draw¬ 
backs. Calculations required for its implementation are actually quite simple and in 
principle, this mechanism can be used to localize Abelian and non-Abelian gauge fields 
on a domain wall in any dimensions, although its use is best justified in supersymmet¬ 
ric (4+l)-dimensional theories with eight supercharges. However, it is neither unique 
nor complete solution of the localization problem, as the discussion still continues. For 
example, it is not yet clear whether the localization mechanism is not destroyed by 
quantum corrections or whether it can be induced by them. We are not going to ad¬ 
dress these issues in this text and we will employ the mechanism in the form it has 
been presented in m- 

We will describe specifics of the gauge field localization problem in several steps. In 
the first section of this chapter we are going to attempt to localize an Abelian gauge field 
in a rather (with benefit of hindsight) naive way. We will investigate an effective action 
of a simple Abelian-Higgs model in (4 + l)-dimensions, which supports a BPS domain 
wall solution. We will show that this model has two degenerate vacua, both of which 
spontaneously breaks 17(1) gauge symmetry. Thus, we will consider a so-called Higgs 
phase (also a superconducting phase) in the bulk, where gauge particles are massive. 
However, it is well known that at the core of the domain wall, the local 17(1) symmetry 
is restored and thus it should support a massless gauge field, propagating along its 
(3+l)-dimensional world volume. This expectation, however, turns out to be incorrect. 
When we identify the physical reason for this setback, it will be easy to motivate a 
proposal made by Dvali and Shifman to overcome this problem. We will describe 
the so-called Dvali-Shifman mechanism in the second section. This will lead directly to 
the the localization mechanism introduced by Ohta and Sakai m, which is the content 
of the third and last section. 
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4.1 Abelian-Higgs model 

Let us consider a simple model with local U{ 1) symmetry in (4+l)-dimensions defined 


by the Lagrangian density 

£ = - ~^F MN F mn + ^d M ad M a + D M HD M H t - 1/ , (4.1) 

4 g- 2 g- 

2 

F = (|i/| 1 2 -u 2 ) 2 + |aF-iLM| 2 , (4.2) 

where er is a neutral scalar field and H = {Hr, Hr) is a row vector of charged Higgs 
fields. The covariant derivatives are given as 

D m H = d M H + i A m H (4.3) 

and Am is a 17(1) gauge field of the field strength 

F M n = 9mAn - 9nA m ■ (4.4) 


The mass matrix M is chosen in such a way to support a domain wall 

"=(o-°™)' (4 ' 5) 

with m > 0. The vacuum manifold consists of two degenerate vacua, with vacuum 
expectation values (VEVs) of various fields given as 

(a, H l , H r ) = (m, v, 0), (—m, 0, v ). (4.6) 


Apart from the local 17(1) gauge symmetry, there also exists a global U[1)a flavor 
symmetry^] which acts on the Higgs fields as Hr — > e ia H r, Hr -a- e~ ia HR. Both of 
these symmetries are spontaneously broken in both vacua. 

Notice that the same gauge coupling g appears not only in front of the kinetic terms 
for the gauge fields and the neutral scalar field a, but also in the quartic interaction 
term for the Higgs fields. This feature is motivated by supersymmetry. It can be 


shown, that the Lagrangian (4.1) can be embedded into a supersymmetric model with 
eight supercharges by adding appropriate bosonic and fermionic fields, which however, 
do not contribute to a domain wall solution. As shown in the previous chapter, this 


1 Phase rotation of Hl and Hr in the same direction is gauged and the remaining global symmetry 

is in the opposite direction. Thus the notation U(1)a- 
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property allows us to construct BPS domain wall solution, which will simplify much of 
the discussion below. Let us stress, however, that one may repeat all calculations with 
more generic coupling constant to arrive at essentially same conclusions. 

To find a domain wall solution, let y be the coordinate orthogonal to the domain wall 
and let all fields depend only on this coordinate. Also let A M = 0, where (i = 0,1, 2,3 
labels coordinates of (3+l)-dimensional world-volume of the domain wall. With this 
setup, the energy density £ = — L can be rewritten as 

£ = ^ (V) 2 + \ D v H \ 2 + y (Vl 2 - v 2 ) 2 + \<tH - HM\ 2 
= ^2 (V + 9 2 {\H\ 2 - v 2 )Y + | D y H + aH- HM\ 2 
+ v 2 d y a - dy(H(al-M)!!^ . (4.7) 


Since the second line in (4.7) is sum of squares, we see that £ is bounded from below 

(4.8) 


£ > v 2 d y a — d y (^H (crl — M)!!^ . 


This is the Bogomol’nyi bound we have encountered in chapter [2] Solution, that sat¬ 
urates the Bogomol’nyi bound, obeys a set of first order differential equations, called 
BPS equations 


d y a + g 2 (\H\ 2 -v 2 ) =0, (4.9) 

D y H + crH - HM = 0. (4.10) 


The BPS equations, which arise in the bosonic sector of supersymmetric gauge the¬ 
ories with eight supercharges, are most easily solved using the so-called moduli matrix 
method (see [56] for a review). In this method utilities the fact, that BPS equations are 


divisible into a vector multiplet part (in our case Eq. (4.9)) and a hypermultiplet part, 


Eq. (4.10). This nomenclature corresponds to types of multiplets in SUSY theories 


with eight supercharges. The hypermultiplet part can be always solved by introduction 
of a “Wilson line” 

S(?/) = Pe^ d! ' (ff+iA “ ) , (4.11) 

where P denotes a path ordering. One can easily verify that decomposition 


H = vS-'Hoe 


iMy 


(4.12) 
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solves the hypermultiplet BPS equation (4.10) identically for any values of Ho- In 
the case at hand, Ho is a complex vector. But in models with larger gauge group, 
this object is a matrix. Since it contains free parameters (also called moduli) of the 


solution, it is called a moduli matrix. The remaining equation (4.9) can be recast as 
the so-called master equation for the gauge invariant quantity Q = |S| 2 

1 


2y ; 


;dy(fl 1 d y Q) = v 2 {l — 14 1 flo) 


(4.13) 


where fio = Hoe 2My H^. Once the value of Ho is fixed, solving the master equation 
amounts to evaluation of Q, from which we can deduce (up to a phase) S and in turn 


H and o. Notice that, since the imaginary part of S never enters the equation (4.13), 
we cannot determine the value of A y . This ambiguity corresponds to our freedom to 
choose a gauge. In particular, we can always choose A y = 0. Since we want the domain 
wall solution as a background solution, we will always use this gauge. 


Also notice that decomposition (4.12) of H into the pair (S,Ho) is not unique. 
For example, a pair ( VS,VHq ) with V £ C/{0} gives the same H. In consequence, 
a physically distinguishable pairs (S', Ho) fall into equivalence classes induced by this 
transformation, called the D-transformation. This equivalence restricts the number of 
physically distinguishable parameters of the moduli matrix Ho. In our case, we can 
use D-equivalence to set one of the moduli of Ho to one, say Ho = (1, e 2myo e ia ), where 
we have rewritten the second moduli into the more convenient form with yo,a € M. 
It is easy to guess what these moduli represent. The yo is obviously a position of the 
domain wall in the y-direction and a is a phase difference of the left and right vacuum 
between which the domain wall interpolates. Using this notation, the master equation 


(4.13) can be recast as 

= v 2 (l - fl -1 [e 2my + e - 2my+4mM )) . 


(4.14) 


Despite its apparent simplicity, no analytic solution of this equation is known. There¬ 


fore, it must be solved numerically. The generic solution is shown in Fig 4.1 (left 
panel). 


In Fig. 4.1 (right panel) we show a comparison between cases with finite and in¬ 
finitely large gauge coupling g. We see that there is almost no qualitative difference 
between these two cases. In fact, as we argued in Ch. [2j the low-energy effective action 
is determined solely from a symmetry breaking pattern and should be insensitive to 
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Figure 4.1: The left panel shows profiles of Hr (solid line), Hr (long-dashed line), and 
a (dashed-line) with finite gauge coupling (g = 0.5). The right panel shows a plot of a: 
dashed curve for finite (g = 0.5) gauge coupling and solid curve for strong gauge coupling 
(g = oo). The other parameters are m = v = 1. 


value of g. Thus, at the level of approximation we are interested in, we can safely send 


g —> oo. In this limit the master equation (4.14) becomes an algebraic equation for 
with solution 

n = n 0 = e 2my + e - 2m y+ 4rn y °. ( 4 . 15 ) 


Using the axial gauge A y = 0 we have S = fl 1 / 2 and a = ^Q~ 1 d y ^. In terms of 
the original fields, we can summarize this domain wall solution in the infinite gauge 
coupling limit as 


Am = 0, 


er = mtanh(2?n(y — r/o)) , 
v 


H l = 
Hr = 


4- e _4m U-yo) ’ 

VQ -2m(y-y 0 ) e ia 
-)- Q-^y-yo) 


(4.16) 

(4.17) 

(4.18) 

(4.19) 


Having this explicit solution, we can now calculate the low energy effective La- 
grangian. For that we use the same procedure as in the chapter [2j namely the moduli 
approximation. We promote all moduli into the four dimensional fields 


Vo -4 2/o 


a 


a(x y '). 


(4.20) 


We plug the solution (4.16)-(4.19) with promoted moduli into the original Lagrangian 


(4.1) and take g —> oo. After performing these steps we select the only surviving two 
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derivative term = D^HD^H^ and integrate it over y to obtain the low-energy 
effective Lagrangian £ e ff 

OO OO 

£ efr = J dyL (2 '> = J d yD^HD^HK (4.21) 

— OO —OO 


Notice that here the gauge fields A^ are non-dynamical and should be eliminated 
r)£ ( 2 ) \ , \ 

_ = 0, =► A, = -^[Hd,H (4.22) 

Inserting this back into and introducing a matrix field !K = H^H, we arrive at the 
very concise form for the effective Lagrangian 


^eff 


2v 2 


dy Tr 




(4.23) 


which will prove useful in the future. After straightforward calculations we obtain the 
result 

£ e ff = 0 d^yod^yo + —d^ad^a. (4.24) 

2 4?n 

Let us make several observations about £ e ff. First notice that 2 mv 2 is precisely the 
tension of the domain wall 


OO 



— OO 


v 2 ct 


- H(<j1-M)H ] 


= 2rmr , 


(4.25) 


where we used Eq. (4.8) with the boundary values taken from Eqs. (4.17)-(4.19). Sec¬ 
ondly, interpretation of both terms in £ e ff is clear. The yo field represents a NG boson 
for spontaneously broken translation symmetry in the y direction, while a is NG bo¬ 
son corresponding to spontaneously broken U(1)a flavor symmetry. Most importantly, 
however, notice that there are no massless modes for gauge fields. On one hand, this is 
hardly surprising, since taking the g —> oo limit turned A M fields into auxiliary fields. 
On the other hand, as we already stressed in Ch. [2j the low-energy effective Lagrangian 
is determined from symmetry considerations alone and, as such, £ e g- should be insen¬ 
sitive to a particular value of g , including g —> oo. 

Despite these observations, one may not be convinced that the lack of gauge fields 


zero modes in (4.24) is a genuine result and not an artifact of the strong gauge coupling 
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limit. Let us, therefore, present a different argument, independent on this limit. We 


will consider the same model as in (4.1) but we will gauge the global U(1)a symmetry. 
Thus, we introduce new gauge fields, say Wm, and we add the corresponding kinetic 


term to the Lagrangian (4.1). We also modify the covariant derivatives of the Higgs 
fields as 


D m H = d M H + yA m H + i W m HQ , Q = 


1 0 
0 -1 


(4.26) 


where Q is a matrix of charges of H = (. Hl,Hr ) under U(1)a- It is easy to see, that 


the background solution (4.16)-(4.19) remains unchanged, with addition of Wm = 0 to 
the list. A derivation of the effective Lagrangian for this gauged Abelian-Higgs model 
follows the same pattern as in the ungauged case and we are not going to present it 
here (for detailed exposition see the Appendix A of |57| ). The result is 

,,2 y 2 

4 m' 


I 277771^ 71^ 

C eff = -—G^CT + ——d^yo + — (d^a + W„) (d^a + W") 


4e 2 


(4.27) 


where G ^ = d^Wy—dyW^. We see, that gauging the global U(1)a symmetry presented 
itself by the shift d^a + W^. If we expand the third term, however, we discover that 
the gauge fields W tl are massive, with the mass y/^". Therefore, they should not be 
considered as a part of the low-energy effective Lagrangian, which contains only zero 


modes. Thus, we returned to the effective Lagrangian (4.24). 

The qualitative reason, why the zero mode of a gauge field cannot be a part of the 
effective Lagrangian in either gauged or ungauged Abelian-Higgs model, is illustrated 


on Fig. |4.2[ We depicted there a brane surrounded by a bulk in the Higgs phase. To 
be concrete, let us identify the gauge symmetry of our model with electromagnetism. 
Then, we may say that the bulk is a superconductor. Due to a well-known Meissner 
effect, electric field lines must terminate on the surface of a superconductor. Thus, if we 
put test charges inside the brane, the force between them will not follow Coulomb law 
~ ~ 2 , but rather Yukawa law ~ e 2 . This is a signal that mediating particles, photons, 
are massive and the force is short-ranged. The mass of the photon is proportional to 
the inverse of the width of the brane. 

In conclusion, we verified that the naive approach, adopted in this section, does not 
lead to massless gauge fields localized on the domain wall. But we have also learned an 
important lesson, that we should not break the gauge symmetry in the bulk. Otherwise 
the Higgs phase outside the domain wall makes the gauge fields inside massive by an 
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brane 



Figure 4.2: Schematic picture of brane surrounded by a bulk in the Higgs (supercon¬ 
ducting) phase. We see that although 1/(1) gauge symmetry (here identified with elec¬ 
tromagnetism for illustration) is not broken inside the brane, due to the Meissner effect, 
the electric field lines are pulled into the bulk, effectively making the electromagnetic force 
short-ranged. The mass of the photon generated by this effect is proportional to the inverse 
of the width of the brane. 


analog of the Meissner effect from superconductivity. In the next section, we are going 


to show that a possible solution of this problem is to consider a dual version of Fig. 4.2 


where we replace a Higgs phase with a confining phase (see Fig. 4.3) 


4.2 Dvali-Shifman mechanism 


An idea, widely recognized as a plausible way how to localize gauge fields, was put 
forward in 1997 by Dvali and Shifman m- This idea is based on a duality consideration. 


That is, instead of a model with the Higgs phase in the bulk as in Fig. (4.2), we should 


rather consider a dual scenario, depicted on Fig. (4.3), where the bulk is in the confining 


phase. The reason is that now the dual analog of the Meissner effect pushes field lines 
into the brane, making the interaction between test charges exactly Coulomb-like. 

In their paper |55j . Dvali and Shifman considered an 517(2) gauge theory in (3+1)- 
dimensions as a toy model. The matter content of this model was arranged in such 
a way, that the domain wall solution produced a spontaneous breakdown of 517(2) 
symmetry down to 17(1) inside the wall. Outside the wall, the bulk remained invariant 
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brane 
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+ 


confining phase 

F ~w.. 


confining phase 
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Figure 4.3: Schematic picture of a brane surrounded by a confining phase. This time the 
dual analog of the Meissner effect pushes the electric field lines inside the brane, making 
the interaction of charges inside it long-ranged and exactly obeying the Coulomb law. 

under SU( 2) gauge symmetry, which is confining above some dynamical energy scale 
A. Dvali and Shifman showed, that if the relevant energy scales of the domain wall are 
much less than A, the massless U( 1) gauge field is trapped on the wall, since all states 
in the bulk are confined and, therefore, massive with the mass of order A. From a point 
of view of an observer inside the brane, it would take an energy exceeding A, to figure 
out that his/hers Universe is actually surrounded by a confining bulk. 

The major drawback of this idea is the fact that it is based on a non-perturbative 
phenomenon, confinement. The (4+1) dimensional Yang-Mills theory is non-renormalizable 
and whether it is confining in the same way as in four dimensions is really a conjec- 
ture0 Therefore, it is difficult to use the Dvali-Shifman mechanism for realistic model 
building of the brane world scenario without a certain amount of hand waving. Despite 
this technical complication, the Dvali-Shifman mechanism is still a sound candidate for 
such a model building (see for example [SS] and references therein). 

2 Despite many decades of research, a generally accepted proposition that four-dimensional Yang- 
Mills theory is confining remains conjectural too. 
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4.3 Ohta-Sakai mechanism 


A mechanism, which realizes confining phase of the bulk at the classical level, was 
introduced by Ohta and Sakai [31] in 2010. The advantage of their approach over 
the Dvali-Shifman mechanism is a perturbative implementation of the confining phase, 
allowing for an explicit model building. The crux of the Ohta-Sakai mechanism is 
best explained using the analogy from classical electromagnetism. In usual dielectric 
medium, the external electric field E is enhanced by induced polarization P, so that 
the total electric flux density D is given as 


D = £qE + P = eE , 


(4.28) 


where £q is the vacuum permeability. The dielectric permeability e is typically greater 
then £0) because typically an induced electric field points in the same direction as the 
external one. However, it has been recognized [59j. |60| that a classical representation 
of a vacuum in the confining phase could be described as a perfect dielectric medium 


e = 0. This idea suggests, that we can reinterpret the situation of Fig. (4.3) from the 
classical point of view as a continuous change of dielectric properties of the vacuum. 
The dielectric permeability e = e(y) varies along the extra-dimensional coordinate y 
in such a way that it rapidly decays away from the domain wall £ —> 0 as y —> Too, 
while staying nonzero near its center. A relativistic description of dielectric medium 
with such a non-trivial dependence on y is described by the Lagrangian m 


L = — - e{y)F MN F 


<MN 


(4.29) 


where Fmn = 9mAn — 9nAm- But this is nothing else than position-dependent gauge 
coupling constant 

e(y )=J F. (4.30) 

In their work m, Ohta and Sakai demonstrated that a spatially varying gauge coupling 
constant does localize massless modes of the gauge fields Am, independently on details 
of the profile g(y). 

In order to localize massless gauge fields on the domain wall in this way, it is 
necessary to introduce some field-theoretical device, which create an appropriate spatial 
profile of the gauge coupling g(y). Surprisingly, Ohta and Sakai showed that such a 
device is already present in a certain kind of theories and surprisingly still, these are 
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exactly those we are interested in! More specifically, they pointed out that generic 
properties of supersymmetric gauge theories in five dimensions naturally incorporates 


/le/d-dependent gauge couplings. Indeed, as we saw in Sec. (3.2) and in particular in 


Eq. (3.69), the kinetic term for gauge fields is multiplied by second derivatives of the 
prepotential. Since for five-dimensional SUSY theories the prepotential can be cubic 
polynomial SB, we see that an inverse of the square of the gauge coupling can depend 
linearly on the adjoint scalars fields E 


1 

si 


<9 2 /(E) = c“E a 


(4.31) 


where a is running through all generators of the gauge group. If we consider a domain 
wall solution and study small fluctuations in the background, adjoint scalars usually 
develops some nontrivial dependence on the y-coordinate. In this way, the gauge cou¬ 
pling ( 75 (E) becomes position-dependent through the dependence on the E(y) fields. 
Moreover, as the term mostly contributing to the low-energy effective Lagrangian is 
given as 

OO 

1 


d y 


-FnuF^ , 


, si(yy^ 

— OO 

the effective four-dimensional gauge coupling is given simply as 


(4.32) 


OO 

\ = [ dy 

91 J 


1 


sl(y) 


(4.33) 


All the remaining work now is to arrange things in such a way that l/y 2 (y) has ap¬ 


propriate y-dependence, such that the integral (4.33) converge to a positive number. If 
this condition is satisfied, we obtain localized massless gauge fields on the domain wall. 

The actual workings of the Ohta-Sakai mechanism is perhaps best illustrated on 
concrete examples. In the remaining of this chapter, let us discuss two particular 
models realizing Ohta-Sakai mechanism. 


4.3.1 Four-flavor model 

The so-called four-flavor model is technically the simplest model where the Ohta-Sakai 
mechanism is realized. The name points to the fact that there are four flavors of Higgs 
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fields. The Lagrangian is 
1 


T 4 . o T) 

*9? 


pMN 


- a(a)G M NG MN + + D M HiD M H\ - , 

*9i 




(4.34) 

+ l^i - HiM\ 2 , (4.35) 

where i = 1,2. Simply speaking, the four-flavor model is a double copy of the Abelian- 


V i = ^[H i H\-c i 


Higgs model we discussed in section (4.1), with additional 17(1) gauge symmetry with 
the field strength Gmn = 9mWn — Correspondingly, the notation we use here 

is almost the same as we used for single Abelian-Higgs model. Notice that the only 
interaction between both copies comes through the field-dependent gauge coupling-like 


term a (a). As we discussed in Sec. (3.2), this term represents second derivatives of the 
prepotential and, as such, it is a linear function of scalar fields a j. Since (apart from the 
a (a) factor) both sectors does not interact with each other, the domain wall solution 


of the Abelian-Higgs model (4.16)-(4.19) is a solution of the four-flavor model as well. 
More precisely, each sector allows this solution independently. If we consider the case, 
where there is a domain wall in both sectors, we have a t = mtanh( 2 m(y — yi)), where y* 
is the position of the i-th domain wall. Moreover, let us also choose the field-dependent 
gauge coupling to be 

a(a) = A(«ti - cr 2 ), (4.36) 

where A > 0. Following the discussion at the start of this section, it is easy to write 
down the effective Lagrangian for the four-flavor model 


where 


Left' — Left' ,1 T 7j e ff 2 ,o G^uG^ , 

4e| 


£eff= 1 ^ L d fl y i d^yi + -^d^oiid^ai 


(4.37) 


(4.38) 


is the effective Lagrangian of the i-th sector (compare with Eq. (4.24)). The effective 
four-dimensional gauge coupling e 4 is given by integration of position-dependent gauge 
coupling a(a(y)) over the extra-dimensional coordinate 


d y 


tanh 


lxj 

^2 = J d y a (<?(y)) = Xm J 

— 00 —00 

1 ( 2 m(y-y 2 )) = A(y 2 - yi). 


V 1 


— tanh( 


(4.39) 
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Notice that to obtain positive (inverse of the square of the) gauge coupling we must 
assume y 2 > y\. In effect, in this particular model, the Ohta-Sakai mechanism does 
not lead to the positive effective gauge coupling e\ in the entire moduli space. This 
feature may not be desired. The problem is that both y\ and y 2 are not parameters of 
the theory, but rather moduli subjected to dynamics. It may happen that during finite 
time of evolution 2/2 — yi hip sign, leading to the ill-defined kinetic term of localized 
gauge fields. This observation suggests that the four-flavor model might be unstable. 
Therefore, one would prefer to have the effective gauge coupling depend only on the 
parameters of the theory and not on moduli of the background solution. In the follow¬ 
ing, we will show that this can be achieved by considering slightly different model, the 
three-flavor model. 

4.3.2 Three-flavor model 

While the four-flavor model can be regarded as the most simple example of an Abelian 
gauge field localization, the three-flavor model, while achieving the same, is more eco¬ 
nomic. In the infinitely strong gauge coupling limit, which we take to obtain back¬ 
ground domain wall solution, each gauge group acts as a constraint on hypermultiplet 
scalars. The Ohta-Sakai mechanism is based on a gauge group containing a subgroup 
U{ 1) x JJ( 1). Thus, we have at least two constraints, making the minimum amount of 
Higgs fields to form a non-trivial background to be at least three. 

The Lagrangian of the three-flavor model reads 

£3 f = --T2^iMNFi IN — a{a)GMNG MN + —^9MOid M <Ji 
4 9t 2 df 

+ \D m H\\ 2 + | {D m ~ iA M )H 2 \ 2 + \(d M + iA M )H 3 \ 2 - V , ( 4 . 40 ) 

V = (<ti — m) 2 \Hi\ 2 + (<ti - a 2 ) 2 \H 2 \ 2 + o% \H 3 \ 2 , ( 4 . 41 ) 

where * = 1,2. Here the covariant derivatives are expanded as Dm = 9m +iWm + i Zm 
and the field strengths as F\mn = 9mWn — 9nWm , F 2 mn = 9mAn — 9mAm and 
Gmn = 9mZn ~ 9nZm ■ In the three-flavor model we choose a(cr) = \u 2 . 

We would like to verify, that this model leads to a positive four-dimensional effective 
gauge coupling, which does not depend on the moduli of the background solution. This 
amounts to find a domain wall solution and, in particular, calculate the profile of a 2 . 
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The BPS equations for the domain wall are (here we take Zm = 0) 


D y H\ + (cr — m)Hi = 0, 

(4.42) 

D y H 2 + (u — cr 2 — iA y )H 2 = 0 , 

(4.43) 

d y H 3 + (a 2 + \A y )H 3 = 0 , 

(4.44) 

\d y oi = Cl - \Hi \ 2 - I H 2 \ 2 , 

(4.45) 

91 


~^d y c t 2 = c 2 + \H 2 2 - \Hs\ 2 . 

9 2 

(4.46) 


It can be checked that a domain wall solution, satisfying above equations in the strong 
gauge coupling limit gf, g\ —> oo, is 


Am = Wm = Zm = 0, 

(4.47) 

Hi - ifc{ -_ , 

■y/l + e’l/eo 

(4.48) 

ia 

H 2 = s/ci -— , 

\/l + eo/e^ 

(4.49) 

H 3 = ^e~ y / 2 , 

(4.50) 

ai = m+ ^d y ln^l + e v /eo) , 

(4.51) 

02 = dyV , 

(4.52) 


where we have denoted eo = e 2m ^ y yo \ and where i] is a solution to the equatioij^] 

ci/c 2 


e~ v = 1 + 


(4.53) 


1 + e 0 e v 

It is clear from Eq. (4.52|) that the effective four-dimensional gauge coupling depends 


only on boundary values of g{y). These can be read off from Eq. (4.53) leading to the 


result 


^2 = f dya(a{y)) = A 


OO 

dya 2 (y) = X j dyd y r] 


= Aln( 1 + |). 


(4.54) 


Thus, in contrast with (4.39), we see that in the three-flavor model the effective four¬ 


dimensional gauge coupling depends only on parameters of the model. Interestingly, it 
can be shown that ln(l + C 1 /C 2 ) is proportional to the width of the domain wall. 

3 The correct root of the quadratic equation (4.531 is selected by demanding e _7) —> 1 as ci —> 0. 


69 















4. LOCALIZATION OF ABELIAN GAUGE FIELDS ON THE DOMAIN 
WALL 


4.3.3 Conclusion 


We have discussed two models as an illustration of the Ohta-Sakai mechanism, where a 
successful localization of Abelian gauge fields took place. In the four-flavor model the 
gauge fields are localized between two non-interacting domain walls and the strength of 
the effective gague coupling naturally depends on their mutual separation. However, as 


Eq. (4.39) shows, the ordering of the walls is important, pointing to a potential insta¬ 


bility of this model. In the three-flavor model, the result (4.54) shows the dependence 


of the effective gauge coupling on the width of the wall. In contrast to the four-flavor 
model, the four-dimensional gauge coupling is manifestly positive. Let us stress, how¬ 
ever, that these two models are particular examples of a whole class of models realizing 
the Otha-Sakai mechanism (see the discussion in [31]). In all these models, however, 
matter fields and gauge fields are localized without any mutual interaction. In the next 
chapter we will discuss how to improve on this feature. In particular, we will extend 
both the four-flavor and the three-flavor models to non-Abelian gauge models and we 
will show how to localize non-Abelian gauge fields together with minimally interacting 
matter fields on the walls. 
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5 

Non-Abelian matter fields on the 
domain wall 


In this chapter we study a non-Abelian extension of the four-flavor model of Ch. [4j 
which we call the chiral model. In this model, two new features arise in the effective 
Lagrangian, which have not so far appeared in any other model presented in this text. 
First, massless non-Abelian gauge fields are localized on the domain wall and second, 
the localized scalar fields are non-trivially interacting under localized non-Abelian gauge 
symmetry. 

Organization of this chapter is as follows. In the first section we discuss the model 
together with the background domain wall solution. Having large gauge and flavor 
symmetry groups and degenerate masses, the chiral model posses multiple discrete 
vacua and multiple domain walls. As we want to preserve most of the flavor symmetry, 
we will take two most symmetric vacua and a very special domain wall solution, which 
leaves most of the flavor symmetry unbroken. Such special solution is called coincident 
domain wall solution, since it can be viewed as composed of multiple, indistinguishable 
domain walls, sitting at the same position. 

In the second section we gauge the unbroken flavor symmetry and show that local¬ 
ized scalar fields interact minimally with localized gauge fields. Both in this chapter 
and in the following chapter, we will pay a lot of attention to the derivation of the 
effective Lagrangian, since the precise character of interactions between localized zero 
modes has phenomenological significance. We devote the third section to these matters. 
As a conclusion, in the fourth section we discuss the geometrical Higgs mechanism, a 
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SU(N) C 

U( 1)1 

U( 1)2 

SU(N) l 

SU(N)r 

U(\) 1A 

U(1) 2A 

mass 

H\l 

□ 

1 

0 

□ 

1 

1 

0 

mil TV 

H\r 

□ 

1 

0 

1 

□ 

-1 

0 

-mil N 

Hi 

adj © 1 

0 

0 

1 

1 

0 

0 

0 

h 2L 

1 

0 

1 

1 

1 

0 

1 

m 2 

H 2 r 

1 

0 

1 

1 

1 

0 

-1 

—m 2 

h 2 

1 

0 

0 

1 

1 

0 

0 

0 


Table 5.1: Quantum numbers of the domain wall sectors in the chiral model. Here □ 
denotes that a particular field belongs to the fundamental representation of the given group, 
similarly 1 denotes singlets and adj denotes adjoint representation. 


unique feature of the effective theory, resembling the similar mechanism of D3-branes 
from string theory. 

Most of this chapter, content and notation alike, is adopted from the work m- 


5.1 The domain walls in the chiral model 


As a natural extension of the four-flavor model of Ch. [4j we consider the Yang-Mills- 
Higgs model with SU (N) c xU(l) gauge symmetry and S[U ( N ) l x U (N)r] = SU ( N) l x 
SU(N)r x U(1)a flavor symmetry. In order to localize the gauge fields, we again 
introduce two sectors £i and L 2 , but only the former is extended to Yang-Mills-Higgs 


system and the latter is the same form as (4.1). The second sector couples to the first 


sector through the field dependent coupling of Eq. (4.36), which after gauging the flavor 
symmetry, plays a role of position dependent gauge coupling, localizing massless modes 


of gauge fields on the domain wall. The matter content is summarized in Tab. 5.1 Since 
the presence of two factors of SU ( N ) global symmetry resembles the chiral symmetry 
of QCD, we call this Yang-Mills-Higgs system the chiral model. 

The Lagrangian is given by 


£ = £1 

£1 = Tr 
Hi = Tr 


+ To , 

-77-2 (HiAfAT) 2 + “2 (Dm'Ei ) 2 + \D M Hi \ 2 
L 2 9l 91 

j (HiH\ - u 2 Tv ) 2 + iSrHi - Hi Ah ! 2 


Vi, 


(5.1) 

(5.2) 

(5.3) 
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with H\ = (Hu, H\ R ) and L 2 being equal to the Lagrangian of Eq. (4.1). Gauge fields 
of U(N) C = ( SU(N) C x U( 1)i)/Zn are denoted as W\m, and adjoint scalar as £ 1 . The 
covariant derivatives and the field strength are denoted as 


Dm'Z 1 = 0Af£i + i [Wi m, £ 1 ] , (5.4) 

DmHi = OmHi + \W\mH] , (5.5) 


and Fimn = 9mW±n — djyWiM + i \W\m-, Wuv]- The mass matrix is given by M\ = 
diag (miljvj — miljv)- Let us note that the chiral model reduces to the (ungauged) 
four-flavor model in the limit of N —> 1, by deleting all of the SU(N) groups. The 
second sector is needed to realize the field-dependent gauge coupling function, which 
we will discuss in the subsequent section. In the rest of this section, we focus only on 
the first sector (i = 1) and suppress the index i. The symmetry transformations act on 
the fields as 

H = (H L ,H R ) U c (H Ll H R )( ULel ° URe -ia) , (5-6) 

£ ->■ UcFUl , (5.7) 


with U c € U(N ) C , U L € SU(N) l , U(N) r € SU{N) R and e ia € U(1) A . 
There exist N + 1 vacua in which the fields develop the following VEV 


H = (H L ,H R )=v^ lN ~ r 

s ™ -c ) ■ 


(5.8) 

(5.9) 


with r = 0,1,2, ■ ■ ■ , N. In the r-th vacuum, both the local gauge symmetry U(N) C 
and the global symmetry are broken, but a diagonal global symmetries are unbroken 
(color-flavor-locking) 


U(N) C x SU(N) l x SU(N)r x U(1) A ->• 

SU(N - r) L+c x SU(r) L x SU(r) R+c x SU(N - r) R x U(l) A+ c • (5.10) 


We obtain the BPS equations through the Bogomol’nyi completion of the energy 
density with the assumption that all the fields depend on only the fifth coordinate y 
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and = 0: 


£ = TV 


-2 \D y E- 9 — (v 2 l N -HH^j +\D y H + T 1 H-HMY 


+ d y |tl- v 2 E - (EH - HM) H ] } 
(EH - HM) H ] } . 


> 3, {tv 


u 2 S- 


This bound is saturated when the following BPS equations are satisfied 

2 

D y E- 9 -(v 2 l N -HH^ =0, 

D y H + EH - HM = 0. 

The tension of the domain wall is given by 


(5.11) 

(5.12) 

(5.13) 


T = 


dy d y {TV v 2 E — (EH — HM) H^ } 


= v 2 TV [E(Too) — E( — oo)] . 


(5.14) 


We will concentrate on the domain wall which connects the 0-th vacuum at y —> oo 
and the IV-th vacuum at y —> —oo. Its tension is 


Trnin. = 2 Nv m . 


(5.15) 


which follows from Eq. (5.14). Notice that T co { n is N times the tension of single domain 


wall of the Abelian-Higgs model (4.25). Thus, we interpret such a domain wall as being 


composed of N different (but indistinguishable) domain walls. In the simplest case, all 
domain walls have the same position, which corresponds to making an ansatz that Hl, 
Hr, E and W y are all proportional to the unit matrix. Then the BPS equations (5.12) 


and (5.13) can be identified with the BPS equations in the Abelian-Higgs model. The 


coincident domain wall solution is 


H l = ve~ 2 e my 1 


N : 


Hr = ve 2 e my 1 tv ; 
E + i Wy = ljv 


where ?/> is the solution of the master equation 

-Ls 2 v> = v 2 (i- e -^e 2m ^y-y°^ 


(5.16) 

(5.17) 

(5.18) 

(5.19) 
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Eq. (5.10) shows that the unbroken global symmetry for IV-th vacuum (Hr = 0, Hr = 
vl n and £ = —mljv) at the left infinity y — > — oo is SU(N)i x SU(N)r +c x U(1)a+c, 
whereas that for the 0-th vacuum (Hi = vIn, Hr = 0 and £ = mljv) at the right 
infinity y ->• oo is SU(N) L+C x SU(N)r x U(1)a+ c - 

The domain wall solution further breaks these unbroken symmetries because it 
interpolates the two vacua. The breaking pattern by the domain wall is 


U(N) C x SU(N) l x SU(N) r x U( 1) A -> SU(N) L+R+C . 


(5.20) 


This spontaneous breaking of the global symmetry gives NG modes on the domain wall 
as massless degrees of freedom valued on the coset similarly to the chiral symmetry 
breaking in QCB 0 


SU(N) l x SU(N) r 


xU(l) 


A- 


(5.21) 


SU(N) l+r 

Since our model can be embedded into a supersymmetric field theory, these NG modes 
( U(N ) chiral fields) appear as complex scalar fields accompanied with additional N 2 
pseudo-NG modes 0 


5.1.1 Localization of the matter fields 

In the following, we will present the low-energy effective Lagrangian, where the mass¬ 
less moduli fields (the matter fields) are localized. This serves the double purpose of 
establishing our notation and determining the spectrum of massless modes of matter 
fields on the coincident domain wall. 

We employ the moduli matrix formalism 

H l = ve my S~ l , (5.22) 

Hr = ve-wS-'e* , (5.23) 

£ + iWy = S-'dyS , (5.24) 

where S E GL(N , C) and Q = SS t is the solution of the following master equation 

~2 d y (0,~ 1 d y Q) = v 2 (In — n _1 rio) j (5.25) 

1 Local gauge symmetry SU(N) C contains global symmetry as a constant gauge transformation, 
which is displayed in the above symmetry breaking pattern. However, Nambu-Goldstone modes only 
come from the genuine global symmetry which is not locally gauged. Hence we do not count the 
SU(N) C transformations. 

2 One of them is actually a genuine NG mode corresponding to the broken translation. 
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where 

n 0 = e 2my 1 N + e~ 2my e^ f . 

We have used the E-transformation to identify the moduli e^, which is a complex N 
by N matrix. It can be parametrized by an N x N hermitian matrix x and a unitary 
matrix U as 


= e*U*, 


(5.26) 


where U is nothing but the U ( N ) chiral fields associated with the spontaneous symme¬ 


try breaking Eq. (5.21) and components of x are the pseudo-NG modes whose existence 


we promised above. 

In the strong gauge coupling limit g —> oo, solution of master equation is simply 
II = Hq- After fixing the U(N) C gauge, we obtain 


S = e x / 2 \j2 cosh(2my — x) 


Let us denote, for brevity 


y = 2 my — x , 


the Higgs fields are then given as 


H l =v 


Hr = V- 


e y / 2 

\/2 coshy ’ 

e~y / 2 

--U*. 


(5.27) 

(5.28) 

(5.29) 

(5.30) 


y/2 coshy 

From this solution, one can easily recognize that eigenvalues of x correspond to the 
positions of the N domain walls in the y direction. Now we promote moduli parameters 
x and U to fields on the domain wall world volume, namely functions of world volume 
coordinates x E We plug the domain wall solutions RL,i?(y; x(x^), U{x y ’)) into the 
original Lagrangian L in Eq.(5.2) at g —> oo and pick up the terms quadratic in the 


derivatives. Thus the low energy effective Lagrangian is given by 


^eff 


dyTr 


d iL H L &‘H\ + d.HRcHH^ - v 2 W^W y 


(5.31) 


where 




+{L»R) 


t 


(5.32) 
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5.2 Localization of gauge fields 


Here we have eliminated the massive gauge field W^ by using the equation of motion. 
Using the solutions for Hr and Hr it is possible to find a closed formula for the 
effective Lagrangian up to the second order of derivatives with full nonlinear interactions 


involving moduli fields x and U. The result is the same as in Eq. (5.46), only with 
covariant deivatives replaced with partial derivatives. Detailed derivation is given in 
the third section of this chapter. 

Here we show the result only in the leading orders of U — In and x: 

L eS = —Ti (d,,^ 8^U + d tl id ll xj +... (5.33) 

When N = 1 and with the redefinitions U = e 1Q , and x = 2 my, this coincides with the 


effective Lagrangian of Eq. (4.24) of the Abelian-Higgs model. As already discussed, U 
is N x N unitary field, which we identify as NG bosons of spontaneously broken flavor 
symmetry and x is IV x N Hermitian matrix of pseudo-NG bosons. We can identfy the 
“pion” field of our model usign the ansatz 


= i Ud^JJ ] , 

J 7T 


(5.34) 


where the pion decay constant f n = \J The effective Lagrangian now reads 

2 

Tr (d^d^T^j + ^Ti-^dfjxd^x'j + ... (5.35) 


5.2 Localization of gauge fields 


Let us next introduce gauge fields which are to be localized on the domain walls. As we 
learned in the previous chapter, the associated gauge symmetry should not be broken 
by the domain walls. Therefore, the symmetry which we can gauge is the unbroken 
symmetry SU (N) l + r itself. 

Let us gauge SU(N)l + r = SU(N)v and let be the SU(N)y gauge fields. The 
Higgs fields are in the bi-fundamental representation of U(N) C and SU(N)y- The 
covariant derivatives of the Higgs fields are modified as 


£>mHh = 8mH\ L + \W\mH\l - \Hi L A m , (5.36) 

DmHir = OmH\r + \W\mH\r - \HirA m ■ (5.37) 


The quantum numbers are summarized in Tab. 


5.2 
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h 2R 

1 

0 

1 

1 

0 

-1 

—m2 
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1 

0 

0 

1 

0 

0 

0 


Table 5.2: Quantum numbers of the domain wall sectors in gauged chiral model. The 
notation is explained in Tab. 


5.1 


We now introduce a field-dependent gauge coupling function, which is a direct 
analog to that of the four-flavor model. 


_1__ A /Tr(Si) _ Es 

2e 2 (E) 2 V Nmi m 2 


The Lagrangian is given by 

£ = £1 + £2 — 


2e 2 (E) 


Tr [G_ 


MN 


qMNI 


(5.38) 


(5.39) 


The £1 in Eq. (5.39) is given by Eq. (5.2) with the covariant derivatives replaced with 


those in Eqs. (5.36) and (5.37). 


First we wish to find the domain wall solutions in this extended model. As before, 
we make ansatz that all the fields depend on only y and = A^ = 0. Let us first 
look at the equation of motion of the new gauge field Am, which is of the form 


D m G mn =J 


t N 


(5.40) 


where Jm stands for the current of Am- Note that the current Jm is zero, by definition, 
if we plug the domain wall solution in the chiral model before gauging the SU(N) R+R . 
This is because the domain wall configurations do not break SU(N)l+r- Therefore, 


Am = 0 is a solution of Eq. (5.40). 


Then, we are left with equations of motion with Am = 0, which are identical to 
those in the ungauged chiral model in the previous sections. Therefore, the gauged 


chiral model admits the same domain wall solutions as those (Eqs. (5.29) and (5.30)) 
in the ungauged chiral model. 
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5.2 Localization of gauge fields 


The next step is to derive the low energy effective theory on the domain wall world- 
volume in the moduli approximation as in the previous subsection. Again, we promote 
the moduli parameters as the fields on the domain wall world-volume and pick up the 


terms up to the quadratic order of the derivative d^. Similarly to subsection 5.1.1 
we utilize the strong gauge coupling limit <?* —> oo. Let us emphasize that we keep 
the field-dependent gauge coupling function e(£) finite. The spectrum of massless NG 
modes is unchanged by switching on the SU(N)l + r gauge interactions^] 

We just repeat the similar computation to those in subsection 5.1.1| Again we shall 
focus on the first sector £1 and suppress the index i = 1 of fields. Since color gauge 
fields becomes auxiliary fields and eliminated through their equations of motion, 
it is convenient to define the covariant derivative only for the flavor SU{N)l+r gauge 
interactions as 

(5.41) 


D„H = dfjH - \HA^ . 


Then we obtain the effective Lagrangian of the first sector as 


L 


i,eff = J dy TV 

— OO 

1 


t + b»H R {b*H R )t - v 2 W^W^ 


2e 2 (£) 


GmnG 


MN 


with 


^=2u 2 L 


D^HlHI - H L (b„H L )t + (L^R) 


(5.42) 


(5.43) 


Eliminating VF /t , we obtain the following expression for the integrand of the effective 
Lagrangian after some simplifications 


£eff “ 2v 2 


dy Tr 


V„H ab WH b0 


(5.44) 


where we defined fields H ab with the label ab of the adjoint representation of the flavor 
gauge group SU(N)l+r+ c ■ The covariant derivatives are given as 


r ^iiHab = djiH a b + i[T ;j ,, H ah ] , H ab = H\H b , a,b = L,R. 


(5.45) 


3 Tree level mass spectra are unchanged even though the chiral symmetry SU(N)l x SU(N)r is 
broken by the SU(N)l+r gauge interactions. 
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5. NON-ABELIAN MATTER FIELDS ON THE DOMAIN WALL 


We will describe the procedure to derive the effective Lagrangian fully in the next 
section. Here we merely state the result: 


where 


£ieff = -Tr 

’ 2 m 


' cosh(£f) — 1 
. " -Cl sinh(^) 


+ U^ U U 


In 


1 + tanh(£i) 

1 — tanh(£.x) 

1 + tanh(£j) 


(D^s) 


cosh(Tf) — 1 
Lx sinh(£x) V. 1 — tanh(£ j) 


(T^x) 


+ \uW \ - inf 1 + tM y > i V J y g ) 


tanh(Tf) \1 — tanh(£.j.) J 


L A {B) = [A, B\ 


(5.46) 


(5.47) 


is a Lie derivative with respect to A. The covariant derivative T> fl is defined by 


V^U = d^u + i [Ap, U ]. 


(5.48) 


The above result suggests that the chiral fields U(x 11 ) and hermitian fields x(x At ) 
are in the adjoint representation of SU{N)l + r. We will prove this assertion in the last 
section of this chapter. 


5.3 Effective Lagrangian on the domain wall 


In this section we derive our main result (5.46) of the effective Lagrangian for the 
gauged chiral model. 


5.3.1 Compact form of gauged nonlinear model 


Our starting point is the Lagrangian (5.42). We adopt the Einstein summation conven¬ 


tion for a = {L, R} and, for brevity, we ignore the kinetic term for the localized gauge 
fields An- Thus, we have 


^eff 


d y Tr 


- v 2 WnW» 


(5.49) 


with the covariant derivatives given in Eq. (5.41) and with the Higgs fields obeying the 
constraint 

Halil = v 2 1n , (5.50) 
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5.3 Effective Lagrangian on the domain wall 


coming from the strong gauge coupling limit. We see that this model is actually a 
nonlinear sigma model. Let us first eliminate the gauge fields W^ to make this fact 
manifest. Using equation of motion we arrive at 




D,H a H\-H a b,Hl 


(5.51) 


and 


D^H = d^H - \HA tl . 


(5.52) 


The effective Lagrangian (5.49) can be simplified by using the following identities 


HaDpHl = d^H a Hl) - D.HaHl , (5.53) 

= -D^HlHb + T)^H ab , (5.54) 


where 

’D^Hab = d^Hab + i [A^, H ab } , H ab = H\H b . 
After some algebra we find: 

WnW* = ^DftHaD^Hl - ^DftHobWHha . 


Plugging above expression back into the (J5.49 ) we conclude 

VftH ab V^H ba 


^eff 


i r°° 

2^L dyTj 


(5.55) 


(5.56) 


5.3.2 Effective Lagrangian 

Now we are ready to perform the integration over the extra-dimensional coordinate. 
The background solution is given as: 


where y 


H = {H l ,H r ) = 


v e^/ 2 v e \ 

\/2 i/cosh (y)' y/2 v /cosh(y)J ’ 


myljsr — x. The composite fields H ab = HlH b are easily shown to be 


Hll = 


Hlr = 


v 2 


2 cosh(y) ’ 


2 cosh(y) ^ Hrl ’ 


v 2 e y + 

h rr = -u—— ub 

2 cosh(y) 


(5.57) 


(5.58) 

(5.59) 

(5.60) 
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5. NON-ABELIAN MATTER FIELDS ON THE DOMAIN WALL 


Using these expansions we arrive at the form 


- 


r°° ( eV 

/ dyTr^—+ 
J —< 


+U^v u u 


' —oo 

*-y 


\ M cosh(y) cosh(y) 


1P- 


cosh(y) 


,U ]r D^(jJ 


y 


cosh(y)/ 




+ UH\U 

1 


cosh(y) cosh(y) 

_ 1 _, _ 1 

cosh(-y)’ cosh(y) 


costr(y) 


(5.61) 


In the following, we would like to carry out the integration over the extra-dimensional 
coordinate y. This can be done in two steps. First, we must factorize all quantities 
depending on y (or on y) to one term inside the trace, effectively reducing our problem 
to match the general form 


dy Tr 


f(my 1 N - x)M 


(5.62) 


where M is some matrix, independent of y and / is some function. In the second step, 
we diagonalize x 

x = P _1 diag(Ai,..., Xn)P 

and use the fact that f(P~ l yP) = P _1 f(y)P. This transformation leads to 

/°° r i f‘°° ^ 

dy Tr f(myl N - diag(A*)) PMP -1 = / dy ^ /(my - \ l )(PMP ~ l ) ll . 

-OO L J J—OO i=l 

For every term in the sum we can perform a substitution y = my — A*. The key 
observation is that in each term the integration is independent on a particular value of 
A j. This allow us to conclude 


f 


1 


dy Tr f(y)M = —Tr(M) / dy/(y). 


m 


(5.63) 


It appears as if we have just made a substitution y = yljy. This is possible, of course, 
only thanks to the diagonalization trick and properties of the trace. In the subsequent 
subsections, however, we will refer to this procedure as if it is just a ‘substitution’, for 
brevity. 


Let us decompose the effective Lagrangian (5.61) into three pieces 

^eff T ~\~ 7 mixed 

and see the outlined procedure for each term. 


(5.64) 
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5.3 Effective Lagrangian on the domain wall 


5.3.3 Kinetic term for U 


First, let us concentrate only on terms containing double derivatives of U, which we 
denote 7jj: 


v 2 f°° ( 1 


cosh(y) 


,U^V^U 


2,-y 


cos. 


Hv) } 


where we have used a fact that inside the commutator it is possible to freely interchange 

P -y P y 


cosh(y) cosh(y) ’ 

since the difference is just a constant matrix. In this way we made T jj manifestly 
invariant under exchange y —> — y. 

Since in the first factor of all y-dependent quantities are on the right side, we 


can, according to our previous discussion, make use of the identity (5.63) and carry out 
the integration: 


v 2 m 

T dy Tr 

* J — OO 


V u U ]c D^U- 


1 


cosir (y) i 2m 


v m 
= —TV 




For the second term, however, we first use the identity: 




k\ 

k=1 

where La{B) = [d.. B\ is a Lie derivative with respect to A. Thus 

\k 

k\ 


(5.65) 


cosh(y) 




/ p y \( k ) 

E L_L £ | (c ,» d , e/) (_) 


Now all y-dependent factors are standing on the right and we can formally exchange 
y —> y. The summation can be carried out easily 


y f eV ^ (fc) 

^ k\ x( j Vcos h(y)J 


The formula for 7\j now reads: 


qV ^x qV 


dsTr 


cosh(y - L £ ) 




(U^V^U) - 


cosh(y) 


U ] D IJ V. 


cosh(y — Lx) cosh(y) 


(U ] T> tM U) r D ll U^U 


(5.66) 


(5.67) 
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5. NON-ABELIAN MATTER FIELDS ON THE DOMAIN WALL 


Since we started with 7jj invariant under the transformation y —>• —y, we should take 
only even part of the above formula (with respect to the exchange —> —£x) as the 
final result 


7u = ^i 


f 


dy Tr 


cosh(£ ;i 


_cosh(y - £ £ )cosh(y) 




(5.68) 


Using the primitive function 
f d y 


1 


In 


1 


J cosh(y — a) cosh(y) sinh(a) 1 — tanh(a) tanh(y) 

we obtain the result to all orders in x 

„,2 


= i^ Tr 


' Inf 1 + tan fj^! ) (U'V“U) 
M tanh(T^) \ 1 — tanh(£.£) J 


(5.69) 


Performing the Taylor-expansion of the function 


1 , /l + tanh(x)\ _ 2x 2 2x 4 4x 6 2x 8 4x 10 12 . . 

tanh(x) n (l- tanh(x) J _ ^ + ^ __ l5" + 945 _ 4725 + 93555 + ^ ^ ^ ^ 


we can easily read off coefficients of terms beyond the leading one. For example, the 
first three terms reads: 






5.3.4 Mixed term 


v m 

-Ti’ 

6 m 

v 2 „ 


90m 


( x,U ]r DJJ x, D^U^U ) 

•( x, x, )+... (5.71) 


Mixed term between x and U is given by 


,,2 poo 


7mixed=-J dy Tr<j U^Ti^U 


1 1 


cosh(y) cosh(y)J [. c °sli(y) cosh(y) 






With use of the identity (5.65) and 


D m /(x) = £t|(T a 


fc =0 


/ (fc+1) (®) 
(A: + 1)! ’ 


(5.72) 


one can prove the following: 


n =2 


[/(x), X> M y(x)] = 1 (D ll x) (/(x)y(x)j - f (n \x)g(x) - f(x)g {n \x) 


(n) 
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5.3 Effective Lagrangian on the domain wall 


We can use this result to factorize all y-dependent quantities to the right and make the 
substitution y = y In 


,2 ro o 


7, 


mixed 


d v J2 


(-i y 


4m J -oo — n\ 

W e -ii 


Tr 


+ 


e~y \ (n) & 


cosh(y) J cosh(y) ' ^cosh(y )) cosh(y) \ v cosh(y) J cosh(y) 

Now we are free to perform summation and integration to obtain 


- 2 


1 


(n) 


1 


7mixed ., Tr 
2m 


u'n„u 1" ( 1 + tsn ^b (m 

A Lx smh(£j.) \_ 1 — tanh(£^) / 


(5.73) 


Performing the Taylor-expansion of the function 


cosh(.x) — 1 , /1 + tanh(x) \ x 3 x 5 I7x 7 31.t 9 n , 

x sinh(x) n V1 - tanh(x )) “ X _ 12 + 120 _ 20160 + 362880 + ^ ^ 

we can easily read off coefficients of terms beyond the leading order in the expansion 


7mixed .. Tl 
2m 


U^V^U[x,V^x] 


24 m 


-Tr 


[x, [x, [x, X ]]] + ... (5.75) 


5.3.5 Kinetic term for x 

Kinetic term for x is given by 

,,2 r°° ( \ 


T~ = — 


- / dyTr<;^—— a)"—L- 

4 J -oo l cosh(y) cosh(y) cosh(y) cosh(y)J 


We are going to need the identity 


Tr 


i 00 (— 1 ) n 

’D,j(x)vy g (x) = J2 


n =2 


(n) 


x)g(x)) - f (n) (x)g(x) - f(x)g( n) (x) } . (5.77) 


With the aid of this we arrive at 


(5.76) 


„2 roc 


7x = 


dy X^^Tr^r 2 ^)^ 


4m J -oo - ,^ 2 n\ 


W e -y 


, ._. ~_ | ( e~y \ (n) e y / 1 \ (n) 1 ' 

\cosh(y )J cosh(y) \cosh(y )J cosh(y) \cosh(y )J cosh(y) 
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where we have again employed the diagonalization trick and the identity (5.63). After 
performing the summation and the integration we obtain 

,2 r cosh(£f) — 1 /1 + tanh(£f)' 


T., = 


2m 


-Tr 


^ fiX r 2 ■ i , ,. 

smh(£ ; i 


In 


1 — tanh(£f) 


(Wx) 


leading to the power series 
v 2 

7$ = -Tr 

2 m 




v m 

H-Tr 

24m 


[x, I> M x] [x, V tt x\ 


+ ... 


(5.78) 


(5.79) 


Putting all pieces together as £ e ff = + “^u + 7 m i xe d, we confirmed our final result 


(5.46). 


5.4 Discussion 

Let us illustrate nonlinear interactions in Ti ie ff up to fourth orders in the fluctuations 
x and U — In 

v 2 / 

£i,eff = — TrlVprtWU + VpxWx + rfVpUlx^xl 


- ^ [D„x, x} [x, T^x] + ^ x] [x, C/tT^C/] + • 


(5.80) 


where "D^U = d^U + i [A^, U] and T)^x = d^x + i [A^, x\. Our result suggests that U 
and x are in the adjoint representation of the SU(N)l +r flavor gauge group. Let us 
now examine the transformation properties of U and x in order to demonstrate that 
they are, indeed, in the adjoint representation. The domain wall solution only preserves 


the diagonal subgroup SU(N)l.\-r+ c . Eqs. (5.6) and (5.7) shows the fields transform 
under the SU(N)l +r+c transformations U as 


H' l = UH l lt f , H' r = UH r U ] , S' = UEU t . 


(5.81) 


Eqs. (5.22) and (5.23) show that 


S' = 11511*, e^ = He 0 U* , fl' = UC2U f . (5.82) 

The complex moduli e^ is decomposed into hermitian part e x and unitary part U in 


Eq. (5.26). Since we can express e 2x = eSe ^, and U = e ^e x , we hnd that they 


transform as the adjoint representation 


2x 


' = He 24 !!*, U' = UUU ] 


(5.83) 
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5.4 Discussion 


Similarly to Eq. (4.39), we can define the (3+l)-dimensional non-Abelian gauge 


coupling e 4 by integrating (5.38) and find 


2 ?= ldy 


1 


= Kv 2 - yi) 


(5.84) 


, 4 ^ 2e2(E) 

where m is the wall position for the i-th domain wall sector. Summarizing, we have 
obtained the following effective Lagrangian 


^eff — ^l,eft T ^2,eff „ 2^^ 

2e 4 


G^G' 11 ' 


(5.85) 


where £ 2 ,eff is given in Eq. 4.24 This is the main result of this chapter. We have 
succeeded in constructing the low-energy effective theory in which the matter fields 
(the chiral fields) and the non-Abelian gauge fields are localized with the non-trivial 
interaction. We show the profile of ’’wave functions” of localized massless gauge field 
and massless matter fields as functions of the coordinate y of the extra dimension in 


Fig. 5.1 



Figure 5.1: The wave functions of the zero modes. DW1 and DW2 stand for the wave 
functions of the massless matter fields of the i = 1 domain wall and * = 2 domain wall, 
respectively for strong gauge coupling limit gi = oo and m 7 ; = 1. The gauge fields are 
localized between the domain walls. 


As we see from Eq. (5.80), the flavor gauge symmetry SU(N)l + r +c is further 
(partly) broken and the corresponding gauge field A^ becomes massive, when the fluc¬ 
tuation cj) = e x U develops non-zero vacuum expectation values. Especially, x is inter¬ 
esting because its non-vanishing (diagonal) values of the fluctuation has the physical 
meaning as the separation between walls away from the coincident case. For instance, 
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5. NON-ABELIAN MATTER FIELDS ON THE DOMAIN WALL 


if all the walls are separated, SU(N)l + r +c is spontaneously broken to the maximal 
17(1) subgroup 17(l) Ar ~ 1 . However, if r walls are still coincident and all other walls are 
separated, we have an unbroken gauge symmetry SU(r) x U(l) N ~ r+1 . Then, a part 
of the pseudo-NG modes x turn to NG modes associated with the further symmetry 
breaking SU(N)l+r+ c —> SU(r ) x U(l) N ~ r+1 , so that the total number of zero modes 
is preserved. These new NG modes, called the non-Abelian cloud, spread between the 
separated domain walls [52]. The flavor gauge fields eat the non-Abelian cloud and get 
masses which are proportional to the separation of the domain walls. This is the Higgs 
mechanism in our model. This geometrical understanding of the Higgs mechanism is 
quite similar to D-brane systems in superstring theory. So our domain wall system 
provides a genuine prototype of field theoretical D3-branes. 

One more issue remains to be addressed. That is the question of sign of the gauge 
kinetic term. In our present model, the positivity of the gauge coupling function is 
assured only when positions of walls are properly ordered (see Eq. ( 5.84| )), namely only 
in half of the moduli space. More economical models such as extensions of the three- 
flavor model of Ch. [4] may not have such moduli and, therefore, the effective gauge 
coupling may be always positive. The purpose of the next chapter is to demonstrate, 
that this is indeed the case. 
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6 

Generalized three-flavor model 


In the previous chapter we have considered a non-Abelian extension of the four-flavor 
model of Ch.[4j We have found, among other things, that the effective four-dimensional 
gauge coupling e 4 of Eq. (5.84) depends on moduli of the background domain wall 
solution and that in some part of the moduli space e\ is negative. This points to the 
potential instability of the kinetic term of localized gauge fields. As this is presumably 
not a desired feature, in this chapter we consider a new model with a strictly positive 
effective gauge coupling, which is a function of parameters of the model alone. From the 
last subsection of Ch.[4]we know that the three-flavor model has exactly this property in 
the case of Abelian gauge symmetry. Therefore, in ths chapter we are going to present 
its non-Abelian extension, which we correspondingly call the generalized three-flavor 
model. 

Apart from positivity of the effective gauge coupling, the generalized three-flavor 
model shares all other properties with the chiral model of the previous chapter with 
additional interesting modifications. In consequence, the structure of this chapter will 
be very similar to the previous one. In the first section we will introduce the model, the 
domain wall solution and establish the main claim of this chapter. The second section 
is devoted to derivation of the effective Lagrangian. Compared to the chiral model, the 
full non-linear interaction between moduli fields in generalized three-flavor model turns 
out to be much more complicated. As we were unable to find a close formula, we will 
discuss the effective Lagrangian using the approximation, where we take the thickness 
of the domain wall to be very small. Interestingly, the effective Lagrangian of the chiral 
model is contained within the effective Lagrangian of the generalized three-flavor model 
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6. GENERALIZED THREE-FLAVOR MODEL 



SU(N) C 

U( l)i 

U( 1)2 

SU(N) l 

SU(N)r 

U(1)a 

mass 

Hi 

□ 

1 

0 

□ 

1 

1 

mljv 

H 2 

□ 

1 

-1 

1 

□ 

-1 

0 

Hi 

1 

0 

1 

1 

1 

0 

0 

E 

adj 0 1 

0 

0 

1 

1 

0 

0 

a 

1 

0 

0 

1 

1 

0 

0 


Table 6.1: Quantum numbers of fields of the ungauged generalized three-flavor model. 
The notation is explained in Tab. 


5.1 


as its zero order approximation. We will conclude this chapter with the third section, 
where we discuss other models, leading to positive effective gauge couplings. 

Most of this chapter, content and notation alike, is adopted from the work m- 


6.1 Localization of matter and gauge fields 


In this section we present a model allowing the domain wall solution with unbroken 
non-Abelian global symmetry. As the second step, we introduce non-Abelian gauge 
fields for the unbroken symmetry and consider field-dependent gauge coupling similar 


to Eq. (4.54) to localize non-Abelian gauge fields. We will then demonstrate that four¬ 


dimensional effective gauge coupling is determined solely from parameters of the model 
and that it is strictly positive. 


6.1.1 Lagrangian with global symmetry and domain wall solutions 

Let us consider a five-dimensional SU(N) C x U( l)i x U( 1)2 gauge theory and N scalar 
fields H\ (Ho) in the fundamental representation with the degenerate mass m (— m), 
together with a singlet scalar field H 3 . The global symmetry is SU(N)l x SU(N)r x 
17(1)^. In addition, we introduce adjoint and singlet scalars E and a associated with 
the gauge group SU(N) C x C/(l)i and U( 1 ) 2 , respectively. We summarize charge as¬ 
signments of matter fields in Tab. |6.1[ 
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6.1 Localization of matter and gauge fields 


We assume the following Lagrangian 

£ = ~ Ti ( Gmn G mk ) - T F MN F M!i + Itt(d«e) ! + H (dua f 

+ Tr \DmH\\ 2 + TV I (D M - iA M )H 2 \ 2 + \{d M + iA M )H 3 1 2 - V , (6.1) 

V = TV |(£ - ml N )Hi\ 2 + TV |(S - al N )H 2 \ 2 + \aH 3 \ 2 

+ ^TVjciljv - H\H\ - H 2 HlY + \<?{c 2 + Tr(H 2 H}) - |iT 3 | 2 ) 2 . (6.2) 

The U(N) C = SU(N) C x U(l)\ gauge coupling and gauge fields are denoted by g and 
an N x IV matrix Wm with M = 0,1, 2, 3,4. The 17(1)2 gauge coupling and gauge fields 
are denoted by e and Am- Covariant derivatives and field strengths are defined by 


DmH \ )2 = c)mH\ 2 + \WmH\ >2 , Dm E = 9 mS + i [Wm, E] , (6.3) 

and Gmn = 9mWn - 9nWm + i [Wm, Wn], F M n = 9mA n - 9nA M - 
The global symmetry and the local gauge symmetry act on the fields as 

H\ —»• e ia U c HiU L , H 2 -> e-^ a+ ^U c H 2 U R , H 3 ^e^H 3 , (6.4) 

E —> U c YiU% , it —)■ a , (6-5) 


where t/ L G SC7(iV) L , E7* G SU(N) R , e ia G U( 1) A and U c £ U(N) C , e i/3 G t/(l) 2 


Let us note that the Lagrangian (6.2) can be embedded into a five-dimensional 
N = 1 supersymmetric gauge theory with 8 supercharges. This fact allows us to obtain 
BPS domain wall solution which we consider in the next subsection. 

Without loss of generality, we can assume the mass parameter m to be positive. We 
also assume c\ > 0 and c 2 > 0. Then, there exist N + 1 discrete vacua, where scalar 
fields develop the following vacuum expectation values (VEVs) 


Hi = 


V^( 


( ljV-r 


0 , 


E = rri 


H 2 = -Jci 

1 N-r 


®N-r 


H 3 = \Jc 2 + rci, (6.6) 


0, 


(7 = 0 , 


(6.7) 


where r = 0,1, ... , N. The local gauge symmetry is completely broken and only a 
subgroup of the global symmetry remains in these vacua. The breaking patterns in the 
r = 0 and r = N vacua are 


U(N) C x SU(N) l x SU(N) r x U( 1) 2 x U{l) A 

-> SU{N) l+c x SU(N) r x 17(1 )a+c j 

0—th vacuum 

-> SU(N) r+c x SU(N)l x U(l) A -c- 

N—th vacuum 
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Let us consider domain wall solutions connecting V-th (0-th) vacuum at left (right) 
infinity y = —oo (y = oo). The coincident domain wall solution preserve the diagonal 
subgroup as the largest globaQ symmetry SU(N)l+r, so the Nambu-Goldstone modes 
associated with this breaking take values in the coset 


SU(N) l x SU(N) r x U(1) A 
SU(N) l+r 


( 6 . 8 ) 


To find a corresponding domain wall solution we assume that fields depend only on 
extra-dimensional coordinate y and that all gauge fields, except W y and A y , vanish. 
The energy density can be “completed to a square” as follows 


£ = -^Tr 

r 


A/ s - 9 - ( Cl ljv - H x h\ - H 2 h\ 


+ Tr \D y Hi + (E - ml]s[)HiY 


+ ^ (dy* - e2 ( c 2 + Tr(H 2 H f 2 - |R 3 | 2 ))) + Tr | D y H 2 + (S - {a + iA y )l N )H 2 y 

+ \d y H% + (cr + iAy)H%\~ + c 2 d y <7 


+ cLTr 


ciE - H] H\ (E - ml N ) - H 2 hI(Y, - al N ) 


(6.9) 


Thus, we obtain the Bogomol’nyi bound for the total energy (per unit volume) 

oo oo 

E = / dy £ > T = / dy [ci3 y Tr(E) + c 2 d y a] = Nmc \, (6.10) 


—oo —oo 

where T is the tension of the domain wall. This bound is saturated when the following 
BPS equations are satisfied 


dyHi + (E + i Wy - mljv).ffi = 0 , 

(6.11) 

d y H 2 + + iWy — (c + iAy) H 2 = 0 , 

(6.12) 

d y Hs + (<j + \Ay)H^ = 0, 

(6.13) 

A/E = ^ 2 (ci1jv - HiH\ - H 2 Hpj , 

(6.14) 

d y o = e 2 (c 2 + Tt(H 2 hI) - |R 3 | 2 ) • 

(6.15) 

1 Local gauge symmetry SU(N) C contains global symmetry as a constant gauge transformation, 
which is displayed in the above symmetry breaking pattern of r-th and TV-th vacuua. However, Nambu- 
Goldstone modes only come from the genuine global symmetry which is not locally gauged. Hence we 
do not count the SU(N) C transformations to preserve the color-flavor-locked vacua ( 6 . 61 . 
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6.1 Localization of matter and gauge fields 


As usual, we employ the moduli-matrix formalism by introducing the moduli fields 
S(y) € GL(N, C) and ip(y) € C 


1 


T, + lWy = S dyS, cr + lAy =-dyip. 


The hypermultiplet part (6.11)-(6.13) can be solved by 


H x = z my S~ x Hl , 
Ho = e&S^H? 
H 3 = e~^Hl, 


2 ! 


(6.16) 


(6.17) 

(6.18) 
(6.19) 


with complex constant N x N matrices and a complex constant H 3 , which 


describe moduli of the solution. The vector multiplet part (6.14)-(6.15) turns into the 


master equations for gauge-invariant Hermitian fields fl = SS^ and ry = Re(?/>). 

dyidyttn- 1 ) = ^g 2 (cil N - {e 2my H^H^ + , (6.20) 

= e 2 [c 2 + e r >Tr(H^Q- 1 ) - e~ y |tf 3 °| 2 ) . (6.21) 

Moduli matrices related by the following ^-transformations give identical physical 
fields 

(S, Hi Hi Hi) {VS, ^ + v, VHl VH°e~Z v , H°el v ), (6.22) 

where V € GL{N, C) and v G C. The equivalence class quotiented by this V- 
transformation defines the moduli space of domain walls. We can use this freedom 
to choose the form of the moduli matrices 


Hi — \fc\\ N , Hi — \fc2 ■ 


(6.23) 


Let us also decompose HS as 


Hi = ^/hle^t/l, (6.24) 

where 0 is a Hermitian N x N matrix and U is a unitary N x N matrix. With this 


choice, the master equations ( 6 . 20 ) and ( 6 . 21 ) become 


dy(dy£in x ) = y g 2 (lN - floll , 

\d 2 g = e 2 (c 2 + cie r? Tr(e 29 i H” 1 ) - e _,? c 2 ) , 


(6.25) 

(6.26) 
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6. GENERALIZED THREE-FLAVOR MODEL 


where fio = e 2my 1-N + e 2 ^e ?? . 

No analytic solution of this system of the differential equations is known in general. 
However, one can study essential features of solutions, if one takes the strong gauge 
coupling limit g 2 ,e 2 —> oo. Eqs. (6.25) and (6.26) reduce to a system of algebraic 
equations in this limit: 


Q = e 2my l N + e 2 V , 

C2 = e~ v C 2 — cie ?? Tr(e 2?i O _1 ). 


(6.27) 

(6.28) 


It turns out that the effective theory describing massless excitations localized on the 
background solution of the equations of this system precisely coincides (at least in the 
lowest order of approximation) with the one obtained from (6.25) and ( |6.26 ) (see the 
detailed discussion in Ref. MS)- It is therefore sufficient just to study solutions of 


(6.27) and (6.28). 


If we use the fact that the moduli cj) can be diagonalized by a unitary matrix P 


4> = mP i diag(yi,..., y N )P , 


(6.29) 


we can recast Eq. (6.28) into a polynomial equation of the order N + 1 for x := e 71 


N n 

, d ^ 1 

x = 1 H- x 


— p 2 m(y-yi) 


C2 “ 1 + eiX 
1=1 


e ?: = e 


(6.30) 


If this equation is solved, one can supply its solution into Eq.(6.27) to obtain H 


In the simplest case, where all walls are coincident cj) = myolisr, we can solve 


equation (6.30) explicitly (eo = e 2m ( y y °)) to find 




= tt — ( e o - 1 + VU “ e 0 ) 2 + 4(1 + IVci 
2eo v 


n = (e 2my + e 2mw e' , )lj V . 


(6.31) 

(6.32) 
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6.1 Localization of matter and gauge fields 




y y 

Figure 6.1: Profiles of e~ v and a in the coincident case. The parameters of the plot are 
given above the picture. Positions of all domain walls are centered at the origin. 


Physical fields can be expressed in terms of 17 and a as 


~ yfc{_ - — - 

yi + e~ 2m yy-y Q ) +r > 

(6.33) 

H '2 \Jci --- , 

yj I-)- e 2rn(y-y 0 )-r] 

(6.34) 


(6.35) 

£ = —d y In 17 , 

(6.36) 

a = d y r \, 

(6.37) 

Wy = Ay = 0 , 

(6.38) 


where we fixed the gauge such that S = 17 1 / 2 and Im(^) = 0. 


This solution is not invariant under the symmetry transformations (6.4) in general. 


In the case of U = ljy, however, the fieds (6.33)-(6.38) do not change under the action 
of the diagonal global symmetry SU(N)l+r+ c - We show the y-dependence of e~ v and 
of a = d y rj/ 2 for the coincident case in Fig. |6.1| 

For a more general case, such as non-coincident walls, the dependence of e _J? and 


a on y is more complicated. Furthermore, the equation (6.30) cannot be solved in 
closed form in general, except for first few values of N. Thus, one has to use numerical 


techniques. In Fig. 6.2 we present an example of five non-coincident walls. 
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6. GENERALIZED THREE-FLAVOR MODEL 




Figure 6.2: Profiles of e~ v and a in the non-coincident case. The parameters of the plot 
are given above the picture. Positions of domain walls are y± = —10, y 2 = 4, = 8, 7/4 = 12 

and 7/5 = 16. 



SU(N) C 

U( 1)1 

U( 1)2 

SU{N) V 

U(1)a 

mass 

Ri 

□ 

1 

0 

□ 

1 

ml tv 

H 2 

□ 

1 

-1 

□ 

-1 

0 

h 3 

1 

0 

1 

1 

0 

0 

E 

adj © 1 

0 

0 

1 

0 

0 

a 

1 

0 

0 

1 

0 

0 


Table 6.2: Quantum numbers of the gauged generalized three-flavor model. The notation 
is explained in Tab. 


5.1 


6.1.2 Localization of non-Abelian gauge fields 

In order to obtain massless gauge fields localized on the domain wall, we need to 
introduce new gauge symmetry, which is not broken in the bulk. As we have seen 
in the previous subsections, the coincident domain wall solutions ( |6. 33 >-(6.38) do not 
break a large part of the global symmetry. Let us gauge SU(N)l+r = SU(N)y and 
denote new gauge fields as Vm- Then the fields H\ and H 2 are in the bi-fundamental 


representation of SU(N) C x SU(N)y and the covariant derivatives (6.3) are modified 
to 

DmH\ 2 = 9mH\,2 + iWjvf-ffi,2 - Mi^Vm ■ (6.39) 


Quantum numbers of the gauged model are summarized in Tab. 6.2 
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6.1 Localization of matter and gauge fields 


We introduce the field-dependent gauge coupling for Vm as 

= Act , 




(6.40) 


where we assume that A is a positive constant A > 0. If we denote the field strength 
for the new gauge fields as Gmn , the Lagrangian for the gauged model is given by 


L = L -TV 

2 sV) 


GmnG 


MN 


(6.41) 


where L is the same as in (6.1), except that the covariant derivatives (6.3) are replaced 


by (6.39). It is not hard to see that the solution (6.33)-(6.38) in the ungauged model 


is equally valid in the gauged model. If we write down the equation of motion for the 
new gauge fields Vm we have 

d M G MN = J N , (6.42) 

where Jm stands for the current of Vm- Since the solution preserves SU(N)y, the 


current vanishes for the domain wall solution (6.33)-(6.38), and Vm = 0 is a valid 


solution to the equation of motion (6.42). Then the other equations of motion of the 


gauged model reduce to those of the ungauged model because of Vm = 0. Therefore, 


we see that (6.33)-(6.38), in addition to the condition Vm = 0, solves the whole set of 


equations of motion of the gauged model. 


6.1.3 Positivity of the position-dependent gauge coupling 


We wish to show that the field-dependent gauge coupling (6.40) assures the positive 


definiteness of the position-dependent gauge coupling for any configurations of the 
domain wall. Since we do not need the effective theory in full, we will derive the rest of 
the effective Lagrangian in the next section. For the moment, it is sufficient to know, 
that the field-dependent gauge coupling l/g 2 (a) is given by its value in the background 
solution 

1 1 

= XdyT] = —A d y In x , (6.43) 


5V) 


domain wall 


g 2 {y) 


where x = e v > 0 is the solution to (6.30|). Differentiating (6.30) we find 


'V ^ 


N 


C2 (1 + eixf 
1=1 


2 m H— d v x 
x 


(6.44) 
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6. GENERALIZED THREE-FLAVOR MODEL 


This leads to the formula 


1 


Aci 


N 

E 


1 + — V 

Go / 


AT 


(6.45) 


2y 2 (y) C2 ^(1 + ejx) 2 / V c 2 ^(l + e;x) 2 , 

which is indeed positive in the whole range of y-coordinate. 

We obtain the effective gauge coupling in 3 + 1-dimensional world volume by inte¬ 
grating l/y 2 (y) over the extra-dimensional coordinate y 


- 2 =x 

9a 


dydyrj = A[y(oo) - y(-oo)]. 


The asymptotic values of rj are found from (6.31) as 
r/(oo) =0, 


rj(— oo) = — In ( 1 + N — 
c 2 


(6.46) 


(6.47) 


The easiest way to see these asymptotic values is to note that in (6.32) we can take the 


limits of (6.30) to obtain 


x = 1 + -V 

Go < 


N 


c 2 f 1 + eix 

Z— 1 




if y 
if y -> 


->• oo , 
—oo , 


(6.48) 


since y is finite at both infinities 0 

Thus, the effective gauge coupling is given as 

1 

9a 


—2 = A In ( 1 + N— 


(6.49) 


It is interesting to observe that the effective gauge coupling is proportional to the width 
of the domain wall ln(l + Nci/c 2 ), which is not a modulus, but is fixed by parameters 


of the theory. This feature is in sharp contrast to that of Eq. (5.84), where the effective 
gauge coupling constant is proportional to the separation of walls in each sector, which 
is a modulus undetermined by the theory. We have now confirmed the stability of the 
gauge kinetic term (by choosing the parameters of the theory as A, m, c\,c 2 > 0). 


Eq. (6.46) shows that the effective four-dimensional gauge coupling constant g 4 is 
determined only by the boundary conditions at infinity. This can be interpreted as a 


kind of topological charge of the r/-kink, whose profile is shown in Fig. 6.3 


Note that we have considered only the BPS solutions so far, and confirmed the 
stability of the model. The anti-BPS solutions are also stable, since they can are 
obtained by the parity transformation y —> —y. 


We can just look at (6.19l with Hg = ^[ 02 . and recall ( 6.61 to obtain the same result. 
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6.2 Effective Lagrangian 



Figure 6.3: Profile of 77 -kink is shown in the left panel for the coincident case. In the 
right panel, plots of Tr(£) (green dashed curve), Tr(£) — a (red solid curve), and a (blue 
dotted curve) are shown. 



6.2 Effective Lagrangian 

In this section we calculate the effective Lagrangian on the coincident background do¬ 
main wall solution. We first consider the simple case, where only U moduli is promoted 
to a field on the domain wall’s world-volume, to discuss possible connections with phe¬ 
nomenology. Then we attempt to tackle the general case, which we, however, will be 
able to do only in certain approximation. 

6.2.1 Preliminaries 

In the spirit of the low-energy approximation, let us pick up only two derivative terms 
as a lowest order approximation for the effective Lagrangian0 

=Tr\D fl H 1 \ 2 + Tr\(D^-iA lx )H 2 \ 2 + \(d ll + iA fl )H 3 \ 2 , (6.50) 

with covariant derivatives given as 

D^H h2 = bplht + = d p H ]/2 - iH l/2 V ll + i W, t H h2 . (6.51) 

Here we have singled out covariant derivatives D tl containing fields associated with 
the gauged flavor symmetry SU(N)l + r. In the strong gauge coupling limit there are 

3 We do not explicitly write down the kinetic term for localized gauge fields until the end of this 
section, for simplicity. 
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6. GENERALIZED THREE-FLAVOR MODEL 


two constrains among the Higgs fields: 


H X H\ + H 2 h\ = c, l,y . 
|R 3 | 2 -Tr(// 2 R 2 t ) = c 2 . 


(6.52) 

(6.53) 


First we eliminate W,, through corresponding EoM 


W/I ~ ~ 2 ci L 


H a D,,Hl - DuH a H\ + AA.JIJl!, =: + -AJlAll , 


Cl 


(6.54) 


where we employ Einstein summation convention for the index a = 1,2. Plugging this 


back into (6.50), we obtain 


£ (2) = Tx{D^H a b^Hl) - cTr(lF /t kF"') + iA^R^R* - d^R-jR* - TrtRoD^R] 

- D At R 2 R 2 t )) - 2A /i Tr(fU"R 2 R 2 t ) + A^(\H 3 \ 2 + Tr(R 2 R 2 ) - ^-Tr[(R 2 R 2 ) 2 ]^ 
+ 5 M R 3 t ^ t R 3 . (6.55) 


Before eliminating remaining auxiliary fields An we first rewrite some terms using the 


constraint (6.52) 


Tr[(R 2 R 2 T ) 2 ] = c 1 Tr(R 2 R 2 t ) - Tr(R 2 R 2 T R 1 R 1 T ), 

-2Tr(lE^R 2 R 2 t ) = iTr(R 2 R /( R 2 - D^H 2 H ] 2 ) + -e a& TV[(R a R M Rt 

^1 

- D,H a Hl)H b Hl\ . 

Using these identities A M can be expressed as 

i R 3 9 M Rj - 5 M R 3 R 3 t + ^s ab Tx[(H a b^Hl - b^H a Hl)H b H ] h \ 


(6.56) 


(6.57) 


An = — - 


2 |R 3 | 2 + ^■Tt(r 2 r 2 r 1 r}) 

Next we use the following identity 


(6.58) 


Tr 


b ll Hab' J Hl - aWnW* 


= —Tr 
2ci 


VnJiab^Jil 


ab 


(6.59) 


where !K ab := H} L H b . Let us stress that covariant derivatives act on H\H b as if they 
were in the adjoint representation: 

^(RtRfc) = R(RtR fe ) + W.HlHb - iHlHhV, = dn(H\H h ) + i Wn, R]rJ . (6.60) 
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6.2 Effective Lagrangian 


Nevertheless, to make this fact explicit, we use for covariant derivatives of composite 
fields different notation: D IL . In terms of composite fields l K ab we can rewritte A^ as 


1 H^h\ - d j JhH\ + ±Tr[ k\ 2 V^ 12 - V^\ 2 ) Jfra] 

2 |^ 3 | 2 + ^Tr[lKl 2 ^ 12 ] 
where we used the fact that 

e ab Tr[(H a D^Hl - b^H a Hl)H b H\\ = Tt[ft{ 2 0yK 12 - V^\ 2 )^ 12 \ 
Putting all pieces together we arrive at the effective Lagrangian 


^eff 


d y 


1 

2c i 


Tr(D M 5f afe D^y +d^H 3 d»Hl - A^[\H 3 \ 


+ ^Tr(^ 2 dC 12 )) 


(6.61) 


(6.62) 


(6.63) 


6.2.2 A simple case 

Let us calculate the effective Lagrangian on the background of the coincident domain 
walls 


H\ = ^e^LT 1 / 2 , (6.64) 

H 2 = ^Te^e^O- 1 / 2 ^ ) ( 6 . 65 ) 

Ho, = v/c^e _T?/2 , (6.66) 


with fl and r] given as in Eqs. (6.31)-(6.32). In this subsection we consider, for simplicity, 
the case where positions of all walls are fixed at y = yo and only relative phases and 
center of mass position are allowed to fluctuate. Thus the moduli are promoted to 
world-volume fields as 


= (/> = y 0 (x»)l N . 


Composites ‘K ab and their covariant derivatives are given as 

*-v 


“Kn — ci 


34i 2 — ci 


e 0 e 


1 + e 0 e v 
el / 2 e ~^ 2 


1 N , 

U\ 


1 + e 0 e v 

df 22 = ci—-— ljv ■ 

1 + e 0 e 


(6.67) 


( 6 . 68 ) 

(6.69) 

(6.70) 
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6. GENERALIZED THREE-FLAVOR MODEL 


v crd fl y 0 l N , 

= —^~ e ~ ri sinh(m(y - y 0 ) - ?7/2)[/ t 9 At ?/o + 

©^22 = ^e^adf.yolN, 


Cl 




2 cosh(m(y - y 0 ) - rj/2) 


(6.71) 
, (6.72) 
(6.73) 


where 


a 


m 


C2 


- = i + ^Mi + eoe 


,-V\2 


Nc i 


-l 


Putting this into Eq. (6.61) we obtain 


(6.74) 


A = 1(7 

M 2iVm 




(6.75) 


and applying it into (6.63) we have 


^eff — 


_C2_ 

Nm 


oo 

dy x _ < ^^ e~ T? Tr[D At 17 t T M I/] + mc 2 J dye^ad^yod^yo 


+ 


C2 


N 2 m? 


d y 




1 — a/m 


e _7? Tr [C/D^tyd] Tr . 


(6.76) 


Let us now compute integration factors in the above expression. This is most easily 
done by substituting x = e~ v and using the identity: 


a_ _ ygr^- 1 )) 

m x + {x-l )(l-^;(x-l)) 

leading to the result 


^eff 


+ 


Cl 
2m 
Nmc\ 


(a + l)Tr[D /i ?7 t I) M 17] + — Tr[UV^]Tr[UV^] 
duyod^yo , 


where 



C2 

Nc i 


C2 

Nc i 



In 1 + 



Let us reinterpret our result in terms of pion field, which we define as 


1 

U 


Ti^fc = i 





(6.77) 


(6.78) 


(6.79) 


(6.80) 
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6.2 Effective Lagrangian 


Also let us define 


jd^ :=iTr . 


(6.81) 


Setting decay constants f n and /„ to have canonical kinetic terms, we can recast 


Eq. (6.78) as (including the kinetic term for gauge fields) 


1 


Nrric\ 


£ e ff = Tr () + -d^qd^rj + d^yod^yo - ^Tr 


with 


u = 


1 

2.94 




ci (a + 1) 
2m 


h — 


Cl 

2m 


(6.82) 


(6.83) 


A new feature of (6.82) compared with Eq. (5.35) is that now f w of the adjoint field is 
larger than of the singlet field by the factor \Ja + 1. 


6.2.3 General case I: formulating the problem 

Let us consider the case where both U and </> are promoted to full four-dimensional 
fields. To obtain the effective Lagrangian, we need to repeat the same procedure as 
in the previous subsection, but here the covariant derivatives acting on functions of 


matrices require more care (see e.g. subsection 5.3.2) and causes difficulty when deriving 
the closed form for the effective Lagrangian. However, we have a convenient parameter 
to expand the effective Lagrangian, the ratio c±/c 2 whose logarithm has a physical 


meaning as the width of the domain wall (see Eq. (6.49)). 

The background solution with promoted moduli fields is in the form 


where 


Hi = ^e my n- 1 / 2 , 

H 2 = , 

H3 = s/c2e~ n/2 , 


Q- 1 = 


e~ v e~ 2r f > 


1 N + e 2rny e -2<t> e -r l 


Cl, 


e ~v = 1 + r^Tr 
C2 


L N 


In + e 2m2/ e -2< ^e -?7 


(6.84) 

(6.85) 

( 6 . 86 ) 


(6.87) 

( 6 . 88 ) 


103 

















6. GENERALIZED THREE-FLAVOR MODEL 


Composites “K a b are given as 


!Kn — ci 

= Cl 


LJV 


1 N + e -^m.y e 2<t> eV ’ 
Q m y Q — ^Q—rt/^ 


1 N + e 2my e~ 2 ^e~ v 


U\ 


94 22 — C\U 


L JV 


liv + e 2m ye- 2( t>e-v 


t/t. 


Let us now list some partial results. The derivatives of = r](y,x tJj ) are 


i \dyV = —Tr 


e 2my e -2<t> 


c 2 L (ljv + e 2rny e 2 ^e~ v Y. 

1 


1 + —Tr 


Q 2my e -2<j) 


c 2 l(l N + e 2my e 2 ^e- y ) \ 


Uri = 


du,4> 


(6.89) 

(6.90) 

(6.91) 

(6.92) 

(6.93) 


2 ^ ' 4c 2 V m/ (cosh 2 (y) M 

where we abbreviated y := (my — t]/2)1n — 4>- Covariant derivatives of composites are 


given as 
'Dfj.Jtu = —ci 


,-2 y 


°° 1 / i \ (k+1) 

< 6 - 94 > 


(l w + e -2»)2 




4 cosh(y) 


ci 

2 


2 cosh(y) 11 
(fc+i) 


1 / 1 ... \ 


V^ 22 = UV^ni-y)^ + VyU ^ U* + [7 ^ , 

±N + e y lM + e zy 


where we have used the identity 


1 1 

Vyf(y) = -^f'(y)d»v - E • 


fc =0 


(6.95) 


(6.96) 


(6.97) 


Using above results and recalling (6.61) we may now compute last two terms ap¬ 


pearing in the effective Lagrangian (6.63) 

„2 


1 


duHz&'Ht = —e v (1 - — Y Tr - 

A 3 16c 2 V m) Lcosh 2 (y) M 




Tr 


1 


cosfr(y) 


d^cj) 


(6.98) 


- V4^(|R 3 | 2 + ±Tr(M\ 2 K l2 )) 


= —e^ (1 - — ) Tr 
16c 2 


m/ Lcosh (y) 


Tr 


cosh 2 (y) 


-V y U^U 


, (6.99) 
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6.2 Effective Lagrangian 


Remaining terms can be organized as follows: 

1 


-Tr 

2ci 




ab 


— ^V(y) “f T'mixiy) 7t/(y)» 


where 


^(y) = jTrjl), 


- 7 -ttt + R At tanh(y)D M tanh(y) }> , 


cosh(y) cosh(y) 


( 6 . 100 ) 


( 6 . 101 ) 


7mix(v) = jTr<| U^U 


7u(y) = ^Tr\ 1 + V^U 

cosh (y) 


' LCU- 1 


|_cosh(y)’ cosh(y)J 


+ 


tanh(y), tanh(y) 


|_cosh(y) 




( 6 . 102 ) 

e-« ) 
cosh(y)J 

(6.103) 


Using the identity 

Trl T> lx f(y)T) li f(y) = ^d^rjd^rjTr 


(/'(y)) 2 ] +Tr [©^(/'(y)) 2 ]^ 


n=2 


+E^ Tr 


(n) 


- 2/^(y)/(y)l (6.104) 


we can write down T^(y) as 


ci r 


^V(y) = Tr~r“ e + t T i ' I — 

4 1-cr/m ^ 4 Lcosh z (y) 


OO 1 


2 ^ n! 
n=3 


1 \ (n) 1 


LV cosh(y)/ cosh(y) 


+ tan. 


(n) 


tanh(y)J 

(6.105) 


Additional identity 


00 (_i \ n / 

[/(^/(y^E^^H® 


n=2 


(n) 


- 2 /^(y)/(y)l, (6-106) 


allow us to expand 7 mtx (y) as 


j (y)- 

A n=2 


E HF Tr { (/, ®f t,£ r 1 M [( 


1 \W 1 


LV cosh(y)/ cosh(y) 

(n) 


+ tan' 


tanh(y) > . (6.107) 
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And finally using 


n= 1 


we can express the last piece as 

1 


(6.108) 


Mv) = jTrja^t/tOT 


cosh 2 (y) 


, EE (~i) n , / e y \( n ) Q-v I 

+ 0E[/ t t/V^-^£l([/ t T /i t/) f——— ) ——— L (6.109) 

A | n! ^ \cosh(y)/ cosh(y)J 


Let us now put all our partial results (6.98)-(6.99) and (6.105)-(6.109) together to 
clearly see the work ahead of us: 


£ e fr = / d y{ TV 


ci 


-4 cosh 2 (y) 




- —e v ( 1 - —)Tr 


16c 2 


m 


£ n<t> 


TV 




+ / dy < Tr 


'D u tfV fi u] + e 57 ( 


-4 cosh 2 (y) 


16c 2 


m 


"IE 2 ;/ 

n =3 __ 


1 \ 00 1 


cosh 2 (y)J Lcosh 2 (y)- 

T) fl U^U I Tr I” L/ 
cosh 2 (y)-I Lcosh 2 (y)J I 

(n) 


cosh(y )) cosh(y) 


+ tan. 


oo 00 

+IE // 2 / 


=y \ ( n ) e -y 


ci 


E 

n =2 


(-i) r 


n! 


—oo 
oo 


dyTr|[/tT / ,17^- 1 (T^)[( 


+ (tan’ 


cosh(y) J cosh(-y) / 
1 \(") 1 


tanh(y) 


A 


cosh(y)/ cosh(y) 


(n) 


tanh(y) 


( 6 . 110 ) 


6.2.4 General case II: method 

In this subsection we would like to solve the following general expression 

OO 

I dyTr[/(y)M 


( 6 . 111 ) 


where M is some matrix, independent of y and /(y) is some function. We proceed as 
follows: first we make a substitution x = my — y/2. Since dx = m^l — cr/nijdy we 


106 
















































6.2 Effective Lagrangian 


obtain: 


dy Tr 


f(y)M 


dx 


= / — Tr 
J rn 

— OO 


f(xl N - 4>)M 


f dx aim m 
+ / . ' . Tr 


ml — cr/m 

While the first part can be solved easily by the diagonalization trick 


f(xl N - 4>)M 


dx m 
— Tr 
m 


r 

f{xl N — 4>)M = / — /(x)Tr(M). 

J ./ m 


( 6 . 112 ) 


we found that the second integral in the above equation is very hard to solve in general. 
To make progress, let us adopt the approximation c\/c 2 1. Since ln(l + Nc\/c 2 ) is 

proportional to the width of the domain wall, we call this approximation the thin wall 
approximation^] Thus, we now only aim to solve 


ci 

4c 2 


OO 



— OO 


dx [ ljv 

— Tr -- 

m |_ cosir (xl at 


</0 


Tr 


f(x 1 


N 


4>)M 


(6.113) 


To evaluate an integral over expression involving two traces we use the following identity 


dx Tr 


/(xlv-0)M Tr g{xl N -4>)N = F{d x )Tr 




Tr 


-x<P N 


, (6.114) 


x=0 


where F(x) is given by the convolution 

OO 

F{a — f3) := J dx/(x — a)g(x — /3). 

— OO 


(6.115) 


Applying this procedure, we conclude this section with the identity 

on 

ro 


dyTr\f(y)M] = —Tr [M] + ~^—F 1 (d x )Tr 
L J m 4 c 2 m. 


^M 


Tr 


x<t> 


x=0 


+ 0^(ci/c 2 ) 2 ) , 
(6.116) 


where 

OO OO 

F° = / d yf(y), F\x):= f dy { {y) -(6.117) 

J J cosh (y — x) 

— OO — OO 

4 To be precise ln(l + Nc\/c 2 ) represents the inner width of the wall. Therefore the word “thin” is 
not really justified as the “outer skin” of the domain wall, proportional to 1/m, can be arbitrary large. 
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6. GENERALIZED THREE-FLAVOR MODEL 


6.2.5 General case III: calculations 

Let us denote 


7u ■= J d y 7 v (y ) = 7° v + 7 l v + ciO 

— OO 

wher^] (x := xIat — </>) 


(6.118) 


ci ^ (_!)" 


:=^- V 

4 m n\ 

71—0 


00 

/dx( 


e x \ (n) e -» 


J U ■ = 


„2 .°°. 
— E 


("l) r 


cosh x / cosh x 

Tv 


■TV 


D /1 C/+17<C5(l7 t D'*17)l , (6.119) 


I 6 C 0 W, n! 

71—0 


dx TV 


cosir (x)- 


X Tvi ©,,[/+!/£$([/+2W)( 


(") e" £ 


\cosh(x)y cosh(x)J 


( 6 . 120 ) 


First expression can be easily integrated (and the result is the same as the corresponding 
expression of the four-flavor model): 




where 


FuM 



cosh(£^) 

cosh(y) cosh(y — L^) 


1 j / 1 +tanh(T^) \ 
tanh(£ ? i) \1 — tanh(£^) J 


( 6 . 121 ) 


( 6 . 122 ) 


Using the technique developed in the previous subsection, the next order in c\/c 2 is 
given as 


7b = 


16 c 2 ?b 


-TV 


UFb(d x , £*) (t/t OT)e^' 


TV 




(6.123) 


x=0 


with 


Fb(x,L (j> )= / d y 


cosh(£ 0 ) 


cosh 2 (y — x) cosh(y) cosh(y — L^) 


(6.124) 


5 We included the pure kinetic term for U into the summation as n = 0 term in contrast to (6.1091. 
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6.2 Effective Lagrangian 


Closed form for this integral exists, but it is complicated. The value of (6.124) can be 
better seen in its Taylor expansion. 




sinh 3 (£^,) sinh 3 (£ 0 ) ’ 


d x Fu( 0 ; £<j>) — 


sinh(3£ fl i) 3 — 4 cosh(£^,) 


+ 


3sinlr(£^) sinlr (T^) 


— (3 + 4 — 


d 2 x Fb( 0,^) = ~ 


3 sinh( 2 £ 0 ) cosh(2£^,) 

— AL; 


sinh 4 ^^) ’ 


2 sinh^T^) sinhAL^) sinh & (£^) 


(6.125) 

(6.126) 
(6.127) 


Let us proceed to the mixed term 7 m i X defined as 


^ mix f d y‘Tmix(y) — T^riix + ^ mix T ClO^(ci/c2) ^ , (6.128) 


where 


__ OO 

<to = _c L ^tir r 

mix 2 m^ n! J 

n =2 J 


d.T 


1 \W 1 


+ ^tanh x^j tanh x 


cosh x J cosh x 
X TV 


T 1 . =_ V 


(-l) r 


802 m n\ 

n =2 


We find the result in the form: 


uH 

dxTr\u^‘DML™- 1 ('Dt 1 (j)) \( , , A , . . A . 

\ ^ 0 V VlAcosh(x)/ cosh(x') 


(6.129) 


1 \ M 1 


+ ^tanh(x)^ tanh(V) 


TV 


LV 


-cosir (V)- 


T° ■ = — TV 
m%x 2 m 


t/t V^UF^L^m) 


(6.130) 


(6.131) 


where 


F Lx(£<t>) = d V^~ 


1 - 


cosh(y)cosh(y — L^) 


— tanh(y) tanh(y — L^) 


cosh(£^,) — 1 /1 + tanh(T^)^ 


£0sinh(£^) \1 — tanh(£0)/ 


And 


•T 1 . = _1—Tr 

° mix o ^ _ 
oc^vn 




Tr 


3 - x<t> 


(6.132) 


(6.133) 


x=0 
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6. GENERALIZED THREE-FLAVOR MODEL 


with 


oc 

,(%,£(/>) = J dy 


cosh(£^,) — 1 


1 


cosh (y) cosh (y - L#) cosh 2 (y - x ) 


(6.134) 


Again, the result of this integral can be obtain in the closed form. First few terms of 
the Taylor expansion reads 

2 sinh( 2 £^) — 2 L ( f > 


F mix^^4>) = 
dx F mix(.®i A</>) = 


L<p sinh( 2 T^) + 2 sinh(£^) 

4 sinh( 3 T 0 ) + 9sinh(T^,) — 12T^ cosh(T^) 


d x F mix (o,^) = - 


3T^ cosh(3£,jj) + 2cosh(2£<^) — cosh(£ 0 ) — 2 ’ 
16 2£^ cosh(2£^) — 3 sinh(2£^) + 4 


sinh(4£^,) + 2 sinh( 3 £ 0 ) — 2sinh(2T^) — 6 sinh(£^) 
Finally, let us calculate kinetic term for <fi field: 


(6.135) 

(6.136) 

(6.137) 


V= / dy7^y) = 7l + 7\ + c l O 


(6.138) 


wherc 0 


oo 1 

= yi / dx 

^ 2 m n\ J 

n = 2 -oo 


oo 


C? A 1 


1 \ (n) 1 / \ H 

+ tanh x tanh x 


( — ) 
V cosh xJ 


cosh x 


x Tr 




scl^n! J 


1 


00 1 


+ ^tanh(x 


00 


tanh(x 


A cosh (x)7 cosh(x) 

1 N 


Tr 


cosir (x)- 


It is easy to check that the result is 


irt = —Tr 
* 2m 




where 


(6.139) 


(6.140) 


(6.141) 


0>(Af> = / dy —y 


1 - 


£2 L cosh(y) cosh(y + £$) 

cosh(£ 0 ) — 1 /1 + tanh(£^)\ 

£|sinh(£ 0 ) H VI - tanh (L$) J 


— tanh(y) tanh(y + £^) 


(6.142) 


“We included the pure kinetic term for <j> into the summation as n = 2 term in contrast to (6.1051. 
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6.2 Effective Lagrangian 


and 


with 


r 2 

—Tr 


8 c 2 TO 


■D^Fl(d x ,L^V^) 




Tr 


OO 

j-i i / r x f i cosh(£ 0 ) - 1 
F ^x,L^)= / dy 


3 -a;</> 


i /^2 


x=0 


cosh (y) cosh (y + L$) cosh 2 (y - x) 

— OO 

where the first few coefficients in the Taylor expansion are given as 

pi rn r \ = 2 sinh( 2 ^) - 2£ 0 
^ ^ sinh( 2 £^) + 2 sinh(£ 0 ) ’ 

r) F 1 (0 L ) - _ 4 sinh(3£0) + 9 sinh^) - 12£^ cosh(£ 0 ) 

1 ^ ^ 3£^ cosh(3£ fl i) + 2cosh(2£^) — cosh^^) — 2 ’ 

16 2£^ cosh( 2 £ 0 ) — 3 sinh(2£^) + 4 


^( 0 ,£*) = 


L 2 ^ sinh(4£^,) + 2sinh(3£^,) — 2sinh(2£^) — 6 sinh(£ fl j) 


(6.143) 

(6.144) 

(6.145) 

(6.146) 

(6.147) 


There are still two terms in (6.110) which remain to be evaluated. These are 

1 


r 2 

C 1 


^ ' 16c 2 


a 


dye 71 1-Tr 

m 


— OO 
OO 




16 c 2 


— / dye^fl — — ^Tr 
co / V mJ 


cosh 2 (y) 
1 


7>u<t> 


Tr 


1 


-cosh 2 (y) 


mJ Lcosh 2 (y) 




TV 


1 


cosh 2 (y) 


V^cj) 


-V^U 


(6.148) 


(6.149) 


In our approximation, these terms are already of the order ciO(ci/c 2 ), therefore we 
can set e~ v ~ 1 and cr « 0. Thus 
2 

T d = -T^F(d x )Tr\e x * t D^]Tr\e-^'D^] + Cl C>((ci/c 2 ) 2 ) 

^ ±bC2m L J L J_ n V / 

„2 


x=0 


% = — ^—F(d x )Tr e x *'D u U*U Tr e~ x ^V^U 
16 c 2 m L J 


+ ciO 


i=0 


(6.150) 

(6.151) 


where 


F(x):= / dy 


xcosh(x) — sinh(x) 
cosh 2 (y) cosh 2 (y - x) "sinh(3x) - 3 sinh(x) 


= 16 


(6.152) 


6.2.6 General case IV: conclusion 

To conclude this section, let us summarize our findings. The low-energy effective La¬ 
grangian up to the second order in ci/c 2 is given as 

£ eff = + 7^ + T« + Til + T; + 7' v + c 1 o((c 1 /c 2 ) 2 ) , (6.153) 
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6. GENERALIZED THREE-FLAVOR MODEL 


where 


£$ = —Tr 
ff 2m 


+ 


cosh(£^) — 1 ta / l + t.nh(^) \ 
^ sinh(£^) V1 - tanh(T^) / V 

L 0 sinh(L^J V 1 — tanhfL^j / 


+ '»„!/*[/l„f ; + taD ^ 'l (t/tB-C/) 


tanh(£^,) \ 1 — tanh(£^,) 


(6.154) 


Interestingly, this is the same effective Lagrangian we obtained previously (5.46). The 
rest of terms of the first order in c\/c 2 are given by 


T U) _ 

J u ~ 


T (i) . 

J mix 


16 c 2 ?b 


r 2 

C 1 


8c2m 


-Tr^C/t UFu^L^V^e^ 
tfV^UF mix {d x ,Lt){Wcj))e x * 


Tr 


3 - x<t> 


x=0 


-Tr 


^—Tr 

^ 8 c-2m 


D^)e^' 


Tr 


Tr 

~X(f) 


x=0 


x=0 


7'j. = - — — F(d x )Tr e x *D u $ 

* 16 c 2 m v ’ Y 


and 


7'u = —^F(d x )TV e^D^U 
16c2?n 1 


d y 


Tr 

Tr 




x=0 

e-^V^U 


x=0 


Fu(x,£(j)) = 

mix (*^j ^(j)) = ^11 

—oo 

CO 

F^x.L^) = J dy 


cosh(£0) 


— OO 
OO 


cosh 2 (y — x) cosh(y) cosh(y — L^) 


cosh(£^) — 1 


— OO 
OO 


F(x) = dy 


cosh(y) cosh(y — L^) cosh 2 (y — x) 
cosh(£<^) — 1 

£ 2 cosh(y) cosh(y + L^) cosh 2 (y — x) 

1 

cosh 2 (y) cosh 2 (y — x) 


(6.155) 

(6.156) 

(6.157) 

(6.158) 

(6.159) 

(6.160) 
(6.161) 

(6.162) 

(6.163) 


All the above integrals can be obtained in closed forms. 


The formula (6.153) illustrates the complexity of the interactions between moduli 
fields <f> and U in the general case. Let us offer some explanation for this complexity. 
The structure, which is new in (|6.153[) compared to (|6.154[), is of the general form 


F(d x )Tr[e^M(y)}Tr[e-^N(y)} 


x=0 


(6.164) 
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6.3 More general models of stable position-dependent coupling 


where M(y ) and N(y ) are some matrix-valued functions, containing either derivatives of 
c l i or derivatives of U. If we assume that </> = mP~ 1 diag(yi,..., yN)P is diagonalizable, 
we can rewrite the above as 
N 

Y - y j ))( PM (y)P~ 1 UPN(y)P~ 1 ) jj . (6.165) 

i,j= 1 

This form suggests that interaction (here represented by function F ) depends on the 
relative size of fluctuation of each pairs of walls. Indeed, notice that if two walls have 
the same position yi = yj i ^ j (meaning that the expectation values of the fluctuations 
are the same), the above form reduces to 


F(0)Tr(M)Tr(IV) 


(6.166) 


which is in a sense trivial, since we already encountered this kind of terms in (5.46) 


Thus, the new kind of complexity in our result (6.153) can be understood as a mani¬ 


festation of the fact, that the interaction does not depend only on various moments of 


the fluctuation as in (6.154), but also on their relative size. 


6.3 More general models of stable position-dependent cou¬ 
pling 


In this section, we wish to show that there are more models with the stable position- 
dependent coupling. We will illustrate the point by extending the model to include 
more fields and more gauge symmetry. 

The position-dependent gauge coupling comes from the cubic coupling between a 
singlet scalar field and field strengths of non-Abelian gauge fields in 4 + 1 dimensions 


£ cubic — C ((T,;) Tl ' 


GmnG 


MN 


(6.167) 


where the function a(cTj) of singlet scalar fields cr^ should be linear, if it is to be embed¬ 
dable into a supersymmetric gauge theory in 4 + 1 dimensions 


a ( cr i) = G (6.168) 

i 

where 7 * are constant coefficients. 
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U{ 1)1 

U( 1)2 

U( 1) 3 

mass 

Hi 

1 

0 

0 

m 4 

h 2 

1 

-1 

0 

m 2 

h 3 

0 

1 

1 

m 3 

Hi 

0 

0 

1 

m 4 

<?1 

0 

0 

0 

0 

02 

0 

0 

0 

0 

c 3 

0 

0 

0 

0 


Table 6.3: Quantum numbers of the U( l)i x 17(1)2 x 17(1)3 extended model. 


Usually each domain wall has one complex moduli: a position and a phase. For 
example, both the massive CP 2 model and the massive CP 1 x CP 1 model have two 
free domain walls, corresponding to the two complex moduli. Although, two free do¬ 
main walls can produce the desired profile of position-dependent gauge coupling by an 


appropriate choice of parameters in Eq. (6.168), they provide the undesired modulus 
for the width of the profile. To avoid this problem, we are led to consider models with 
a single complex moduli. The simplest one of such models is the three-flavor model in 
Ref. [31] (or see Ch. [4]), where three scalar fields are constrained by the two Abelian 
gauge symmetry U( 1) x U( 1). The generalized three-flavor model of this chapter is an 
extension of this model to non-Abelian gauge group: U{ 1) X 17(1) —> U(N ) x 17(1). 
The next simplest possibility is to consider four scalars constrained by three Abelian 
gauge symmetry 17(1) x 17(1) x 17(1), which we call extended models. We give quantum 


numbers of fields of a typical extended model in Tab. 6.3 In the limit of strong gauge 
couplings, the gauge theory becomes a nonlinear sigma model whose target space is 
given by an intersection of three conditions as 


(c 2 x cp 1 ) n (h 2 x cp 1 ) n (c 2 x cp 1 ) ~ cp 1 . 


Here H 2 denotes hyperbolic complex plane. 

The vacuum condition is given by 

|Pi| 2 + |P 2 | 2 = ci, -|R 2 | 2 + |P 3 | 2 = C 2 , |P 3 | 2 +177 4 | 2 = c 3 , 

Pi(cti - mi) = 0, H 2 (ai — <t 2 - m 2 ) = 0, 

H 3 (cr 2 + a 3 - m 3 ) = 0, H±(g 3 - m 4 ) = 0, 


(6.169) 
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(6.171) 

(6.172) 










6.3 More general models of stable position-dependent coupling 



\Hi\ 

\h 2 \ 

m 

\H4\ 

71 

72 

73 

(1> 

0 

Vci 

V C l+2 

V C 3-l-2 

”*2+3-4 

”*3-4 

777-4 

( 2 ) 

x/cl 

0 

y[C2 

\/ c 3-2 

7711 

m 3 _4 

777-4 

(3) 

\/ c l+2 

V~ c 2 

0 

y/C3 

mi 

***1-2 

7714 

<4> 

y/Cl+2-3 

V C 3-2 

\/C3 

0 

m\ 

***1-2 

***2-1+3 


Table 6.4: VEVs of candidate vacua: We use abbreviations like C 3 - 1-2 = C 3 — C\ — c 2 . 


where Cj is the Fayet-Iliopoulos parameter of U(\)i and m a is the mass for Hj. All 


possible solutions to these equations are shown in Tab. 6.4 There are four solutions 
but only two of them are valid solutions for any choice of real parameters of c+ When 
we choose ci > 0 , C 2 > 0 and C 3 > ci + C 2 , we are left with the vacua ( 1 ) and ( 2 ) in 
Tab. IOI 

The moduli matrix formalism is powerful enough to give generic solutions of the 
BPS equations of this nonlinear sigma model. Especially, we are interested in the kink 
profiles of <ti, 02 and 03 . They can be expressed as derivatives of real functions rji 


(J i = \dyVii (* = 1,2,3), 


(6.173) 


where the real functions r/i are determined by the following algebraic conditions 


e -rn+ 2 m 1 y e -r) 1 +r)2+2m 2 y-2a _ ^ 
_ e -yi+V2+2m 2 y-2a _|_ e ~ri 2 -y 3 + 2 m 3 y _ ^ 
e -V2~V3+2m 3 y _|_ e -y3+2m 4 y _ 


(6.174) 

(6.175) 

(6.176) 


with a real constant a. The parameter a is (the real part of) the unique modulus of the 
solution, corresponding to the position of the domain wall. As we expected, we find 
only a single modulus. The width of the domain wall is not a modulus, but fixed by 
the theory. 

Since we want a configuration that m —> 0 at both spacial infinities, we can choose 
m 3 = 7714 = 0- Then cr 2 = 73 —y 0 at y = ± 00 . Several numerical solutions are 


displayed in Fig. 6.4 From Fig. 6.4, we clearly see that <72 > 0 and <73 < 0 for all the 
values of y. As a position-dependent gauge coupling, we can choose a two-parameter 
family 


£ C ubic = - ( 72^2 - 73<73) Tr 


GmnG 


MN 


(6.177) 
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c-j — 1, c 2 — 1 j c 3 — 3 j — *1, m 2 34 — 0 


Ci = 10 3 , c 2 = 1, c 3 = 10 3 +1, = 1, m 2 3 4 = 0 




Figure 6.4: The kink profiles of (red solid line), 07 (green dashed line) and <73 (blue 
dotted line) for two different sets of model parameters. 


where 72,3 can be any non-negative real numbers. This class of models can be easily 
made to localize non-Abelian gauge fields and minimally interacting matter fields by 
extending two of the U( 1) factor groups to (possibly different) U(N ) gauge groups 
with N scalar fields in the fundamental representations, similarly to our model in 
this chapter. A new interesting feature of the extended model (and its non-Abelian 
extensions) is that two possible profile of singlet fields 07 , <72 can provide a different 
profile for different non-Abelian gauge groups such as SU( 3), SU(2) and 17(1) with 
their associated matter fields. 
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Conclusion and discussion 


In this thesis we have investigated the possibility that the brane-world scenario might 
be realized in the field-theoretical framework by a topological soliton. As a model of the 
brane we used a simplest soliton possible: a BPS domain wall. A natural consequence 
of this choice was that the laboratory, where our models were constructed, was (4+1)- 
dimensional field theories embeddable into N = 1 SUSY gauge theories with eight 
supercharges. Within this class of theories we showed that localization of non-Abelian 
gauge fields, together with minimally interacting matter fields, is quite natural, once 
it is realized that the gauge symmetry, which is to be localized on the domain wall, 
should not be broken in the bulk. 

One of the interesting features of our models of Ch. [5] and Ch. [6] is the geometric 
realization of the Higgs mechanism. Since we used as a background a coincident do¬ 
main wall, representing the most symmetric point in the moduli space of solutions, we 
managed to localize a large gauge group. It is quite possible that quantum or other 
corrections will result in departing from this coincident point and some gauge fields 
becomes massive, eating corresponding NG bosons in the process. Thus, we see an 
analogy to the usual Higgs mechanism based on purely geometrical basis. Amongst the 
possible future investigations, it would be worth-while to study non-coincident solution 
to further clarify this geometrical Higgs mechanism. 

We have also noticed that a part of moduli fields in the low-energy effective La- 
grangian resembles the pion in QCD. Similar attempts have been quite successful using 
D-branes [62]. We believe that our methods can provide more insight in various aspects 
of low-energy hadron physics. 
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Moreover, at the end of the chapter [6] we explored possible generalizations of our 
stabilization mechanism by including more fields and more gauge symmetries and found 
a class of more generic models with an added flexibility for model building with different 
localization profile for different gauge groups. 

To build realistic models of brane-world with our scenario of localized gauge fields 
and matter fields, we should address several questions. Perhaps the most important 
question is to obtain (massless) matter fields in representations like the fundamental 
rather than the adjoint of localized gauge fields. One immediate possibility is to use 
the localization mechanism of fermions in a kink background. It has been found that 
zero modes of such fermions are localized in such a way to give automatically chiral 
fermions [3]. 

Secondly, we should devise a way to give small masses to our matter fields in order 
to do phenomenology. Since some of our matter fields are the NG modes of a broken 
global symmetry, we need to consider an explicit breaking of such global symmetry. 

Thirdly, another question is to study possibility of supersymmetric model of gauge 
field localization. We need to settle the issue of possible new moduli in that case, which 
we discussed in Ref. m- Moreover, we should examine the mechanism of supersymme¬ 
try breaking. For example, we can consider a model similar to the one we introduced 
in Ch. [5j but where different halves of supercharges are preserved at each sector (BPS 
and anti-BPS walls). In such a model the SUSY is completely broken in the system as 
a whole. It has been proposed that the coexistence of BPS and anti-BPS walls gives 
the supersymmetry breaking in a controlled manner [63]. In our present case, BPS and 
anti-BPS sectors interact only weakly. If we choose flavor gauge fields for each sector 
separately, we have only higher derivative interactions induced by massive modes. If we 
choose the diagonal subgroup of each sector as flavor gauge group, we have a more inter¬ 
esting possibility of the massless gauge field, propagating as a messenger field between 
two sectors. 

Another interesting possibility for a model building is to localize gauge fields of 
different gauge groups with different profiles. This situation is often proposed in recent 
brane-world phenomenology, for instance in Ref. (61] . 

Finally, let us examine similarities and differences of our domain walls compared 
to D-branes. The most interesting similarity of our domain wall with D-branes is the 
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realization of geometrical Higgs mechanism, where massless gauge fields in the coinci¬ 
dent wall become massive as walls separate. On the other hand, there are differences as 
well. D-branes in string theory are defined by the Dirichlet condition for fundamental 
string attached to it, but no such condition is visible in our domain walls. Domain walls 
similar to ours have been constructed in Ref. {65j, where a probe magnetic charge is 
placed in the bulk Higgs phase. The magnetic flux from the probe magnetic monopole 
is carried by a vortex which can end on the domain wall. They observed that this wall¬ 
string junction configuration resembles an open string ending on D-branes in string 
theory. However, this phenomenon is in theories with one dimension less, namely the 
world volume of domain walls has only (2 + l)-dimensions (fundamental theory is in 
(3 + l)-dimensions). In our model based on a theory in (4 + l)-dimensions, monopole 
(codimension three) is a string-like soliton and might possibly be a candidate of some¬ 
thing similar to the fundamental string. It is worth pursuing a possibility of composite 
solitons consisting of domain walls and other solitons, such as monopoles, in order to 
clarify more similarities of solitons with D-branes. Another interesting issue of the 
effective action on D-branes is the non-Abelian generalization of the Dirac-Born-Infeld 
(DBI) action. There have been studies to obtain first few corrections to Yang-Mills 
action in string theories. Since we are at present interested in up to quadratic terms in 
derivatives, we have found the ordinary quadratic action of Yang-Mills fields together 
with the action of moduli fields interacting minimally with the Yang-Mills fields in ad¬ 
dition to the nonlinear interactions among themselves. Our result is trivially consistent 
with the quadratic approximation of DBI action, but does not give us informations on 
non-Abelian generalization of DBI action. In order to shed light on that issue, we need 
to compute higher derivative corrections to our effective Lagrangian. 

We plan to consider all these issues in future studies. 
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